Problème aux limites de Laplace — système de coordonnées toroïdales — p 
Solution de l'équation de Laplace en coordonnées toroïdales (n,Ÿ,@)_ avec conditions aux limites 
Dirichlet, Neumann ou Robin par la méthode de séparation des variables 
Le système toroïdal (n,0,@) (3D) est défini par le changement de variable : 


_ aSinh(n)Cos(p) 
E Cosh(n) -— Cos(0) 
= ASinh(m)Sin@) | wee 0< n < +00 
Cosh(ņ)-— Cos(0) 
"E a Sin(0) 
Cosh(n) — Cos(0) 





—r<0<r 0<p<2r 











où a fixe la dimension du rayon de l'anneau du tore 
n désigne la distance à l'intérieur de la section du tore 


Ô l'angle polaire dans la section du tore 
@ l'angle sur l'anneau du tore 


Typiquement nous recherchons la solution du problème aux limites intérieur en coordonnées 
sphéroiïdales allongées sur un espace bornée dans les limites : 

ne [0,+c0] ; Qe [7x ] > PE [0,27] soit 1€ [mn] ; de [9] > PE [pP] 
A noter que l'intérieur d'un tore en coordonnées cartésiennes correspond à un problème extérieur 


de Dirichlet (en coordonnées toroïdales), soit 7 S Ino] ; 9era] ; pelo27 l dont voici 
l'illustration avec l'anneau fin figurant une valeur radiale n tendant vers l'infini, et un rayon de plus 


en plus faible. 





Tandis que le problème sur le domaine ” S lon] ; Selrx] ; pelo2r] correspond à un 
problème extérieur au tore en coordonnées cartésiennes. 
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Le Jacobien en coordonnées sphéroïdales allongées et la métrique sont : 
x _, Cosh(n)Cos(p)(Cosh(n) - Cos(0)) - Sinh(n) Cos(o) _ : Cos(p}(1-Cosh(n)Cos(0)) 
ôn (Cosh(n) -Cos (0) (Cosh(n)- Cos(0)° 
ox Sinh(n)Cos(o)Sin(0) 
= 7 
00 (Cosh(n)-Cos(0)) 
x 7 Sinh(n)Sin(@) 
op (Cosh(n)- Cos(0)) 


ðy _ Cosh(n)Sin(gXCosh(n)- Cos(8))— Sinh(n) Sin(o) Le Sin(@Y1- Cosh(n)Cos(0)) 



































ôn (Cosh(n)-Cos(®) Y (Cosh(n)- Cos(0)) 
Oy sr Sinh(n)Sin(@)Sin(0) y ai Sinh(n)Cos(o) Oz Er Sinh(n)Sin(0) 
00 (Cosh(m)- Cos(0)) ôp  (Cosh(n)-Cos(6)) ôn (Cosh(n)- Cos(0)Ÿ 
& _, Cos(OYCosh(n) - Cos(0))- Sin’ (0) _,_Cos(@)Cosh(n) 1 êz 
00 (Cosh(n) - Cos(0)) (Cosh(n)-Cos(0)) Ôp 
Ox ôx x 
ôn 00 ôp 
yeda E | a(Cosh(n)- Cos(0 X° Sinh(n) L a Sinh(n) 
ôn 09 ôp (Cosh(n) - Cos(0)) (Cosh(n) - Cos(0)) 
Oz Oz Oz 
ôn 09 ôp 
5? = 4 -~ (dn? +d0? + Sink? (ndo?) 


(Cosh(n)— Cos(0)) 


Formule d'inversion du système de coordonnées toroïdales 


Il s'agit de passer du système cartésien vers le système toroïidal. Pour cela plaçons-nous dans un 
plan méridien p, z, p étant la distance à l'axe de révolution z. Il est facile de montrer que le système 
toroïdal est défini de la manière suivante par ce que l'on appelle une fonction de révolution, 


fonction de la variable complexe 1 * ið. 


z+ip=ia Coran 227°) 


Pour obtenir la formule désirée il suffit d'inverser cette fonction, soit : 


n +i0 = 2Cotanh (PTE | = 2AreCotanh| PIE) 
a a 


p=ÿx+y 
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Iso-surface du système de coordonnées toroïdales 


Voici deux iso-surfaces n=ctse du système toroïdal, comme nous l'avons évoqué les anneaux 
intérieurs ont une valeur n plus grande. 





Tale < 1 


Voici trois iso-surfaces Ô=n/5, 37/4, -n/3 du système toroïdal, ce qui permet de modéliser des 
unions et des intersections de sphères : 











La jonction entre les sphères (iso surface 9= Cste) se réalise à la limite n>. 
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sphères : 


Voici des iso-surfaces Ÿ=n/5, 37/4, ce qui permet également de modéliser des intersections de 











Là encore la jonction entre les sphères (iso surface Ÿ= Cste) se réalise à la limite n°. 
Valeurs particulières du système toroïdale 
La valeur n=0 correspond aux valeurs cartésiennes : 


x=0 y=0 





= CAm) avec -mn<0<T 
1— Cos(0) 
Cela décrit finalement l'axe z des coordonnées cartésiennes. 


La limite n> correspond aux valeurs cartésiennes : 
x=aCos(ọ) y=aSin(g) z=0 avec 
Soit une cercle de rayon a dans le plan (x,y). 


0<p<27r 
_ a Sinh(n)Cos(@p) 





Les valeurs Ü=+n correspond aux valeurs cartésiennes : 





= LD A z=0 avec O<ņ<+% 0<p<2x 
1+ Cosh(n) 1+ Cosh(n) 
Ce qui décrit l'intérieur d'un disque de rayon a centré à l'origine dans le plan (x,y). 
La valeurs 9=0 correspond aux valeurs cartésiennes : 
Pia a Sinh(n)Cos(o) 
Cosh(n) -1 





_ aSinh(n)Sin(o) za e 
Cosh(ņ)—1 
Ce qui décrit l'extérieur d'un disque de rayon a centré à l'origine dans le plan (x,y). 





O<mp<+o 0<p<2r 
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Comportement des iso-surfaces =o 


Les iso-surface =, forment des sphères, comme suit : 
EE ad Sink? (n) É i ceo) a? (1-Cos(6,)Cosh(m)) 




















(Cosh(n) - Cos(0,)) Sin(6,) }  Sin°(6,) (Cosh(n) - Cos(6,)) 

2 2 Cos(6,) F a” . 2 .… 12 E 2 
x +y iE : ee) . ua eo (0,)Sinh? (1) + (1- Cos(8,)Cosh(m} ) 
n a 2 TEE LR 
= PE CSS PCR (0) + Cosh" (n) 2Cos(0,)Cosh(n)) SmO) 

2 2 | eD) a? 
>x +y +|z-a = > 
Sin(0,) Sin“ (0,) 


Ces sphères ne sont pas complètes, limitées par le disque de rayon a centré à l'origine dans le plan 
(x,y). La fonction Sin(ĝo) étant croissante pour les valeurs Ÿ=0, €[0,x/2], on voit que les rayons des 
sphères s'amoindrissent successivement. Chaque sphère est alors englobée par la précédente, 
jusqu'à la sphère centrée à l'origine de rayon minimum a pour o= n/2. Comme indiqué sur la 
figure suivante : 











Au delà pour 9 €[7/2,7], les sphères sont successivement plus grande, mais le centre s'éloignant 
vers des valeur négatives z->-c, les sphères ne s'englobent plus, mais de rayon grandissant de 
centre à l'infini négatif et de rayon infini lorsque ĝo -> n (dégénérant en l'extérieur du disque de 
rayon a pour ), comme suit dans la figure suivante : 
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Pour ce qui des angles négatifs, on observe le même comportement en fonction de la croissance de 
la valeur absolu des angles entre [-7/2,0] CE 
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L 1 r 1 
D © © 
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Et entre [-n,n/2] : 


-1 


—0.8-0.60.40.20.0 


Problème aux limites de Laplace — système de coordonnées toroïdales — p 


Comportement des iso-surfaces n=no 


Ces iso-surfaces représentent des anneaux concentriques d'équation cartésienne implicite : 


e E a° Sinh’ (7) + a° Sin? (0) + a° (Cosh(n) — Cos(0)Ù _ _2@ Cosh(n) 














(Cosh(n) - Cos(0)) Cosh(n) - Cos(0) 
yera a” Sinh’ (N) n a Sinh(n) D PE a > x + y’ 
(Cosh(n)-Cos(0)Ÿ ` (Cosh(n)- Cos(O)) (Cosh(n)-Cos(0))  Sinh(n) 


2 2 
yx +y -24 Cosh(n) APE) + 
Sinh(n) Sinh(n) 
>x + y’ +z + a? = gy EN Jx? +y 
Sinh(n) 
Lorsque n=0 l'anneau se réduit à l'axe z et lorsque n>% l'anneau se réduit à un cercle de cente O 


>x +y +z’ +a’ =2aCosh(n) 





2 dre... 12 — A s, 
et de rayon a: * +Y =a" zZz=0_ Ente ces deux extrêmes le rayon de l'anneau ne fait que 
décroître. De ce fait l'espace extérieur à l'iso-surface n=noest n<noet l'espace intérieur par n>no 


Le dessin suivant illustre le comportement de ces iso-surfaces : 


n=+0.5 
— +1 
— j+15 

f=+2 

n=+0.5 


+15 
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Singularité du système de coordonnées toroïdales et régularité du gradient 


Le jacobien du système de coordonnées toroïdales ne présente pas de valeur nulle sur un domaine 
bornée, et s'annule à l'infini. En conséquence, il n'y a pas de problème de régularité du gradient: 


GradtT (,0,9)- (Cosh(ry) — Cos(@)) (2222), (7028), pi (mea, 
a 














ôn 60 Sinh) | 64 


Le Laplacien en coordonnées toroïdales 3D (n,ĝ,ọ) s'écrit sous la forme : 


Ô Sinh(n) ÔT(n,0,@p) N 
ôn | Cosh(n)-Cos(0) ôn 


0 
+ Sinh 
in os 





(Cosh(n)- Cos(0)) 
a Sinh(n) 








AT(7,0,@)= rare), 


Cosh(n) - Cos(0) 00 
A (Cosh(n)- Cos(0)) TN, 0, p) 
a Sinh(n) ôg’ 
Lorsque le problème ne dépend pas de la variable ọ, soit est indépendant de la position sur le 
pourtour de l'anneau, alors le Laplacien en deux dimensions (n,0) s'écrit : 
© | Sinh(n) rer) 
(Cosh(n)- Cos(8)) | 2n \ Cosh(n)—Cos(0) ôn 


a” Sinh(n) ô 
+ Sinh 
inh(n) 50 | 

















AT(n,0)= 





Cosh(n)-Cos(0) 00 


L'équation de Laplace dans le système de coordonnées toroïdales n'est pas simplement séparable, 
soit s'écrivant sous la forme du produit de 3 fonctions des coordonnées indépendantes. Mais en 
revanche elle est R-séparable. Cela signifie que proportionnellement à une fonction commune des 
trois coordonnées, elle devient simplement séparable. Il existe actuellement deux procédés de R- 
Séparabilité pour l'équation de Laplace en coordonnées toroïdales. 


1 Tuer) 
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Forme traditionnelle de la R-Séparabilité en coordonnées toroïdales (R-Séparabilité-1) : 
T(7,0,p) = R(7,0) H(7)O(0)#(p) 
R(n,0) = JCosh(n) — Cos(0) 
H") + Sen À + É E -}# ()=0 
AT(n,0,p) = 0 = 16"(0)+ «,@(0)= 0 
F'(p)+a;F(p)=0 








en posant les cas suivants 


Cas(l) a,=n et &, = m => 





Cosh(n) E é 
Sinh(n) JE? -1 





Posons € = Cosh(n) = = Sinh(n)=Cosh°?(m)-1= VE? -1 Sink? (N) = €? 
n 





a H6 dé dH(é) _ FE 
dn dn dé 


d'H(£) _ dé z(20)- E 2E Er Ft). (e JHO , g dH(E) 


























dn? dndë\ dn DERE dé? dé 
2 d’'H() | dH (8) | 2_ m 2 2) H(&) dH(&) [> 1 n? | 
>(£ -—1) P LÉ n e E?) a er (r ’ D JA 


=> H (n) = AP" (Cosh())+ BO" ,(Coshn)) 


@"(0) + n°@(0) = 0 = O(0) = ACos(n0)+ BSin(n0) 
P'(o)+ mY (p) =0 > Y(p)= ACos(mp)+ B Sin(mçp) 
Cas(2) a,=n et æ, =0 => le système ne dépend plus de p 
T É - 20 = 02 H= AP, ,(Coshtm))+ BO, , (Cosh(n) 

@"(0) + n°@(0) = 0 = @(8)= ACos(n0)+ BSin(n0) 

P'(E)+m P(p)=0= Y(p)= A+Bp 

Cas(3) a, =0— valeur propre nulle (lorsque le système ne dépend plus que de la coordonnée radiale) 


a Cosh(n) 
FO Se 


H"(n) + 








HG) + 2 HG) = 0 





= H(n)= AP ,(Cosh(n))+ BQ ;(Cosh(n))= T(n,0) = \Cosh(n) - Cos(6) (4 P ,(Cosh(n))+ BO , (Coson) 


En réduisant à deux dimensions avec un problème ne dépendant pas de l'angle ọ, il suffit de 
prendre les solutions des deux derniers cas. Cette forme de R-séparabilité est particulièrement 
adaptée aux problèmes dépendant uniquement des variables n, Ÿ, en se réduisant à l'utilisation 
des fonctions de Legendre de première et deuxième espèce. 
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On peut également a, <0 également possible. En gardant un solide de révolution azimutale 
complet (soit une constante de séparation azimutale, en ọ, égale à m entier), il donnerait des 
équations de la forme, 


Cas(4) @,=-T° et æ, = m > 

F Cosi) A Jro =0 
Sinh(n) Sinh (n) 

=> |0"(0)-1°0(0) = 0 & @(0) = ACosh(t0)+ B Sinh(t0) 

P"(p)+ mY (p) =0 < Y(p)= ACos(mp)+ B Sin(mo) 


2 


H" (n) 





mol +T’ 





z = Cosh(n) 


Posons À Sinh(n)= 42° —1 
2 _ 2 _ d 
d Vr-id d vz lad| Vz i z 











dn dz dn? dz 
1 m 
=> R" (n) + Cotanh(mR'(n) +| =+ 7? Rm) = 
(n) + Cotanh(n)R'(n) É T a (n) 


2 





olz? -jeere ier + koso 


1- 7° 
2 " ! il 2 m’ 
ell-z )R (z)—2ZR'(z) - 11° a z 
-zZ 


RM) =A P*" (Cosh(n))+ B O*! (Cosh(n)) 





ko =0=R(2)=4P% (z)+BO* (z) 


——+iT 
2 


= 4@(0) = ACosh(r0)+ B Sinh(rO) 
W(o) = ACos(mp)+ B Sin(mp) 


Ce choix de solutions est adapté à des problèmes aux limites inhomogènes en dimension angulaire 
Ÿ (sur des intersections, union et différence de sphères). On utilise les fonctions coniques de Mehler. 
Dans ce type de problème, notamment lorsque la géométrie implique toutes les valeurs de la 
dimension radiale, les transformations de Mehler-Fock sont particulièrement commodes pour 
construire formellement les solutions sous forme de représentations intégrales, par exemple sur 
des configurations géométriques comme des lentilles sphériques. 
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Nouvelle forme de la R-Séparabilité en coordonnées toroïdales (R-Séparabilité-2) : 


T (1.0.6) = R(n.0) U(n.6,p)= R(n.0) H(MO(6)F(p) _R(n.0) = = 


2 2 2 
AT(1,0,9) = 0 & AU(,0,9) = CE er, RES 2), Up) UE o 
7 p 











& (Sinh (E"MOO YE) + OH MEE) oH (MO) + ZOTO -o 











> (sinop) (Z 5 ae be nn F'(E) _ m 
H(n)  @(@)) 4 Y(E) 
S H"(n) 1 _ (a Je re @"(8) SE 
Hn) (Sinh) 4 E(6) 


(Sinh mY (H" 7) - a, H (n))+ È z a, Ho =0 


AU(n,0,@) = 0 = {6"(0)+ «,@(0) = 0 
P'(p)+a;F(o)=0 





en posant les cas suivants 


Cas(l)a,=n eta;=m => 


(Sinh MY (H"(7)-n°H(n))+ G z me} (7) = 0 





dn (Sinh(m) (Sinh) Y 
dH(E) a AO O a PH), ra E 
Sa N a A a 


(CE D ae -ER r)i- Jon =0 


Posons  £ = Cotanh(n) 2e } (Cosh) — (Sinh) = £?-1 











> 2 

E 
(e egw +2 E ( m? ao =0=> H(n)= AP" (Cotanh(n))+ BO" ,(Cotanh(n)) 
©" (0) + n°@(0) = 0 => O(0) = ACos(n0)+ B Sin(n8) 
P"(p)+ mY (p) =0 => F(p)= ACos(mp)+ BSin(mp) 
Cas(2)a, =0 et a; = m° 











(e ED +26 Am + G nè Jr =0>H(n)=AP_,(Cotanhn))+ BQ  ;(Cotanh(n)) 


@"(0)=0— O(0)= A+ BO 
PE) + n'Y(p)=0= Y(p)= ACos(mp)+ B Sin(mo) 

Cas(3) a, =0 et a, =0 

(Sina HG) + ZR < 0 

=> H(n)= AP (Cotanh(n))+ BQ ;(Cotanh(n)) et @(0)=A4+B6 et W(o)=A4+B@ 


Cette forme de R-séparabilité est particulièrement adaptée aux problèmes ayant une dépendance 
à l'angle azimutal op. 
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Tout comme dans le cas de R-séparabilité-1, un quatrième cas de valeurs des variables de 
séparation peut-être envisagés, notamment lorsque æ <0. En gardant un solide de révolution 
azimutale complet (soit une constante de séparation azimutale, en @, égale à m entier), il 
donnerait des équations de la forme, 

Cas(4) a,=n eta;=-T => 


(Sinh (4m -nH ))+ G + 20 = 0 


= |(@"(0) + n°@(8) = 0 = @(8) = ACos(n0)+ B Sin(n0) 
P'(E)- Tr Y(p)=0—= Y(p)= ACosh(rp)+ B Sinh(rp) 





z dé 1 _(Cosh(m)) -(SinhY _ s _ 
Posons  Ë = Cotanh(n) An Ska} Qu é —1 
dH(Ë) _ {z2 ,\dH(é) CAE) (R yË LGR MELAS] 


(ee -7 20 a JE mo. FR 











> 2 
& 
2_ THN) dH (n) 

= (£ spe dE J 


+2 

a 

(q-e) EED 2 0 | 
+ 


n 
2 2 





jron =0 








+r’ i 
4 4 
l 2 n? E LE +n tn 
Te ATE = Ho) =0> H(§)= AP LEE 9, 6) 
H(n) = AP; | (Cotanh(n)) BOT l (Cotanh(n)) 


=> 4 ©(0) = ACos(n0)+ B Sin(n@) 
W(@) = À Cosh(rp)+ B Sinh(rp) 
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Contraintes et conditions imposées pour les solutions de l'équation de Laplace pour les 
problèmes aux limites dans le système de coordonnées toroïdales. 


D'après les formules de changement les coordonnées toroïdales se comportent à l'inverse des 
coordonnées cartésiennes : 
- a Sinh(n)Cos(@) z aSinh(n)Sin(@) >= a Sin(0) 
Cosh(n) — Cos(0) Cosh(n) - Cos(0) Cosh(n) - Cos(0) 


avec 








O<np<+o -r<0<r 0<p<2r 
n—o—x=aCos(g®) y=aSin(g) z=0— cercle de rayona 
a Sin(0) 

| 1-Cos(0) 
Dès lors la description d'un problème intérieur en coordonnées toroïdales n correspond à un 
problème extérieur au domaine décrit en coordonnées cartésiennes, et inversement. Typiquement 
un problème d'équation de la chaleur à l'intérieur d'un anneau torique représente un problème 
"extérieur" en coordonnées toroïdales, et inversement un problème d'électrostatique extérieur à un 
anneau chargé représente un problème intérieur en coordonnées toroïdales. 





n=0—x=0 y=0 => axe z 


La première contrainte qui s'impose à ces problèmes est d'abord que les solutions soient finies en 
tout point du domaine étudié, paramétré dans le système de coordonnées toroïdales. La deuxième 
contrainte est liée à la forme géométrique du domaine. Si par exemple le domaine est un anneau 
torique complet, il convient que les solutions soient 2r-périodique tant en Ÿ que ọ. 


T(n,0,@) fini,ne comporte aucune singularité, continue et dérivable 
T(n,0,@) périodique en 0 de période 2x soit T(n,0,0)=1(n,0 +7x,œ) 


OLGAT re en 0 de période 2x soit OT(n,0,p) _ OT(n,0 +7,¢) 
op 00 00 


T(n,0,@) périodique en o de période 2x soit T(n,0,@)=T(n,0,0 +27) 


CHU A périodique en o de période 2x soit OT(n,0,p) _ OT(n,0,p+2x) 
0 0p 0p 








La périodicité des solutions dans les angles Ŷ et @, liée à une géométrie d'anneau complet (et non 
de section d'anneau dans les deux angles) imposent nécessairement que les paramètres n et m des 
solutions soit des entiers naturels, tant dans la R-Séparabilité-1 que la la R-Séparabilité-2. De plus 
contrairement à la sphère, il n'y a pas de contraintes de bornes supérieures de m comme pour la 
cas d'une sphère où par exemple m <= n ( cas R-Séparabilité-1) . 
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On peut ajouter d'autres contraintes liées à la répartition des conditions aux limites. On peut 

imaginer par exemple des conditions aux limites indépendantes de @, qui soit paires en È ou bien 

impaires en Ÿ, ou bien des conditions aux limites dépendantes de , qui soit paires en ọ ou bien 

impaires en œ. Ce qui signifie que les solutions exhibent également les mêmes propriétés : 
Conditions aux limites dependant pas de @ 


T(n,8,@)= T(n,0) 
Conditions aux limites paires en 0 Conditions aux limites impaires 


T(7,0)=T(n,-0) T(n,0)=-T(n,-0) 


Conditions aux limites générales paires en 0 


T(7,0, p) = T(n,-0,p) 
Conditions aux limites générales paires en Q 


T(n,0,p)=T(7,0,-œp) 


Conditions aux limites générales impaires en 0 


T(7,0, p) = -T (ņn,—0,p) 
Conditions aux limites générales impaires en @ 


T(n,80,ọ) m= -T (n,0,—ọ) 


Pour la contrainte de finitude de la solution, il convient de choisir une base de fonctions propres qui 
n'induit pas le choix de fonctions divergentes sur le domaine étudiée. Seules les parties radiales 
(dépendance en n) induites dans la résolution R-Séparable de l'équation de Laplace présente un 
risque de divergence. Ce sont des fonctions de Legendre de première et deuxième espèce. Par 
exemple sur un problème bidimensionnelle torique (n,9), dans la R-séparabilité normale, les 
solutions s'écrivent sous la forme d'une série : 


70) = Can Cot] > De (Cosk(n))+ 8,0. 


n=0,+00 





Cost) (as Cos(n9) + B,,Sin(n o») 


Et dans une configuration tri-dimensionnelle du problème aux limites (R-Séparabilité-1) : 
T(N,0,9) = Cosh(n)-— Cos(0) x 


{Z 


n=0,+00 m=0,+00 





[ape P” (Cosh(n))+ B} PO ET) 


et dans la re les solutions restent dépendantes des trois coordonnées : 


T(n.0.0)= [eee , 








x | [ar  (Cotanh(n))+ B",0", (Com) |4 Costa 9) + B*Sin(nO)\A?Cos(mp)+ B:Shtmo)] 
n=0,+00 m=0,+00 MT MT 


En réalité la R-Séparabilité-2 n'introduit pas de nouvelles fonctions propres à partir desquelles 





développer une solution, car il se trouve que les “harmoniques toroïdales" sont liées par les 
formules de Whipple : 





tn") > 
r|» a :) z Sinh) 
2 
(m-n) 
2 | 2 
(mens 2) A anh) 
2 


Ce qui prouve que les deux développements en série sont identiques ! 


5 
P” (Cosh(n))= | Q"_; (Cotanh(n)) 





2 
9": (Cosh(n))= ) Pey (Cotanh(n)) 


Problème aux limites de Laplace — système de coordonnées toroïdales — p 


La divergence des solutions dépend donc du comportement des fonctions de Legendre aux points 














connus de singularité, qui sont ici z=+1 ou -1 ou z=+% : 
| Singularité P, (z) | Q, (z) 
(raie) -y (u+ D) x 
2 fitu +ne-D] fed D) 
-y+ PE Da- ye- 
2“T| u+ -) 
Z= +0 R Valu +1) y 
Var (u +1) f(a : 3) 
Coo 2y (— = Cotan(zu) R 
| Sin(zu) iy É + Hfi _ (u+1) ci D) 2 + Log 127) 2Sin(xu) 
T i 2 2 
z=-1 -y (-u)-y(u+1)-2y fi-en) 
7 Sin(xu) ED (1 T j 2 
ENS ET 
“ii Fr _ a) Qu + l - Costau] - Ge D) 








Soit en simplifiant avec les valeurs limites : 
Singularité P, (z) Q, (z) 
z=l fini + 00 
z = +0 + 00 Jini 


z=-1 fini +œ u entier sinon fini 
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Pour les deux formes de séparabilité, le paramètre z change de comportement, pour 


z = Cotanh(n) = 


la fonction est croissante en n, par contre pour 
décroissante en n. Le choix des solutions finies peut se résumer par le tableau suivant : 


z = Cosh(n) 


Cosh(n) 


Sinh(n) — la fonction est 





Problème intérieur 
coordonnées cartésiennes 


en 


Problème extérieur 
coordonnées cartésiennes 


en 





Domaine cartésien 


Intérieur d'un anneau torique 


Extérieur d'un anneau torique 




















Coordonnées toroïdales Problème extérieur en Problème intérieur en 
coordonnées toroïdales coordonnées toroïdales 

Domaine de la variable n [No +) [0,no] 

Domaine de la variable |[Cosh(no),+°) [1,Cosh(n0)] 

Cosh(n), R-séparabilité 1 

Domaine de la variable |(1,Cotanh(no)] [Cotanh(no),+°e) 

Cotanh(n), R-séparabilité 2 

Fonctions radiales finies R-\Q,:2(Cosh(n)) P,12(Cosh(n)) 

séparabilité 1 

Fonctions radiales finies R-\P,:,(Cotanh(n)) Qn-12(Cotanh(n)) 





séparabilité 2 














Problème aux limites de Laplace — système de coordonnées toroïdales — p 


Problème intérieur de Dirichlet à deux et trois dimensions: soit à rechercher la solution de 
l'équation de Laplace à l'intérieur d'un tore dont les conditions aux limites de Dirichlet dépendent 
des angles 6 et ọ, puis du seul angle 8 : 


AT(n,0,9)=0 (n,0,p)eQ= h > 1,0 € |-z,zļo € [o,2x |} „-„ fini 
T(7,0,p),., = f(6,p) oubien T(,8,p),, =f(8) oubien T(7,6,p),, = f(e) 


Le problème en coordonnées toroïdales se présente donc comme un problème extérieur, du fait 
du paramétrage du système de coordonnées. On rappelle les contraintes de ce problème exposées 
précédemment : 

T(n,0,@) fini,ne comporte aucune singularité, continue et dérivable 





T(n,0,@) périodique en 0 de période 2x soit T(n,0,0)=1(n,0 +7x,@) 


GPO en 0 de période 2x soit TAN,O, p) 2 ÔT(n,0 +x,@p) 
op 00 00 
T(n,0,@p) périodique en ọ de période 2x soit T(n,0,@) =T(n,0,p +27) 
_— périodique en @ de période 2x soit CRE O). are +27) 
p p p 








R-séparabilité 1 : problème intérieur dépendant des variables n et 6 


Cela conduit au choix d'une partie radiale décroissante en n et s'annulant à l'infini, et des parties 
angulaires possédant des valeurs propres n et m entières, soit une solution qui se développe sous la 
forme d'une série : 


T(n,0,@)= JCosh(m) - Cos(0) 





(z Q (Coshm) X Y O” (CA Cast) + B, int), „Cosmo D, Sn) 


n=1l,+00 m=1,+00 2 


Lorsque les conditions aux limites ne dépendent pas de ọ. La solution s'écrit : 


roo- En Ca | 8,0 | (Cosh(n))+ X Q à COTE) 





n=1,+00 
Limitons nous au cas sans dépendance à w. Lorsque 7 conditions aux limites angulaires sont 
paires, alors la solution se développe en série : 


J0,p)=f(@0) et f(-0)= f0) 
reno) = CHEN C0) B, Q (Cosh(m))+ X. 4A „Cos(n9)Q (Gant) 


n=1,+00 





Dans le cas de conditions aux limites impaires, la série prend la forme : 


f(6,@)=f(6) et f(-0)=-f(0) 
T(n,0)= /Cosh(n)- wO E A, Sin(n9)Q (Goshen) 


n=1,+00 





Problème aux limites de Laplace — système de coordonnées toroïdales — p 


Compte tenu des conditions aux limites, nous avons alors la solution générale pour des conditions 
aux limites de Dirichlet paires, f(8)=fo(-8), sous la forme : 


T(n.0),., = a de 2 4,Cos(n9)Q (Cosh) 





0,(0)=1> |©, =27 ©,(0)=Cos(n0) => 
fa (0) 





z fas fa (0) Cos(nẸ) 
Cosh(1,)-Cos(6) 














A = 





À, 





Ie = 





o 


z Jo(0) af Jo(0) 
Bo = fao JCosh(1,)-Cos(@) ni JCosh(1,)-Cos(6) 


B, = = fas Jo (0) Cos(n9) -2ùf 49 LO Cos(n#) 
JCosh(1,)-Cos(6) JCosh(l,)- Cos(0) 


Q ,(Cosh(n)) Q (Cosh(n)) 
g 2 + > B, 2 
20 ,(CoshG,)) „re Q, (Cosh(l,) 
































ĮJCosh(n) -Cos (0) 





T(n,0) = Cos(n39) 


Pour développer plus avant des cas particuliers de conditions aux limites on utilise le résultat des 
représentations intégrales des fonctions de Legendre de degré demi-entier. Sachant que selon 
Lebedev "Special Functions and their applications", 1965, section 7.10 formule 7.10.10, ainsi que 
Louis Robin « Fonctions-Sphériques-de-Legendre-et-Fonctions-Sphéroidales, Tome III », 1969, page 
113, également cité dans un l'article de S.C.Loh, 1959 « The Calculation of Electric Potential and 
the Capacity of a Tore by means of Toroidal functions » formule (6) page 699: 


o 1 (Cosai)= [de Cosin) "1 fgg Cos(ng) 1f Cos(n0) 








V25 (Cosh(n)- Cos(8)}> 242; (Cosh(n)-Cos(6)P 22 (Cosh(n)- Cos(8)}> 
= o, (Cosh(m)=Q  ,(Coshn)) 


pour n = 0,1,2,3,... 


Il vient pour une condition aux limites di simplifiée: 


fo(0)=1 B, = 








fao TO a “2920 (Cost) 


: Cos(nà) 
B, - jao JCosh(1,)-Cos (0) 








=2V20 1 (Coshü, ) 








T(n.0)- 242 Cohen Cos(0) 3 Q ,(Coshn))+ > Cos(n9) Q à (Ca) =] 


n=1,+00 
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Normalement cette expansion en série doit converger uniformément vers 1 car dans ce cas la 


solution est triviale. Pour le montrer utilisons le résultat établi ci-après sur le problème extérieur 

de Dirichlet (problème du potentiel électrostatique extérieur d'un corps au potentiel donné). La 

solution s'écrit (voir plus loin) : 

[Qi (Cosh(1,)) 
2 


2 P ,(Cosh(l,)) 








P ,(Cosh(n))+ 





T(n.0)=V,2V2 VCosh(n)-Cos(0) 
7 Q ı(Cosh(l,)) 


+ : Cos(n9)P (Cosh 
27 La) 7 Ce) 


n—— 
2 





La solution résulte du développement en série sur un système complet (mais non orthogonal) de 
fonction toroïdales, sur un problème aux limites qui prend la valeur Vo en n=l, en portant 
justement l'expression à la position en coordonnées toroïdales n=l, il vient avec une valeur de 
potentiel V,=1: 


ns JCosh(1,)-Cos(6) | 








T 


To (CoshG,)}+ X 0, (Con cast 


72 n=1,+00 2 


La valeur l, pouvant être quelconque on en déduit immédiatement pour toute valeur de n : 


i= Ce Cfg (Coshon)+ X, 0 ,(Coshtn)Coscn9) 
TT Tz I= 


n=1,+0 3 








On peut d'une certaine façon qualifier cette égalité de « partition de l'unité » pour les fonctions 
toroïdales. 


Problème aux limites de Laplace — système de coordonnées toroïdales — p 


Problème intérieur de Dirichlet sur un tore avec une condition aux limites impaire en à : 





Selon des calculs similaires à ceux de l'exemple précédent, la solution générale pour des conditions 
aux limites de Dirichlet impaires, f.(8)=-f,(-8), sous la forme : 


=n TO), = 1O = JE O] Z 4,Sin(n9)Q (cart) 


1 
n=1,+00 2 





©,(0) = Sin(n0) = | 





O, 




















8,= [as fa (0) Sin(n9) -243 fa (0) Sin(n9) 
"7  „Cosh(l,)-Cos(0) ẹṣ JCosh(,)-Cos(8) 








_ VCosh(n)-Cos(0) B, , 
T(n,0)= - pr 2 Con) Sin(n DO, : (Cosh(r)) 


Prenons maintenant une condition aux limites simplifiée dans laquelle la surface est portée à la 
valeur Sin( ©), Il vient : 

















B= 2f 49 Sin(9)Sin(n9) Sin(n 9)Sin(9)= Cos((n —1)9)- Cos((n +1)9) 
o JCosh(I,)-Cos(0) 5 
s Em Cos((n —1)9)- Cos((n+1)9) _ ÿ; (c ni k c _ ) 
ma | JCosh(,)-Cos(8) gs osh(l,) Gi osh(l,) 





Q ,(Cosh(n)) 


n= 





T(n.0) = JE AD z 


P fo. (Cosh(l, ))- Q 1 (Cosh(l, )] Da) Sin(n9) 


== 
2 
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Prenons également une condition aux limites simplifiée dans laquelle la surface est portée à la 
valeur Sin(mŸ), Il vient : 

















B, = 2[ 48 SORA) Sin(n 9)Sin(m 9) = Cos((n - m)8)- Cos((n + m)3) 
o JCosA(I,)-Cos(0) 2 
a , (Cosh(l,))- (Cosh(l | 
f Cos((n—m)9)- Cos((n+m)9) _ ol o> a) 9, i 0) do 
= B,= [49 TEE L 
0 JCosh(1,)-Cos(6) fo, a (Cosh(,))-Q non) si m<n 








JCosh(m)- Cos(6) Qı (Cosh(n)) 
T(7,0)= V2 = >. CM (Cosh(l, ))- 2, (Cosh(1, )] o (Ct) Sin(n9) 


Peu 
2 

Prenons maintenant une condition aux limites simplifiée dans laquelle la surface supérieure est 

portée à la valeur +1 et la surface inférieure à la valeur -1. Il vient : 


n=1,+00 


























f,(@)=1 6ef0,7] _4 Sin((2m +1)9) 
f,(@)=-1 80e a FR 2 (2m +1) 

8 1 4 o Sin((2m+1)9)Sin(n9) 4 1 4 Cosfn-2m-—119)-Cos((n+2m+1)9) 
B, = fas = fas 

A mTio (2m+1)y Cosh, )-Cos(0) x, (2m+1); JCoshG,)- Cos(0) 

442 1 
E r TONN (CoshQ,))- PR (can, ) 

h(L,))- h(l 
Taj, o) - 42E -Cos0) Con (Cesç,)) oo $ ) 2n an Sin(n9) 
HAS 2 un. (2m +1) Q, (Cosh(l,)) 


De même pour le problème électrostatique extérieur la solution s'écrit : 





(o ,(Cosh(l,))- Q , 


|[n-2m-1|-= n+ PTR 
2 2 


(can) P ,(Coshn)) 
5 Ha 








Ta- 4V2./Cosh(n)- Cos(0) 
r m’ n=1,+00 m=0,+00 (2m + 1) P 1 (Cosh(l,) 
n= 


Problème aux limites de Laplace — système de coordonnées toroïdales — p 


Solution général au problème de Dirichlet sur un-tore complet pour des conditions aux limites en 
Ÿ admettant un développement en série de Fourier 


Compte tenu des résultats obtenus précédemment, en supposant que les conditions aux limites 
admettent un développement en série de Fourier : 





VCosh()- Cos(8) Q 1 (Cosh(n)) 











£ (4, Cos(n9)+B, Sin(n9)) Jo(0) 


T(7,0)= 2, Q (Cosh(,) =; [49 JCosh(1,)-Cos(6) 


„= fas Jo (0)Cos(n9) „= fas Jo(O)Sin(nd) 
JCosh(1,)-Cos(6) /Cosh(1,)-Cos(0) 














Si f,(0)= D af Cos(m9)+af Sin(m9) 
m=0 


af = — | dO f,(0) af =— L d0 f,(0)Cos(n9) af =— Li d0 f,(O)Sin(n 9) 
27 ”, RE T 


= 


A (F af Cosimo) + af Sinoma) (Cost) 
























































+= [as E 35 fas (af Cos(m9)+ af Sin(m9))Cos(n 9) 
ood JCosh(l,)- Cos(0) TA JCosh(l,)-Cos(0) 
E Cos(m9)Cos(n9Q) 5 Sin(mÿ9)Cos(n4) 
=Y af [43 f 
| [JCosh(l,)- Cos(0) a jas [JCosh(l,)- Cos(0) 
r jas Sin(m9)Cos(n9) | ja Cos(mÿ)Cos(n9) =2fa3 Cos(m9)Cos(n9Q) 
JCosh(,)- Cos(0) ka j= aa Toil, )-Cos(0) 
sä Do Jas Cos(m9)Cos(n9) D as jas Cos( n— ts 
JCosh(,)- Cos(0) £ 2 JCosh(1,)- Cos(0) 


= À, = >: an 0, à (Cosh(,)) A, = E D (9 (Cosh(l,))+ E (cout) 


m=0,+00 








(F at Cosimo) + a Sim 8) tn R 
d9 m=0 






































pz Í _ S [as (af Cos(m9)+ af Sin(m9)}Sin(n 9) 
A JCosh(1,)-Cos(6) re JCosh(1,)-Cos(6) 
RS Cos(m9)Sin(n8) > Sin(m9)Sin(n3) 
= £ |do 1 
Le | JCosh(1,)-Cos(6) Ža jas JCosh(1,)-Cos(6) 
r fas Cos(m8)Sin(n8) Ar [as Sin(m8)Sin(n8) = 2| 49 Sin(m9)Sin(n9) 


JCosh(1,)-Cos(6) JCosh(1,)-Cos(6) JCosh(1,)-Cos(6) 


= B, =2$ a; sJas Sin(m9)Sin(n9) _ a |as Cos(n - ml9)-Cos((n + m)9) 
JCosh(i,)-Cos(0) te "n Cosh(1,)-Cos(6) 














> B,= FD D in E (Cosh(l,))- Eli (Cosh) 


Problème aux limites de Laplace — système de coordonnées toroïdales — p 
Il vient la solution suivante : 


af = 1 Í d0 f,(0)Cos(n9) af = 2 Í d0 f,(0)Sin(nd) 
T $ T i 


Q ue 
Y afo 1 (Cosh, )}=—2 











m=0, +0 m Q \Cosh(l, I 
Q i (Cosh) 
re,o) = VEP O à. 2 2 ato, (Cont), P (Cosh, Ja TIA Cos(n9)+ 
Q (Cosh) 
B B l 
2 2% (o, (Cosa) Q a (Cosha, > Cosh(,) Sin(n9) 


Extension au problème d'un demi-tore dont la base est portée à la valeur 0 





Soit le problème sur le demi-tore : 

AT(n,0,9)=0 (n,0,p)eQ = h € l, [0 € [0,7 Lo € [,2x ]} 

T(7,0, PB = f(8) T(7,6, Pia, a 
pel0,27] pel0,2x] 

D'après la construction des solutions sur un tore entier avec des conditions aux limites impaires : 

ATGL0,p)=0 (n6,p)eQ=hel,+0l0el-7x,rlee[027]) 7(7,0,0)... fini 


T0 PW g =) f(-0)=-f(6) 
p ] 


pel0 27] 





„-„ fini 

















Q a (Oa) 
5 7(7.6) - JCosh(n) - Cos(8) S 8, 
T 


n=1,+0 'Q (Coma, ) Sin(n2) > T(n,0) =0 
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En raison de l'unicité de la solution du problème de Dirichlet même dans un domaine non connexe 
comme le tore, on a trouvé une solution pour le problème de départ. Si le demi tore est porté à la 
température constante To : 


AT(n.0,p)=0 (0,0,p)EQ =f e|, +0}0 [0x1 [027] 7T(n,0,p) 





. fini 














T C0) =T, T(n,0,p)1ln+-]=0 
pelo,27] pel0.27] 
Ce (Cosh(l,))- 0, (Cosa) Q (Cosh(n)) 
=> T(n,0) = AT, Cohen Cos(0) >. pa 2 Gr) 2 0 Tor) 
x Sin(n9) i 





Pour une fonction limite quelconque en introduisant le développement en série de Fourier : 


af =l ao j,(0)Sintn9) =? [40 FOS) 
T *, TS 


À = fao f,(O)Sin(n 9) 





Q (Cosh) 
T(n,9)= V2 y Cosh(n) Cos(8) afo a (Cosh(,))-Q (Cost, )] 72 Sin(n9) 
m n=1,+00 m=0,+0 LE 72 7 


im Q , (Cosh(1,)) 


n—— 
2 
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Problème intérieur de Dirichlet sur un tore complet et creux: soit à rechercher la solution de 
l'équation de Laplace à l'intérieur d'un tore creux dont les conditions aux limites de Dirichlet sont 
fixés en paroi intérieure et extérieure à des valeur différentes : 


AT(n,0)=0 m0)eQ=hel,.1,l0el-7;,r} 
TO  =7 Te _ =P 


6e[-7,r] oef-r,r] 





On ramène ce problème à un problème homogène de Dirichlet sur la face intérieure, à savoir : 
T(1,0)=T; +T, (7,0) 

AT,(7,0)=0 (n,6)eQ=fhel,,,.) 

TO a, =0 TO a, =B-T 


Pour trouver la solution on peut facilement mettre à profit la formule de partition de l'unité 
appliqué aux valeurs aux limites du domaine : 


a J- 
TTN d Cosh) =O io (Cosap) Y o (Cosh, kostn) 


1 
72 n=1,+00 2 


n=l 











T, -T =(T, -T,)2V2 





EC Ho  (Cosh(l,,))+ 0 (con Dkcstn)] 


1 
72 n=1,+00 72 


Problème aux limites de Laplace — système de coordonnées toroïdales — p 
La solution étant de la forme (conditions aux limites paires) : 


B, Lo  (Cosh(m))-Q (cart) + 


a Ai 
2 2 





T,(1,0) = JCosh(n) - Cos(0) 
+Y ACatn| 9  (Cosh(m))-Q (Can) 


n=1,+00 "z 


T,0,0), = 7-7 


B, Lo (CosA(,,))-0 , (con) + 


2 2 





=> /Cosh(1,,)-Cos(0) 
F2 ACatn| 9 (Cosk(,,))-0 


n=1,+00 


(Cosh(l,, )] 


= (T, -TP\2 TP EO) Q ,(Cosa(,,))+ X Q (Cost )Xostn) 
T n=l, +0 "7 


SL 
2 








1 
2 -3 3 


B, Lo ,(Cosh(l,»))- Q (con) + 


2L aL 
2 2 


1 
26 79 1 (Cosh,)}+ 








<> = (T, = T) 
RY ACotn| 9 a (Cosh,,))-Q ; (cat) di 2 9, 1 (Cosh,2))Cos(n8) 
Q :(Cosh(i,:)) 
B, = (T, T, ) 12 2 
TT 
(e, (Cosh(,,))-Q ı (Cosh(l, )] 
dé 5 | O , (Cosh) 
4=(T, RE 2 z 
T 
Lo (CosA(,,))-Q nn) 





Il vient donc la solution du problème : 
T(n,0)=T, + 


(o, (Cosh(m))-0 , (con) 
2 Q ; (Cosh(l,, 1 





2 2 
+ 
25; 





O :(Cosh(,,))-Q ,(CoshG, )] 


+(T, or) JCosh(n)- Cos(0) 
TT 


n—— n—— 
2 2 


Lo (Cosh(m))-Q ,(CosA(L,,)) 








+ $ Q (Cosh(l,,)) 


n=l+o "73 


Cos(nd) 
Lo  (Cosh(,,))-0 (con) 


n—— n—— 
2 2 
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Soit encore : 








r'(9,0) = 7,243 ACosh@n)= Cos(6) | 79 1 (Coshtn)+ 
m 2 


=1,+00 n 


(e 1 (Cosh(n))- 
(osha, „)) 2 


Q (Cosh(n))Cos(n 9) + 


o (Euh) 


+ 





fo 1 Cosh(l, Rim 





0 (Ent) 


+6, -2) 22 JOoGn CO) 
| | 
o 


a (Cosh(m)-Q (Cosh(l,, )] 








n=l,+o 72 


+Y 0 ARR 
Q 


T(n,0)= 


_1 (Cosh(i,:))- Q i (Cosh(l, )] 


Cos(nd9) 


Q ,(Cosh(l, fo, (Cosh(m))-Q ,(Cosh(l, )] 





ZTO à (Coshm)+ (7-1) 





_ 2424 Cosh(n)- Cos(0) TO 


i (Cosh(n))+ 
T 2 





+ >. Cos(n9) 


n=1,+00 


2 





+(7,-T)0 ;(Cosh(l,,)) 


(o, (Cosh(,,))-Q , (Cosh(l, )] 


(0, (caen) 0, | (cart) 





Lo (Cosk(,))-Q ; (cart) 


Il reste à prouver que cette formule donne une fonction qui ne dépend pas de l'angle 8 ! ce qui 


n'est pas une mince affaire ! 


+ 
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Problème intérieur homogène et inhomogène de Dirichlet sur un tore complet et creux: soit à 
rechercher la solution de l'équation de Laplace à l'intérieur d'un tore creux dont les conditions aux 
limites de Dirichlet sont inhomogènes en paroi extérieure et homogène en paroi intérieure 


AT(n,0)=0 (n,0)eQ= in € Pts . € -zrl 
TD =) T0) 


RTS a 


Compte tenu de la condition homogène de Dirichlet, la solution se développe comme suit : 

P (Coshn)) Q :(Cosh(n)) 

T(n,0) = /Cosh(n)-— Cos(6) 2 2 ( 
>, ı (Cosh(l,2)) O à (CosA(L,,)) 


n=0,+00 








A,Cos(n9)+ B,Sin(n0)) 


Et en introduisant le développement en série de Fourier de la fonction limite : 


af =} fao f,(6)Cos(n8) af =} fdo f,(6)Sin(ns) 
me TI’, 





VCosh(n)- Cos(0) . 


T(n,0)= 2 








(Cosh) Q :(Cosh(1,:)) 





s 2 
Le dy nA (Cosh(l, )) P, (Cosh(l,)) 0, (Cosh) + 


2 








P, (Cosh(n)) ee) 


P ,(Cosh(l,)) O ,(Coshtl,»)) 


P (Coshn)) Q (Cosh) 


is 
w| 








n=0,+00 m=0,+00 








= 

x > > fe 

+ Am (o, (Cosh (L ))+ Fa (Cosh (Z, )] P (Cosh (lp )) Q 3 Cosh (l, )) Cos (n Q) F 
a 


P ı(Cosh(n)) Q à CEA) 


n— 





P (Cosh) © (Cash) 


2 


f h(l )- h(l 72 cm : 9 
+ D 2 af (o, Cos a,)) 9, 1(Cos a) P (Cost) o E J ) 


(CoshG,,)) Q p (Cosh(L,,)) 


j-i 
2 





Seny 











Problème aux limites de Laplace — système de coordonnées toroïdales — p 


Problème intérieur homogène et inhomogène de Dirichlet sur un demi tore complet et creux: soit 
à rechercher la solution de l'équation de Laplace à l'intérieur d'un tore creux dont les conditions 
aux limites de Dirichlet sont inhomogènes en paroi extérieure, homogène en paroi intérieure et 


homogène nulle sur la tranche 





D'après la construction précédente de la solution sur le tore entier creux, lorsque la condition aux 
limites inhomogène est impaire en angle 8, alors la solution est nulle sur la tranche. On peut donc 


utiliser cette solution pour ce problème sur le demi-tore. : 
AT(n,0)=0 (7,0)eQ= in € [4:10 € [0,x |} 























T u” K hO) T oM =0 TOO) = 
af = rfa fOSin(n9) = T(n,0) = Ÿ a LONG), 
P ,(Coshm)) Q (Cosh) 
P, (Cosh) ©. (Con) 
, EE 4/0, (Cent), (Cona) 2er EE: E Sin(n9) 
P (Cosh) 0 (CG) 





Problème aux limites de Laplace — système de coordonnées toroïdales — p 


Problème intérieur homogène de Neumann et inhomogène de Dirichlet sur un tore complet et 
creux: soit à rechercher la solution de l'équation de Laplace à l'intérieur d'un tore creux dont les 
conditions aux limites de Dirichlet sont inhomogènes en paroi extérieure et extérieure à des valeur 
différentes : 


aT0)=0 MDQ- helher) TO), = ET, 0 
de q Ke] Nm 
Compte tenu de la condition homogène de Neumann, la solution se développe comme suit : 


P, (Cosh(n)) Q (Cosh(n)) 


T(,0)= Cosh -Cos Y | = 72 — (4,Cos(n9)+ B,Sin(n0)) 
are ,(Cosh(l,2)) PL  (CosA(I,,)) 
nm" nm" 








Et en introduisant le développement en série de Fourier de la fonction limite, on trouve la solution 
suivante : 


af = A Í d0 f,(8)Cos(n9) af -o Í dO f,(0)Sin(n3) 
A T Tz 





VCosh(n)- Cos(0) 7 





T(n,0)= 














2 
P(Cosh(m)  Q :(Cosh(n)) 
Ê p (Cosh) Zo 1 (Cosh(l,s)) 
X a*Q  (Cosh(,)) CEE tir ai 
m=0,+0 #2 a (Cosh (la )) Q, (Cosh (Li )) 
£ P ,(Cosh(,)) -Z0 ,(Cosh(,)) 
3 TE. a 


P (Cosh(n)) Q : 1 (CoshGn)) 
: ô 





ð 
a (Cosh(1,,)) à (Cosh, ,) 


xl+ D OT ato, (Cost ))+0, (Cort) ea =} 





n=0,+00 m=0,+00 P 1 (CosA(1,,)) Q à 1(CoshG,)) 
2 








Lo 


1 (Cash 
on © osh(l,,)) 


a 7 (Cosh(1,,)) 


P (Cosh(n)) Q à COS CD) 





26 
P \Cosh(l — c h(l 
on al os (l,2)) on © osh( „2)) 





‘+ È afo ,(Cosh(l,))-Q (Cosa) Sin(n9) 
ae m0, Fee rte P (Cosh) O, 1 (Cosain) 

2, 2 

ôn »; 








 (CosA(L,,)) a 1 (CosA(,:)) 





Problème aux limites de Laplace — système de coordonnées toroïdales — p 


Divers exemples tirés de la littérature 


Homogénéisation de l'équation de poisson pour un système de coordonnées orthogonales de 
révolution 


On exprime une solution particulière de l'équation de Poisson avec un terme source constant, dans 
le cas d'un système de coordonnées orthogonale de révolution (autour de l'axe z). Alors la 
coordonnées z ne dépend pas de l'angle azimutal @, tout comme le rayon p : 


Intérieur Poisson AT, (a,B.p)=Q Tala, B.p)=Ta(a, B) 





Extérieur Laplace A Tla, Ê, p)=0 Te Ê., p)= T. (a, 8) 
Ô ð, 
Tala BP)... = Tayla, Bo) aa, z Tula, P, o) = p T (a, p, p) 
2; 2 2 
wee OT OT a a A R 
ox 0y Oz yla, b.p) plap) 


Dans ce cas, il existe toujours une solution particulière de l'équation de poisson que l'on eut 
utiliser pour homogénéiser cette équation en une équation de Laplace. La forme de la solution 
particulière est la suivante : 


Er 4 29) 


Ton =Q = 4A= 0 A4= ST (x2)- 2 (y +y?) 


Si l'on retranche de la solution intérieure cette solution TA alors le problème aux limites 
devient de Laplace et si à l'origine le problème ne dépend pas de l'angle azimutal, la résolution du 
problème devient classique. Toute solution de la forme : 


Tyan(x:92)= Ap° = +»)+B= AT, =Q, 


part 


T (x, y;z )= Be +y +z 2)+ B= AT = Q, 


Pp part 


où B est une constante quelconque est également solution particulière de l'équation de poisson. 


Problème aux limites de Laplace — système de coordonnées toroïdales — p 


Exemple d'électrostatique : problème extérieur de Dirichlet, tore dans un champ électrique 
uniforme le long de l'axe z, porté au potentiel 0, exemple tiré de l'ouvrage de N.N.Lebedev, 
1.P.Skalkaya, Y.S.Ufliand « Problems of Mathematical physics » (problème n° 494 page 236) 


_ aSinh(n)Cos(p) _ aSinh(n)Sin(o) = a Sin(0) 
Cosh(n) - Cos(0) Cosh(n) - Cos(0) Cosh(n) —- Cos(0) 
avec 0<nm<+o 0<0<2m 0<p<2r 
Axez—m=0 z= ce = cn = a Cora £) 
EO 25) 


Cercle limite = n > +% x=aCos(p) y=aSin(p) = rayon a 














Pour un tore d'iso-surface n= n0, si les dimensions du tore sont définies comme suit, r : rayon de la section de l'anneau, 
l: rayon du tore, rayon du cercle au centre de la section, Alors, il vient : 








RIRE a Sinh(l,) a Sin(0) 
P =, x + y = Z'— 
Cosh(l,)— Cos(0) Cosh(l,)— Cos(0) 
a Sinh(l,) a Sinh(n) _ aSinh(I,) 


z=0=> Cos(0)=+1 p=—T— Pun = N= 
Cosh(l,)+1 Cosh(l,)+1 Cosh(l,)—1 


a Sinh(l,) aSinh(l,) _ 2aSinh(l,) — 2a 











Pmax — Pmin — E S 2 =; =2r 
Cosh(l,)—-1 Cosh(l,)+1 Cosh'(l,)—-1 Sinh(l,) 
ja Pii Pa < ENGCON COSG) 
= = =a 
2 Cosh’ (l,)—1 Sinh(l,) 
Cosh’ (l )—1 F2 
r = — a >P -r? SE E ne => Sinh(l Na Cong 
Sinh(1,) Sinh’ (l,) " : De 


D'autre part la solution doit se comporter à grande distance comme le potentiel d'un champ électrique uniforme, soit : 


Lim T(n,0)= Ez = Eja LAC, 
240 Cosh(n) - Cos(0) 
La solution doit donc présenter le même niveau de symétrie que le terme Sin(8) . on développe donc la solution avec 
les termes d'une condition aux limites de type impaire. La solution du problème extérieur ajoute au développement en 
série pour une condition aux limites impaires, la forme asymptotique du potentiel à grande distance (par évidence 
solution de l'équation de Laplace par retour en coordonnées cartésiennes): 

P_,(Cosh(n)) 
T (9,0) = Ea Sin(0) P JCosh(n)- Cos(0) S 8 n- 


A Sin(n9) 
Cosh(n) - Cos(0) m aa Paa (Cosh(l,)) 
Le | 











Avec : 


Pa fas f,(0) Sin(n9) -afas f,(0) Sin(n9) 
y JCosh(l,)-Cos(0) œ „/Cosh(l,) - Cos(0) 














La condition aux limites à la surface du tore au potentiel nul devient alors : 
Sin(0) E JCosh(l,)- Cos(0) 




















T(n,0 =E B Sin(n9)= 0 
á Vo, Cosh(l,)— Cos(6) T 2 POr 
oem Sin(0) a 2f d9 HO Sin(n9) 2E,a]d9 Sin(0)Sin(n9) 





Cosh(l,) — Cos(0) o  „/Cosh(l,)— Cos(0) o (Cosk(I,)- Cos(0)} 
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En utilisant la formule d'intégration définie : 














1f Cos(n0 1 7 Cos(n0 
Q 1 (Cosh(m))= gle ae) ie D d0 20) T 
g 2 (Cosh(m)-Cos(0)} 22o  (Cosh(m)-Cos(0)} 
= 5 d0 Cos(n8) I pour n = Q,1,2,3,... 
242»  (Cosh(n)- Cos(0)} 
et en intégrant par partie, il vient : 
f Sin(0)Sin(n9) f 1 


d =-2{49 | 
o (CosA(I,)- CosO 0 | (Cosh(l,)- Cos(0)} 








Sin(n9) 


Sin(n9) f nCos(nẸ) 
f d9 


(Cosh(,)-Cos(0)} |, * (Cosh(l,)- Cos(0)}: 











e Pong) =2V2nQ 1(CoskG, ) 
o  (Cosh(l,)- Cos(0)} 
= [as Sin(0)Sin(n9) =24/2n0 (Cosh, )) 


0 (Cosh(l,)- Cos(0)} 
Ce qui donne comme solution, finalement : 


Q (Cosh(l,)) 


pails 





Sin(@) 


T(n,0)= Ea ID à À Ja [Cosh PA Ee T) P e 





A profil sur l'axe z est le suivant (n=0) : 
| Q (Cosh(l,)) 
T(0,0)= E,a Sin(0) — 4E 4 2A Cos(0) £ "P Cona) Sin(n9) 


° 1-Cos(0) ~x tre 





Q (Cosh) 


T (0,80) = E,a Coran $) Siaa Sin( 2) n Sin(n9) 
2 mT 2e VE (CosA(I,)) 
2 


Avec les données l,r, dimensionnelles du tore, la solution s'écrit : 
Q (Cosh(l,)) 


8E VE -r° n-i , 
T(0,0)= EN’ -r Coral $ | A KOPU en 


a o (Cost, ) 
= T(0,0)= EN -r Coran $) — Es) » "P (Cond) Sin(n9) 
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Exemple d'électrostatique : problème extérieur de Dirichlet, tore dans un champ électrique 
uniforme le long de l'axe z, porté au potentiel VO 


A la solution précédente, il suffit d'ajouter par principe de superposition, la solution du tore au 
potentiel fixe VO sur la surface : 


Q (Cosh(l,)) 


Sin(0) 4 — -1 | 
none PDG OPEN 


n== 
2 








O ,(Cosh(l,)) 


JCosh(n)-Cos(0)| -3 
+72 7 o 1 (Cosk(,)) 


t 2 ET) Cos(n9) P (Goshe) 

















Exemple d'électrostatique : problème extérieur de Dirichlet, tore dans un champ électrique 
uniforme le long de l'axe x ou y, porté au potentiel 0 


4 Sinh(n)Cos(@) _ a Sinh(n)Sin(@) 22 a Sin(0) 


Cosh(n) - Cos(0) Cosh(n) - Cos(0) Cosh(n) - Cos(0) 
avec 0<n<+o 0<0<2m 0<p<2x 











D'autre part la solution doit se comporter à grande distance comme le potentiel d'un champ électrique uniforme, soit : 
a Sinh(n)Cos(@) 

s Cosh(ņ)— Cos(0) 

La solution doit donc présenter le même niveau de symétrie que le terme constant 6 soit Cos(8). De plus on ne retient 

de la dépendance azimutale que le terme en Cos(b). On développe donc la solution avec les termes d'une condition 

aux limites de type paire. La solution du problème extérieur ajoute au développement en série pour une condition aux 

limites paire, la forme asymptotique du potentiel à grande distance (par évidence solution de l'équation de Laplace par 

retour en coordonnées cartésiennes): 





Lim T(n,0)= Ex = 


P' ,(Cosh(n)) 








T(n,0,p) = Ea SAGESSE) + yCoshtn) TOONE) D ô,B, = Cos(n 9)Coslo) 
Cosh(n) — Cos(0) 27 ie Pi (Cosh(l,) ) 
RS 


6, =1 6,=2 n>0 


Ja (0) Cos(nÿ) 
JCosh(,)-Cos(0) 








=2 Í d9 
Avec : 
La condition aux limites à la surface di tore au potentiel nul devient alors : 


Sinh(l,)Cos(p)  „JCosh(l,) - Cos(0) E 
nọ, 7 E04 Cosh(l,)- Cos (0) 2x 2 PS PRE? 





T(n7,0,@) 




















Sinh(l,) fa (0) Cos(n9) _ 2E,aSinh(l | d9 Cos(n9) 


9(0)= -E;a „=2 3 
fol ) Co sh(L, )— Cos(0) + l JCosh(1,)- Cos (0) 0 (Cosh(l,)- Cos(0)} 
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En utilisant la formule d'intégration définie : 


z Q° ,(Cosh(n)) 
d9 Cos(19) me 2/2 ne 
0 (Cosh(n)- Cos(9)} Sinh(n) 
Ce qui donne comme solution: 
B,=4V2E;aQ" ;(Cosh(l,))  &=1 8,=2 n>0 
2 





Sinh(n)Cos(o) 
°” Cosh(n)-Cos(0) 





T(n,0,p) = 





L2Esa V2 /Cosh(m) - Cos(0) 3 





2 


Q' ,(Cosh(l, )) 





m E N P' ‚(Cosh(l,) 
2 


Si le champ est dirigé selon l'axe y, il vient : 
ô =1 ô, =2 n>0 


Sinh(n)Sin(o) 
°” Cosh(n)-Cos(0) 








T(n,0,p) = 


25a J2 /Cosh(n) - Cos(0) 





j P' ,(Cosh(n)) Cos(n9)Cos(p) 


Q' (Cosh(l,)) 





2 
m Er P' ,(Cosh(l,) 
2 


Pr; (Cosh(n)) Cos(n 9)Sin(o) 
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Et si le tore est porté au potentiel VO<>O , il vient : 
Champ éléctrique axe x 


Sinh(n)Cos(o) 
g Cosh(ņ)— Cos(0) 





V2/Cosh(m)- Cos(0) 5 z o 


m P' ‚(Cosh(l,)) 











T(n,0,@)=4+2E,a P (Cosh(n)) Cos(n9)Cos(p)+ 
2 


O ,(Cosh(l,)) 

P 
Q  (CosA(l,)) 

PET a n (Cosh) 


ı(Cosh(n))+ 





+y, V2,/Cosh(n) - Cos(0) 
T 











Champ éléctrique axe y 


Sinh(n)Sin(o) 
°” Cosh(n)-Cos(0) 








Q" ,(Cosh(l,)) 
Le: P, 
P' (Cosh(l,)) "=; 








V2 /Cosh(n) - Cos(0) 5 


T 


(Cosh(n)) Cos(n9)Sin(o)+ 


T(n,0,p)=<+2E;a 


O , 1 (Coshü, ) 


P Vont ) 


ı(Cosh(n))+ 








" V2 /Cosh(n) - Cos(0) 





T 








72 


+25 


n=1,+00 


O \(Cosh(,)) 


FC) Cos(n9)P. AA) 
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Exemple d'électrostatique : problème extérieur de Dirichlet, potentiel extérieur à la surface d'un 
tore conducteur parfait induit par une charge placée à l'origine des coordonnées et donc porté 
au potentiel 0. Exemple tiré de l'ouvrage de N.N.Lebedev, I.P.Skalkaya, Y.S.Ufliand « Problems of 
Mathematical physics » (problème n° 495 page 237) 


La solution additionne le potentiel crée par la charge à la solution général du problème aux limites 
avec une fonction limite paire en 9. Il ne reste plus qu'à identifier les termes du développement par 
la condition aux limites à la surface : 



























































T (1.0) = q _ _4(Cosh(n) - Cos(@)) _ qy Cosh(n)- Cos(0) 
Fun JL +y +z ay] Sinh? (N) + Sin? (0) a]Cosh(r) + Cos(0) 
F - fao (0) F >| Jg LO) Cos(n9) 
° ẹ  JCosh(,)-Cos(@) " œ  „fCosh(l,)- Cos(0) 
Cosh(n)— Cos(® 
raa Ces woe i Je Coro) X P eeng) COPE, (Coshon) 
Si T(n,0)=T,(n.0)+T.(n.0)= T.(n.0),, =-T,0,0) , =-4 LOD O ig 
n=l, n=} a |Cosh(l,)+ Cos(0) 
f Cos(n0) B f Cos(n0) : ey 
l JCosh{n) — Cos(0) 220 OS a JCosh(n) + Cos(6) 220) Q i (Coston) 
__af l __4 
TAS n JCoshG,)+ Cos(@) 2929 (Cosh) 


y 
œ 
I 








q% Cos(n9) o q . 
LA l di JCosh,)+ Cos(@) aa 1) 9, (Cash) 
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Ce qui donne la solution : 


q\Cosh(n) — Cos(@) 
a]Cosh(n) + Cos(@) 











o , (Cosh) 


T(7,0) = 2P «Cost, De ot) 


q 3 jz VCosh(n)— Cos(6) 
a m 








Q (Cosh(l,)) 
j Cos(n9)P ; (Cosh(n)) 


+ 5 (-1) 


n=1,+00 








P (Cosh (l) 








Ge gla 





q-|Cosh(n) — Cos(0) 
a\]Cosh(n) + Cos(0) 








0, (oska, )) 


= T(n,0)= P Cona, Ja (Cosh) 


q 3 VCoshn) — Cos (0) 
a 7 Oo ,(Cosh(1,)) 


+2 $ -1 - C0) GTR) 














n=1,+00 








Le profil sur l'axe z est le suivant (n=0) : 
































g\1-Cos(0) 

ajl + Cos(0) 
T(0,0)= TON Q (Cosh, )) O (Cosh(l,)) 

t2 = ne +2 > ( 1)" 72 Cos(n9) 

P Cona, I) na P, (Cosh(l,)) 

/ j s2) 0, ,(Cosha, )) Q_:(Cosh(,)) 

= T(0,8)=% Tarf $ 2 > P (Ces, 5 ot 1) P (Cosh) Cos(n9) 
à sin(9 0 : O ,(Cosh(l,)) O, (Cosh(l,)) 

=> T(0,0) = PF ran $) 2 P (Co, j) eA 1) > Cosa) Cos(n9) 
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Exemple régime stationnaire équation de la chaleur en présence de source de chaleur uniforme: 
problème intérieur de Dirichlet, à l'intérieur d'un tore dont la surface est à une valeur nulle avec 
une répartition uniforme de source de chaleur. Exemple tiré de l'ouvrage de N.N.Lebedev, 
1.P.Skalkaya, Y.S.Ufliand « Problems of Mathematical physics » (problème n° 498 page 238) 


Il s'agit d'une équation de poisson. On utilise la solution particulière pour homogénéiser l'équation 
de Poisson. 


Qa Sinh’ (n) 

4 (Cosh(n)- Cos(0)Ÿ 
On retranche cette solution particulière pour homogénéiser l'équation de Poisson et la transformer 
en équation de Laplace : 


AT, 4 = -0, > Talay) = -2e +» )= 








Cosh(n) - Cos(@) B, 
Tion (1,0) = - 3 O Cona) A) PER ı(Cosh(n)) 


2 





T(n,0) = Tanx y, z)+T, hom (n, 0) 


Qa? Sink? (l,) 
T ,0 = 0 Trom 379 D Le 0 = 
(n J > Tomy 2j, par (7 Jia 4 (Cost, )- Cos(0)} 





























„Qa? Sin”) _ fa Jo(0) Cos(n9) 
=)= 4 (Cosh(l,)- Cos(0)¥ e l aR [Cosh(l,) -Cos (0) 
B, = Qa’ Sinh’ (1) T Cos(nĝ) 
4 0 (CosA(, )- Cos(0)} 
De plus  Q",(Cosh(n)) = Sin (nm + J d0 CU r Mmentier 
"z 27 2/0 (Cosh(m)-Cos(8))" "7 
e 1) ur Qm-i} F Cos(n0) 
o" (Cosh(n}= ED sinn”) L= fao | 
Re: V2 (m1) 4 (Cosh(n) - Cos(8)}"*z 








3Sinh (7, Jz Cos(n0) 
> m=2 0° 1 (Coshll, )- Li : 
z 8V2 à (Cosh(n)-Cos(0))} 








= fao — LD 82 o (co) B, = 2220 2 (Ci) 
0 (Cosh(n)- Cos(0)} 3Sinh (n) F2 l 3 
Sink? (n) 
4(Cosh(ņ)- Cos(0)) 
o°, (Cosh, )) 
2 Q ACen 


T(N,0) = -Q,a° i gi Cosh(l,) 
— — y Cosh(ņ) — Cos (0 
DREO T a 


+ LÈ o (Con) 2: Q (Cosh(n))Cos(n9) 














Problème aux limites de Laplace — système de coordonnées toroïdales — p 643 


Exemple gravitationnel: problème extérieur de Dirichlet, champ gravitationnel d'un tore de 
densité massique constante p. Exemple tiré de l'ouvrage de N.N.Lebedev, I.P.Skalkaya, 
Y.S.Ufliand « Problems of Mathematical physics » (problème n° 499 page 238) 


Il s'agit d'une équation de poisson. On utilise la solution particulière pour homogénéiser l'équation 
de Poisson. 


2 Sinh’ (n) 

(Cosh(n) — Cos(0)Ÿ 
On retranche cette solution particulière pour homogénéiser l'équation de Poisson et la transformer 
en équation de Laplace : 








AT sort = ATP > T ar (x,y,z) = pe + x)= -mpa 















































Re Cosh(n) - Cos(0) À À 
T°" (n,0 _\ 0 Cosh n Cosh(n))Cos(n9 
w” (7,0) > A ACTA) os DD o (Cost) CA osh(n))Cos(n 9) 
2 2 
Sinh’ (n) 
T ,0 =T Ji Thom ,0 E 2 Thon ,0 
int (n ) “a Yz z)+ int (7 ) npa (Cosh(n)- Cos(8)) ii int (7 ) 
JCosh(n) — Cos(0) B 
T..(n,0) = o (Cosh K h) )Cos(n9 
o (1,0) f P(Cosh(, j osh(n))+ R D Poa- Cosh) osh(n))Cos(n9) 
E 
0 0 
T0), =La, et se _ OT: (0,0) 
n n=l on n=l 
puisque O) TaD) Cosh, )-CoslO) aTa (7.0)  _ Coshl,)-Cos(O) aT., (7.6) 
On nai On a ôn ya} a ôn jar 
O°, (Cosh(n)) 
. 2 _ De. 
pose Sinh? (n) Le 82 {/Cosh(n)- Cos(@)| -; + Y © ,(Cosh(n)}Cos(n9) 


(Cosh(n)- Cos(0)) 3 T 2 ar E 
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Par application des conditions de continuité on aboutit au système d'équation linéaire par 
identification des termes en Cos(nô) : 
Cosh(n)- Cos(0) 























Posons A(n) = 
T 
©, (Coshtn))= -ap ÈS 0? (Ces) 2, (Coshtn))= -roa ÉS2 0: (CosHn) 
l in=0 
Tm (1:0) = A(n) XE, 2,(Coshn)Cosn9) e, =42 "7 
PRE 1 sin>0 
Q (Coshn)) 
Ta (1,0) = A(n D e, Q,(Cosh(n))+ 4, O (Cost) Cos(n9) 
> 2 
P_(Cosh(n)) 
T (1,0) = A(n Sr) SE 
Continuité > Q, (Cosh(1,))+ A,=B, Continuité de la dérivée première D _ = a 
To Lu, To lu, 
Q \(Cosh(n)) 
— À MX e, Q,(Cosh(n))+ A, Q Z Cosh(,) Cos(nà) + 
O (Cosh) 
+ aP Q, (Cosh(n))+ 4, © C0) Cos(n9) = 
' B, i 
= À HZ. Cash) EP. CAGE) FA) ? FC. (Ca) £,Cos(n8) 
Q (Cosh(n)) x 
Q, (Cosh(n))+ À n o 7 (CosAt, ) = P Cord, e- (Ca) 
0 | (Cost) P 1'(CosA(,)) 

> Sinh(l,)| Q,'(Coshli,)+ 4 T Cot, 5j = Sinh(1,)B, AU) 


Q, (CosA(1,))+ 4, =B, 


(Cosh(,))  P ,'(CosA(,)) 


Q ; 1 
7 1a, '(Coshli, )+ 4, —— sp 
ý ANNE (Cosh(,))  " P_,(Cosh(l,)) 


n 
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On résout aisément ce système d'équations linéaires : 





O ,'(Cosh(l,)) P_;"(Cosh(1,)) 
B,=Q,(Cosh(,))+ 4,  Q,'(Cosh(l,))+ A "o at, a = (Q,(Cosh(l, 7) 


al 
2 


= 
~ 








1'(Cosh(, U 7 P '(Cosh()) 
Kenu a 


| E 


Ei 
Cost) Q Ci) 


Pa 
2 








P, '(Cosh(l,)) Q "(CoshG,)) 
SU ne: Cet, ) © Cost, D Sin GP. 1 (Cosh, D r(Cosh, )) 
et Q (Cosh, ))P_ ;'(Cosh(,))= = — T ao a ) Gn )) 
P P (Cast, o, (Cosh(, )-Q,(n)P ;'(Cosh(i,)) 
Sinh’ (l,)P Cosna ) 4 (Cosh(,)) P (Cosh(1,)) 


A, = Sinh?(1,)Q Gant). 1 (Cosh), (Cosh(,))- P a (Cosh4, ), (Cosh(, )] 


=> 


B, = Sinh? (l,)P (Conty |o i (Cosh), (Cosh(,))-Q PTT 
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Les solutions intérieure et extérieure s'écrivent donc : 


a,= P, (Cosh(l E ,'(Coshli, )- R. à (Con, o? 1 (Cosh, ) 
b,=Q (Cor) ,'(Coshl,)-0 ,'(Coshg, »0° à (Covi, ) 


1 sin=0 
Posons ô, = 2€, = | 
2 sin>0 





T 0,0)= -pa 22 [Cost )- Cos(0) D afo, ı (Cosh(n))+ Sink’ (l)a, o (Co) Ent 


T..(n,9) = -pa° TA )- Cos(0) > ôb, P CAGE EE) 


BC) 


Pr 1 (Cosh(, )}= PÔ 2 (Cosha, ))= En 


0" et ) 








0° (Cash) = Sint? (n° (Cash) 
O ,(Cosnn)) = Sin (9°) (Cosh(r)) 
2Sinh(n )Cosh(n 2°, (Cosh (n )) + Sinh? (n o“: (Cosh(n ) 


= O° ,'(Cosh(n))= : Sinh(n) i 


2Cosh(1, 2 ï (Cosh (4, ) + Sinh(l, Tos à (Cosh(1, ) 


Sinh°(1, ) 








=> Q? ;'(CosA(s, )= 


2, 


EF (conte) 2costl1,)0: (Cos, J+ sinali, JO, (Cost, ) _ Simi,)P! ,(Cosali, o? , (Cesu, ) 
P (Cosh | 2Costl1,)0: (En, J+ sinali, JO , (Cost, ) - sinat, JO, (Cosh(,)o? ,(Coshl,) 








Tu (1,0) =-pa° 42 VCosh(n)- Cos(0) $, 5 0? (Cosin )+a,0 (Ces) ]Cestns) 


n=0,+00 








12307; 0)= — pa” 





2 VCosh{n)-Cos(0) X. Sbp P COHERE 


n=0,+00 
Dans un article récent de 2018 « Exact gravitational Seal of a homogeneous torus in toroidal 
coordinates », l'auteur M.Majic utilise une autre solution particulière de la forme : 


> Sinh? (n) + Sin? (0) _ _, Cosh’ (n) — Cos?’ (0) _ ; Cosh(n)+ Cos(0) 
(Cosh(n)- Cos(8)) ne (Cosh(n)- Cos(8)) Dé Cosh(n)- Cos(0) 
Cosh(n) 
Cosh(n)- Cos(®) 





xX +y +2 = 





xX +y +z +a’=2a? 


Arpa * Cosh(n) 
3 Cosh(n)- Cos(0) 





AT = 40 > T (ny) + y 42 a?) 


part 


=> 
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Ce faisant on, a abaisser d'un degré le dénominateur de la solution particulière. La décomposition 
de Fourier de cette solution particulière est la suivante 


T0) = Conte — Cos(0) 5 



































E,B,Cos(n9)Q (Cosh(n)) ô, =2e, 
Z joje Arpa ” Cosh(n) nee Arpa “Cosh(n) F jo Cos(n9) 
3  Cosh(n)-Cos(0)  ‘” 3 A (Cosh(n)- Cos(0)} 
De plus Q",(Cosh(n))= Ep" Sinh” (n) Cm 1} fao La r mentier 
a v2 rm 5 (Cosh(n)- Cos(0))""3 
1 Sinh(n) $ Cos(n0) 
=>m=1 Q' ,(Cosh(n))= d0 z 
E 442 l (Cosh(n)- Cos(0)} 
f Cos(n0) 42 16xpa °Cosh(n) | 
d9 -= Cosh B, = - V20' (Cosh 
k (Cosh(m)-Cos(o)p SA 2 ES 3Sinh(n) 9, 1 oh) 
D HO) = 7e *Coshln ) jz VCosh(r)- Cos(8) Dao 1 (Cosh(a DCos(n9)Q ,(Coshon) 
á 3Sinh(n) m n=0, +0 
O' (Cosh(n))= Sinh(n)O” , (Cosh(n)) 
D' 2 2 
A O? (Cosh(n)}= Sinh (no 


l (Cosh(n)) 


2 ar an D y% 
\désignent les dérivées première et seconde 











> T,,:(7,0)= PP *Cosh(n l pacan Ca à 2:07 KERN )CosnDO s (Cash), 
,(Cosn()= PEN? 


Cosh(n X (Cosh(n))= CRE Co |  (CosA(n)) 
(Coshty))= SP 2 ego, (Cosh{n)}+ Cash(}® (Cosh( ) 
Q, (CosA(n j= 8mpa Lf 1 



































OL CoH) SA o: (can) 
Ou bien en utilisant l'équation de Legendre 
: o, (fn -1) -200,6 : O, (Cosh(n) 7 ] = 2Coshtn)o (Coshtn)) 
o” (E)= = CURE nm g ii (Coshtn))= a Sinh’ (n) i 
(Cash) = E72 PACE mfo CO E 0000) 
Q, '(Cosh(n)) = | os 0 i (Cosh(n))+ ne ro, l (Cosh(n ] 
©, (Con) SRE CO (Co) 
Q,,'(Cosh(n)) = | Lee lo i (Carn: Saren ne re (ant) 


Problème aux limites de Laplace — système de coordonnées toroïdales — p 648 


Solution sur une lentille sphérique symétrique assimilable à un problème homogène angulaire 


Supposons que l'on complète la lentille sphérique trouée définie par le domaine : 


A (o0) ef] e za) 





0.0 


Et que le domaine complémentaire représente la zone en bleu sur le schéma ci dessus. 


Q = fo) nel, }8e Fr] 


On définie le problème aux limites suivant sur la réunion des deux domaines : 








AT(,0)=0 (7,60) =fel0+el0efr-a,7r+a) à = à on. 
TNO o] = Lo T(7,0) neli, 0] | j 
TOY pa s T T(7,0) melo.+e] | j 


O=r+a0 O=r+a0 


Soit donc une valeur TO sur la surface de la zone centrale de la lentille, O ailleurs sur la surface 


extérieure. 
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Sur le domaine: Q,-{7,0)/ne [m00 efr -a,,x +} le problème aux limites est 
homogène en angle avec une fonction limite constante à définir, il se ramène donc à la catégorie 
de problème de cette section. Le problème aux limites sur le domaine 
Q= {7,0)/n € [0,n, | € [z -QT +a l} se ramène à un problème homogène angulaire en 
retranchant la solution constante comme ceci : 

T(7,0) = T + To (n.0) 

Ta (1,0) =0 Ta (n.0) =0 


nelo 0 ] 
O=7+a0 





ne[o “no ] 
0=7-a0 


Il devient alors évident que l'on a déjà rencontré ce genre de problème aux limites dans différentes 
configuration géométriques : système polaire, sphérique, sphéroïdal allongé et aplati. On va donc 
utiliser le même procédé de construction. 


Le problème aux limites se résume donc au suivant : 
T1(7,0) = T g Toa (1,0) 


Ta (7,0) z 0 Toa (n,0 Ju no] z z 


T,:(7,0) = 0 To(n, 0) ia 

ndomol n0 

0=7-20 s- z+a 
Avec comme contrainte le respect des conditions de continuité et de continuité de la dérivée 
première à l'interface des deux domaines : 


To + Too 0) = Toz llyo0) 
Cosh(l no) — Cos(0) aTa (l n0? 0) _ Cosh(I 0) Cos(0) OT, (l 10? ADN aTa (l 10? 0) _ To n0? 0) 
c ôn c ôn ôn ôn 


La solution dans les domaines respectifs se développe en série de sinus comme suit : 
P (Cosh(n)) 




















E ĮJCosh(n) - Cos (0) Ar | E 
Toa ,0)=T, + a 2 A, P. Con Pa (0+a,—x) 
Q (Cosh(n)) 
_ „[Cosh(n) — Cos (0) a- f G 
To: (1,0) = =. 2 B, TE Conio) A (0+a,-x) 


Il faut également trouver une formule de décomposition du terme constant TO (soit l'unité), à 
partir d'un problème aux limites simple sur [0,In0] sans le calculer Ines nous avons: 


























Cosh(n)- Cos(0) Sin(A,(8 +a,-x)) 
1= U (n)Sin(A, (0 +a,- ,(Cosh(n d9 0 

A 21 ()Sin(2,(0+a,-x)) U,(Cos D-T Tosha) CO) 

2a 
s m í Sin(,0) 
Sin(22,@)=0= À, = U,(Cosh(n d9 i 
EEE Em Aos a JCosh(n)- Cos(0 + 7-a) 
2% x 

U,'(Coshtn)=-} [ a9 SA) 


0 (Cosh(n)- Cos(0 + x -a, » 
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Ce qui donne le système linéaire suivant en appliquant les deux conditions de continuité : 
P (Cosh(n)) 











Cosh Cos(0 
Ta (7,0) = vCos cp- AD gle FES LA (Ci, ) 


ba (0+a,-x) 




















Q ı(Cosh(n)) 
Cosh Cos(0 An : Cosh(n) -— Cos(0 
To (1.0) = Je CEA, Di o ACTU APR) AN a 2e) 


P _1(Cosh(n)) 
T(7,0)= A0) È | U,(Cosh(n))+ À "P CA) Sin(2,(0+a,-x)) 


n=1,+00 
An- 





Ses ps 





0, _ (Cosh(n)) 
T(n,0)= A) X. B Sin(4,(0 +a,- 7)) 


2 
n=1,+0 4 Q (CosA(,) 
Ta 


Continuité simple = U,(Cosh(n))+ 4, = B, = 4, =B, -U, (CosA(,)) 
Pi | s'(Cos, ) Q '(Cosh(l,o)) 














Continuité de la dérivée = A, P TO T, -) +U „'(Cosh(l,))= B, o F (Cosh T) 
2 "2 
P sos) O 2, Eii '(Cosh(i 7A) P i '(Cosh(la)) 
PD 2 2 
SRG 7) A Le 
"2 2 n55 


De plus P "Cosh(n), a (Cosr()- T ,(Coshll, )0, E '(Cosh(l,))= Ee Wronskien 
: a : d 


4, Si) Cosita) v, (cos), (ce )-0, à CAC) 


Soit i 


B,- siml), (Emi). (os), (cos,)-0, (cos, (cils) 


2 
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Injecter dans les solutions cherchées, il vient : 


a, =U, (Cos(,,)0, : '(Cosh(1, a))- Oi ,(Coshll, JU, '(Coshll,o)) 
b, =U, (Cosa(s,.)P, à (CosA(l,,)-U, (Coshll, o)P, (Cosato) 


























Ta, (7,0) = L Ce) > (v. (Cosh(n ))+ a,Sinh° (l,o )P. (Cosh) |sinla,(0 +a, _ x) 
T(7,0) = Sinh? (l, lens SOE) > bO 1 (Cosh(n)}Sin(a, (0+a,-x)) 
“a Sin(,0) 
U,(Cosh(x) )= | jé JCosh(n)- Cos(0 + x - à) 
U,’ (Cosin) =} f as Sin(1,9) 


0 (Cosh(n)-Cos(0 +7 -a, DE 
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Problème intérieur de Dirichlet à trois dimensions avec dépendance azimutale : équivalence des 
solutions en R-Séparabilité-1 et R-Séparabilité-2 


La solution de R-séparabilité-1 en trois dimensions donnent pour une condition aux limites paires 
en Ÿ et p: 


Tn.0,0) = JON = Cas > 420 (Cash }Costn Come) 








T(n,0,0),. = f(0,p)= Er) CR | D >, 4,97 (Cosh, }enstnshcax(mo)| 


n=0,+00 m=0,+00 


2 Cos(nd) 
Ann = e JCosh(l,)- Cos(0) à 





i do f (0,9) Cos(m) 





O” (Cosh(n)) 
_ VCosh(n)-Cos(0) Anm n-i 
PAPE r? 2, 0 du) O, Cora, DE RTEINOE 











Si la fonction limite ne dépend que de w, il vient : 


Cos(n9) 2€ 
fas JCoshll,)—Cos(0) G laer Gin ki 


fao Cos(nĝ) 
JCosh(l,)- Cos (0) 








Anm =B,Cn > B, = 


Or B, = 








= 22 20 1 (Cosh(, )) 





Q (Cosh(,)) Q” , (Cosh()) 


a JCosh(n)-— Cos(0) n-i i 
T(n,0,ọ)= 24/2 7 2o 2 CS ETES O", (Cont. j) Cos(n9)Cos(mop) 








La solution de R-séparabilité-1 en trois dimensions donnent pour une condition aux limites impaires 
en Ÿ et p: 


Tn.0.0)= CON = CoO `y >: 4,07 , (Coship)Sintn BISG) 


n=0,+00 m=0,+00 








T(n,0,p),., = 0.0) JC) CO | >. >, DA 0) 


Sin(n 9) f | 
= | d9 d 0,p)S 
An = J to | p f(0.p)Sin(mp) 








O” ,(Cosh(n)) 
_ „Cosh(ņ)- Cos(0) 4,» n- | | 
MERS = an Eiri) O Cong) COO 
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Si la fonction limite ne dépend que de ọ, il vient : 


[us SRE Ca Jde r@) Sin(mp) 








Anm = ByCn > B, = 





Q” 1 (Cosh()) 


>; >; B,Cu 2 Sin(n9)Sin(mp) 
venais l+ Sno 8,0) O” (Cosh(,)) 
2 








T(7,0,0) = dCosh(n) = Cos(0) 


R-séparabilité 2 : problème intérieur dépendant des variables n, et ọ : 


D'après le tableau des solutions finies la solution s'écrit pour des conditions aux limites paire en f 
eto: 


T(n,0,p)= ul DD AP , Cotanht))Costn 9)Costmo) 








n=0,+00 m=0,+00 








T(n.0.p),., = f(0.,9)- ee D OX 4,,P",(Cotanh(l,)Cos(n9)Cos(mp) 


Sinh(l ) n=0,+00 m=0,+00 


Cos(nd9) 
A = Íd 
ti ] JCosh(i,)- Cos(0) 5 








1 do f (6,6) Cos(mp) 





P” ,(Cotanh(n)) 
2 jo 9)Cos(mop) 








Cosh(n)- Cos(0) jt Ain 


T(n,0,0)= 
ERP T° Sinh) | pon (148, (1+ 84,9) P", (Cotanh(,) 
#3 


Si la fonction limite ne dépend que de ọ, il vient : 


Cos(n9) 2x 
Jas JCosh(, )— Cos(0) Cn = Jde FE) Cos(mop) 


Cos(nd) 
JCosh(1,)-Cos(0) 


Anm = ByCn > B, = 














Or B, = - do =2V20 (Cosh(,)) 





O, (Cosh(l,))Cņ P" ,(Cotanhn)) 


JCosh(n)-Cos(0) [Sinh(l,) "i E 
TG:8,9)= 292 m Snh)| tante (+8,69)1+6n0) P” Corana) OO 








avec la formule de Whipple 


r|m -n+ -) 2 y 
Q", (Cosh(l, )) = = | | P", (Cotanh(, )) 


(men P 2) m Sinh(l,) 
2 
fe or emm 


AGDE Cos(0) Cm x 














xCos(n9)Cos(mop) 
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En appliquant une fois de plus la formule de Whipple sur le terme : 





Piment A 
o” ,(Cosh(n))= ) - ) P" (Cotanh(n)) 
E rm +n+ 5) z Sinh(n) E 


ne O ,(CosA(,)) Q",(Cosh(n)) 


1 
Cosh(n)- Cos(0) ps = 
rm none” O” (Cosh, )) 0+8, o+, o) 


n 
s 
2 


Il vient : 





T(n,0,p)= 








Cos(n9)Cos(mop) 


Ce qui est exactement la formule obtenue avec la R-Séparabilité-1 ! 


La solution pour une condition aux limites impaire en Ÿ et ọ se développerait comme suit : 


rono,py= AD Ce 2 2 AP (Conan) Sin) 














CIE ER) D D 4,,P",(Cotanh(l,))Sin(n9)Sin(mp) 


T(n.0,0),., 00 | Sinh(l,) Oo m0 : 
n n=0,+00 m=0,+00 2 








f Sin(n9) 2€ | 
A, „= (d9 d 0.) S 
” | Tert) -C | P fa (0, p) Sin(mp) 





? Pry (Cotanh(n)) 


Sinh(l 
A. >» > mm 2 Sin(n9)Sin(mo) 
Sinh) | vente A +8,0)1+8,,0) P" (Cotanh(,)) 
2 


éd VCoshtn)- Cos(0) J 


T 








Si la fonction limite ne dépend que de ọ, il vient : 


T Sin(n 9) ay ; 
A,m = B,C, > B, = |d C,=l|d Sin(mp) 
L JCosh(,)-Cos(6) l Pa SMON 











Pa (Cotanh(n)) 
72 JS 9)Sin(mp) 








JCosh(n)- Cos(0) Fe 5 B,C, 


T ,0, = 
nE 7° Sinh) pona (148, o) +8,0) P, (Cotanh(l,) 
tre, 





rfr -n+ | > 2 
O” ,(Cosh(l,))= 2 | | P” | (Cotanh(1,)) 


mn r{min,) m Sinh(1,,) 
2 
rm -n+ | Prg (Cotanh(n)) 
i JCosh(n) -Cos (0) B,C, 2j) m- Sin(n9)Sin(mo) 


m? V7 Sinh) | ent (+ +80) {men 5) Q"  (Cosh(,) 


m 
1 
n—— 
2 


1 








T(n,0,@) z. 
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En appliquant une fois de plus la formule de Whipple sur le terme : 





rm n+) > pA 
Q” ,(Cosh(n))= - 4 ) P” (Cotanh(n)) 











T Ho rsang) 
2 
Il vient : 
Q" ,(Cosh(n)) 
T(n.6.0)-% CAGE CC B,C, "3 Sin(n9)Sin(mp) 
x natremo (1+8,0)1+ n0) O” ı (Cosh(,)) 
2 


Ce qui est encore exactement la formule obtenue avec la R-Séparabilité-1 ! 


Bon, en fait point de magie, si l'on tourne la formule de Whipple comme ceci : 


r|m- n+) 2 y 
Q", GARS : = à (rs) A 
72 


j y 
| | P” (Cotanh(l,)) 


) m Sinh(l,) 


x o Con E Si nh(l, Re ot) 
Q” ,(Cosh(l,)) Ņ\ Sinh(m) P" , Coi, )) 
nz m-5 
C'est une formule redoutable reliant la R-Séparabilité-2 à la R-Séparabilité-1. 


P E (Cotanh(n)) 
2 a 9)Cos(mop) 








ie n+ ; 
ora (Cosh(l, ))= i 
T i 
(m3 














T(n.0.0) = VCosh(n)-Cos(0) |Sinh(l,) Ann 
a m’ Sinh(n) n=0,+00 m=0,+00 (+6, 0) +0,,0) P” , (Cotanh(, ) 


1444444444444444444444244444444#4444444444443 





Q" , (Cosh(n)) 


n,m - 2 Cos(n 9)Cos (mp) 
sonate (L+ 6,0) 0,0) Qi (Cosh(l, )) 


2 


T(7.,6,0) = EE Cos(0) A 











4 P", (Cotanh(n)) 


Sinh(l 4 
T(7,0,p)= 3 4 0) y > Aam 2 Sin(n9)Sin(mop) 
z Sinh) | „one (+8, 1+8, 0) P” (Cotanh(l,)) 


14444444444444444444442444444444&4444444444443 
U 


JCosh(n)-— Cos(0) 











Q” ,(Cosh(n)) 
pe 9)Sin(mop) 








T( 0 ) À /Cosh(n)- Cos(0) > > Anm 2 
4 i 4 m’ n=0,+00 m=0,+00 ad + Ô,.0)4 + Ô mo) e (Cosh(l, ) 
2 
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Problème extérieur de Dirichlet : soit à rechercher la solution de l'équation de Laplace à l'extérieur 
d'un tore dont les conditions aux limites de Dirichlet dépendent des angles O et @, puis du seul 
angle 0 : 


AT(m0,p)=0 (m.0,p)eQ={0<n<1,} T(7,0,6),. fini 
T(7,0,6),. = f(6,p) 
ou bien 


T(n,0,9),., = (0) 


Le problème en coordonnées toroïdales se présente donc comme un problème intérieur, du fait du 
paramétrage du système de coordonnées. On doit également tenir compte des contraintes de ce 
problème : 
T(n,0,@) fini,ne comporte aucune singularité, continue et dérivable 
T(n,0,@0) périodique en 0 de période 2x soit T(n,0,@)=T(n,0 +27x,@) 
DEGRÉ) es en 0 de période 2x soit Geo) CPGE PE 0) 
0p 00 00 
T(n,0,@0) périodique en ọ de période 27 soit T(n,0,0)=T(n,0,p +27) 
CPGE) ces en de période 2x soit OT(,6,) L OT(n,0,p+2x) 
0p 3p 0p 








R-séparabilité 1 problème extérieur de Dirichlet 


Cela conduit au choix d'une partie radiale décroissante en n et finie pour la valeur n=0, et des 
parties angulaires possédant des valeurs propres n et m entières, soit une solution qui se développe 
sous la forme d'une série : 


T(n,0,p)= JCosh(n)-Cos(0) 





(s, P (Cosh) X, Y P” (Cosh Cast) +8, Sin), „Cosmo)+ D, Sin) 


n=1,+0 m=1,+00 


Lorsque les conditions aux limites ne dépendent pas de ọ. La solution s'écrit : 





T(n,0)= /Cosh(n)- DUT P (Coshn))+ J, P 1 (CoskGK4, Cos(n9)+ B, Sintno)) 


n=1,+00 
Limitons nous au cas sans ne en ®. ee les conditions aux limites angulaires sont 
paires, alors la solution se développe en série : 


f(6,p)=f(0) et f(-0)= /f(0) 


T(n,0)= /Cosh(n)- UE P\(Cosh(m))+ J, A,Cos(n9)P. 0) 





n=1,+00 


Dans le cas de conditions aux limites impaires : 


f(@,p)=f(0) et f(-0)=—/f(0) 


T(n,0)= /Cosh(n) - CO! © A Sinn) P (Cain) 
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Compte tenu des conditions aux limites, nous avons alors la solution générale pour des conditions 
aux limites de Dirichlet paires, f(8)=fo(-8), sous la forme : 


0©,(0)=1> |0 =2x 









































O, n 
TN,0) = SO) = SEJ CO] 4 P, (Cosh(l,))+ > 4,Cos(n9) P (Cosh) 
f fo(0) f fo(0) 
Fe f JCosh(1,)-Cos(6) l [JCosh(l,)- Cos(0) 
f (0) Cos(n9) f f, (0) Cos(n9) 

B, - as =2f 

JCosh(,)-Cos(0) œ  „/Cosh(l,)- Cos(0) 
T(n.6)= \/Cosh(n)- Cos(0) B, 





72 ‚Cosma, A a 00 ı (Cosh(n)) 


Avec le résultat 
T 27 
0 1 (Cosh(n}= -p fao — OT 2 a — OI — 
"3 2% (Cosh(m)-Cos(8)} 2V20  (Cosh(n)-Cos(0)} 
_ do Cos(n0) | 
202 (Cosh(n)- Cos(0)} 
pour n = 0,1,2,3,... 














Il vient pour la condition aux limites paires : 


Ja(@0)=1 








= fao JCosh(l, - Cos(0) =2V20 ,(Cosh(,)) 


D Í d Cos(nà) 
"  JCosh(,)-Cos(6) 








=2/2Q ,(Cosh(l,)) 





Q (Cosh(l,)) Qı (Cosh(l,)) 
_ r VCosh(n)-Cos(6)| *- i 
T(n,0)= V2 - P Cost, je CAR Cash) Cos(n9)P , eu) 


Normalement cette expansion converge vers une fonction ue et tendant vers 0 lorsque 
n->0. 





La solution générale pour des conditions aux limites de Dirichlet impaires, f.(8)=-fa(-8), se présente 
sous la forme : 
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2 


®,(0) = Sin(n0) = |@, |” =x 





O, 











T(7,0) 





n=l; 


= 1,0) JET O | D ASntn3)P, (Cost) 


8,= [as fa (0) Sin(n9) af as f(O) Sin(n9) 
" a  „Cosh(l,)-Cos(0) œ  „Cosh(l,)-Cos(0) 




















_ VCosh(n)-Cos(8) B, | 
T(7,0) = 7 2 P (Col) Sin(n9) Taekon) 


La solution de R-séparabilité-1 pour des conditions aux limites paires en trois dimensions donnent : 





Tn.0,0) = [Cos CA O| Dr A, Pr (Com) (n 3)Costmgp) | 





T(7,0,6),., = f(0,p)= OR] > >. A, mP” (Cosh(, Eostn 9Costmp) 


n=0,+00 m=0,+00 


É Cos(n9) f 
= |d9 d 0,p)C 
An J TOHID- C Í p f (0.p) Cos(mp) 











JCosh(m)- Cos(0) A P”, (Cosh(n)) 
T(n,0,p)= d - £ >. nm 2 Cos(n9)Cos(mop) 
x no mate (l+ Snol + 67,0) P”, (Cosh(,)) 
2 








Si la fonction limite ne dépend que de ọ, il vient : 














f Ci Ke! 2m 
Anm =B, ,Cn > B, = fas TET T Cn = Jao) Cos(mọ) 
f Cos(nĝ) 
ee y =2420 ,(Cosh(l 
"ne | JCosh(1,)-Cos(0) 20, l os (,)) 





Q (CoshG,)]C, P”,(Cosh(m)) 


2 — Cos(n9)Cos(mp) 
n=0, +0 À. a + Ô,.0 JA + Ô mo) P", (Cosh(l, )) 
2 








z VCosh(n)-Cos(0) 
T(n,0,p)= 242 iE 


Problème aux limites de Laplace — système de coordonnées toroïdales — p 659 
R-séparabilité 2 : problème extérieur dépendant des variables n et ọ : 


D'après le tableau des solutions finies la solution s'écrit pour une condition aux limites paire en Ÿ et 
@ : 


AAEE n en 2 2 AmQ, 1 (Cotanh(n)})Cos(n 2207) 














T(1,0,6),, = (0,0) = Fee E È 4,0" :(Cotanh(l,)Cos(n9)Cos(mp) 


Sinh(l ) n=0,+00 m=0,+00 


À fas Cos(nd9) 
"a  JCosh(,)-Cos(O) % 








í do f (0,p) Cos(mp) 





” i(Cotanh 
TE a, Q (Cotnh(n)) 
me : — Cos(n9)Cos(mop) 
Sinh) | ne nee (1+8,0)1+ 8,9) Q” :(Cotanh(l,)) 
2 


baie JCosh(D= Cos(®) J 


T 








Si la fonction limite ne dépend que de ọ, il vient : 


Cos(n3) 2x 
Jas JCosh(l, )—Cos(0) Cu = [de ro Cos(mop) 


Cos(nà9) 
JCosh(1,)-Cos(0) 


Anm =B,Cn > B, = 











Or B, = fao =2V20 1 (Cosha, ) 





o (Cosh, ), 0"; (Cotanh(n)) 
FRIC) 











T(n,0,@)= 24/2 CO- Cos(0) er n) 
z ShO) | 22. 2. + BA On) o" (Coran) 


avec la formule de Whipple 


1 
rfnim+s) 2 2 
P", (Cosh(l,))= - | | Q" (Cotanh(l,)) 


m(n- PE 2) r Sinh(l,) m- 
2 
r|» 7 ;) O ,(Cosh(l,)) C,Q" _,(Cotanh(n)) 


Cosh(nņ)- Cos(0) 5 2 2 Cos(nļ)Cos(mę 
nS) | nono en ss ) P"i(Cosh(,)) (+6,,)1+6,,0) 
2 2 














T(n,0,p)= S] 


En appliquant une fois de plus la formule de Whipple sur le terme : 





r|» +m+ ) 4 

2 2 

P", (Cosh(n))= | de ) 9": (Cotanh(n)) 
2 rfa- Pn 2) m Sinh(n) 3 


Il vient : 
Q, (Cosh(,)) C, P” 1 (Cosk@n)) 





T(.0,0)= © 





Cos(n9)Cos(mop) 





De p Cona, )) (+6, rS, 5 


Ce qui est exactement la formule obtenue avec la ! 
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La solution pour une condition aux limites impaire en È et ọ se développerait comme suit 


T(n,0,p)= eee >. 2. À, 20% , Coranhtp)Sinn9 Sinne) 














T(1.0,6), 6.0) RD O > È 4,0" (Cotanh(,)JSin(n9)Sin(me) 


Sinh(l ) n=0,+00 m=0,+00 

















Sin(n9) 
À, „= de f,(6,p)S 
” | Tert )- Í P fa (0,9) Sin(mp) 
- Q” ı(Cotanh(n)) 
_ „Cosh(n)-Cos(0) it Apm m1 | 
ee T° SAOD | one Oro 8m) O (Coma GC) 


Si la fonction limite ne dépend que de ọ, il vient : 


Sin(n9) G = [de fa (o) Sin(mp) 


Anm = B,Cn > B, = fas JCosh(,)- Cos(0) 











T(7.6,0) = TUE Cos(0) 








Et i oi (Cotanh(n)) 


+ B Cn — Sin(n9)Sin(mop) 
Sinh) |, ne A +6,0)1+ 8,0) Q", (Cotanh(l,)) 
2 


Comme toujours il y a une totale équivalence des formules de R-séparabilité-1 et R-séparabilité-2. 
La formule de Whipple combinée comme suit: 


P" (Cosh) Fe TA Q" (Cotanh{m) 


P” (Cosh) Sinh(n) Q” Ca) 








Il vient : 











o" ı (Coranht) 
VCosh() - Cos(0) [e Anm 2 jes 9Cos(me) 


T,0,p)= 
CRO m SIAOD | none CFE, + 8,6) O" OA 


1444444444444444444444244444444%8#4444444444443 
U 





P” (Coshn)) 


ĮJCosh(n)-— Cos(0) A,n -1 
T(n,0,9) = : 2 C 9)C 
GP) m 22, UE i) P" (Cash) PEE E 
"3 











o" (Corah) 
di jS 9)Sin(mp) 








_ VCosh(n)-Cos(0) |Sinh(l,) Anm 
EANO T° SA 2.22. (+6,6)1+6,,0) 2” | (Cotanh(,) 


1444444444444444444444244444444444444444444 943 
U 





P” ,(Cosh(n)) 


JCosh(n) - Cos(0) À,» > , i 
T(n,0,9) = > > > 2 S Q) S 
CERO z? nom A+ 8,9) + 8,9) P”, (Cosh(l,)) SPORE) 
F2 
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Remarques sur les fonctions propres sur des tores complets pleins en fonction de la 
décomposition par symétrie des conditions aux limites 


Dans la présentation de quelques unes des solutions des problèmes aux limites inhomogènes de 
Dirichlet ou de Neumann (non traité ici), nous n'avons pas introduit la décomposition des fonctions 
limites sous forme de symétrique dans le plan polaire (soit pour un angle azimutal @ constant). 
Toutes les solutions présentées, pour autant que les calculs soient justes, sont évidemment dans 
une des catégories de symétrie Y+X+, du seul fait de la construction des fonctions propres 
angulaires à l'aide des fonctions sinus et cosinus. 


On décomposerait par exemple une fonction limite en ces quatre composantes de symétrie, dans 
un problème sans dépendance en @ (utilisation de la R-Séparabilité-1) 


fO) = frs: (0) + frs x -(O)+ frx 0) + fr x_(O), sachant que f(0)= f(2x +0) 
f (0) -ZOI -0+ fO) + F6 -7) 











4 

fra 0) POLE D ICO PEN 
fx. (0)= f(0)- f(x —-0) ee _ f(8-x) 
fr = OSSI EE?) 


On vérifie aisément que les fonctions construites vérifient bien les propriétés de symétrie : 
Jrxx O= fix rt —-0) et fr, x:(0)= fi, (0) 

Jrix-(O)= fix, 070) et frx O) = — fr, xO) 

fr-,x:(0) = -=fr x+ (z-0) et Jrs xa (0)= fr, x0) 

Sr- x-O)=-Sr x-0) et fr, x0) =-fr x0) 


Le tableau suivant indique les éléments pour la construction d'une solution par principe de 
superposition utilisant les fonctions de la R-Séparabilité-1 : 




















Symétrie Valeur |Fonctions Normes des Fonctions radiales 

angulaire | propres | propres fonctions propres R-Séparabilité 1 

en Ÿ À, angulaires et) angulaires 

normes ®,(8) | ||®,(8)] 1? 

Y+X+ 2n Cos(2n8) n=0->2n, n>0->Tt P  (Cosh(n)) Q, (Cosh(n)) 
"5 Re | 

Y+X- 2n+1  [Sin((2n+1)8) n E i (Cosh(n)) o (Cosh(n)) 
2 Fa : 

Y-X+ 2n+1 |Cos((2n+1)8) m P  (Cosh(n)) o (Cosh(n)) 
niz de , 

y-X- 2n Sin(2n8) T p  (Cosh(n)) O, (Cosh(n)) 
Lors n-5 l 
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On pourrait compliquer à plaisir la construction en introduisant la décomposition symétrique 
similaire en angle azimutal dans le cas où la fonction limite présenterait une double dépendance en 
Ô et @, on obtiendrait 16 composantes en tout en introduisant les fonctions radiales de toroïdales 
associées d'ordre m supérieur, m étant la valeur propre azimutal utilisée. 


Cependant il peut être utile d'exhiber la construction des fonctions en R-Séparabilité 2 avec le 
tableau correspondant pour une fonction limite qui n'a pas de dépendance en ©. 


SO) = frx lO) t frix (PO) + frx Ot frx), sachantque  f(p)= f (27 +) 
je f(e)+ f(x -p)+f(p)+f(e-7) 














Sraz O A 
po O EDE 
ad O O 
pa Oe 


Par ailleurs comme le développement en série de la solution implique nécessairement un terme 
dépendant de venant du terme de R-Séparation-2 suivant, la dépendance finale de la solution en 
Ÿ et ọ est inévitable : 


T(n,6, p), = f@)= YY H(4,)9,(0)Y, (p) = A eL R(7,0) = JE EE) 











L'emploi des fonctions angulaire en Ÿ est donc obligatoire, mais restreint à la seule valeur À,=2n, 
car la fonction limite est par essence de symétrie Y+X+ en Ÿ (fonction constante à l'égard de Ÿ). Le 
tableau suivant indique les éléments pour la construction d'une solution par principe de 
superposition utilisant les fonctions de la R-Séparabilité-2 : 
































Symétrie | Valeurs | Fonctions propres | Normes des | Fonctions radiales 
angulaire | propres |angulaires ®,(0)0,(p) fonctions R-Séparabilité 2 
en @ An Àm propres 
angulaires 
|1Onlp)] 1? 
Y+X+ 2n,2m |Cos(2n8)Cos(2mvp) n=0->2n, pe  (Cotanh(n)) o”  (Cotanh(n)) 
n>0->Tt 2 
m=0->27, 
m>O->Tt 
Y+X- 2n,2m+ |Cos(2n8)Sin((2m+1)æ) |n=0->27, P” \(Cotanh(n)) O?” ,(Cotanh(n)) 
1 n>0->Tt 2 2 
T 
Y-X+ 2n,2m+ |Cos(2n8)Cos{((2m+1)p) |n=0->2n, | P” (Cotanh(n)) Q” ,(Cotanh(n)) 
1 n>0->n 3 Fe 
T 
Y-X- 2n,2m |Cos(2n8)Sin(2mp) n=0->2n, P™  (Cotanh(n)) Q” ,(Cotanh(n)) 
n>0->7n ee 72 
T 
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Remarque sur la solution du problème aux limites inhomogène de Neumann ou mixte de Robin 


Concluons ce chapitre sur les problèmes aux limites sur un tore complet avec des conditions aux 
limites inhomogène de Neumann ou mixtes de Robin à la surface du tore. Le calcul n'est pas 
analytiquement plus complexe à représenter sachant qu'il faut tenir compte de l'expression du 
gradient en coordonnées toroïdales : 

Gaata te CAC Coste) OAR PTE). E 07620) 

c "an 00 ? Sinh(n)  ôọ 
y Cosh(l,)— Cos(0) ÔT(n,0,@) 
c Ô 








= f(8,e)= 


n=l; 


Robin & + BT(n,0,0) 


OT(n,0, 
ILE + BT(,.0.0) 








n=l 


ôTN,0,p) 


On 


Cc 


Neumann = 
Cosh(l,)— Cos(0) 








€ OT(n,0,p) 
= f(6,p)= : 





f(0.p) 





n=l; 
JCosh(n) —Cos(0) R- séparabilité — 1 


R(7,0) = JE — Cos(0) 
Sinh(n) 


n=l; 











R — séparabilité — 2 











Cosh(l,)- Cos(0) dH(,) dR(I,,0) 
rins yy? (2a, an S a i 6,0)4, (p)= f (0.9) 
” o" |+ BR(,,0)H(1,) 
dH (l) dR(l,,0) ____cf(6.p) 
Neumann = ZE 20,0) dh + aU lo, (8)T,,(p) = Cosh(l,) = Cos() 


Cela complique évidemment les calculs, mais ce n'est pas infaisable pour un problème en 
particulier. 


= f(0,œ) 
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Solutions des problèmes aux limites homogènes angulaires 9 sur des sections « coniques » 
toroïdales 


Jusqu'à présent nous n'avons aborder que les problèmes aux limites sur des domaines où les 
angles 6 et b balayent tout l'intervalle [0,27] (tore complet). Considérons d'abord une restriction 
de valeurs sur l'angle 6. Dans ce cas la frontière du domaine est constituée par les iso-surfaces 
angulaires 8=Cste et une au moins des iso-surfaces radiales n=Cste. Si l'intervalle radial est complet 
n€[0,°] alors la solution du problème de Dirichlet homogène est trivialement nul (domaine fermé, 
simplement connexe, principe du maximum et du minimum). Il en est de même pour une 
problème homogène de Neumann ou un problème mixte Neumann, Dirichlet sur l'une ou l'autre 
des iso-surfaces 8=Cste. 


Soit à rechercher la solution de l'équation de Laplace à l'intérieur d'un domaine défini comme suit : 
= (70,p)/n E IMAC s (9,0.1e z [0.27] Et sur ce type de domaine on considère des 


conditions homogènes sur les deux iso-surfaces angulaire 6. De ce fait on reste dans la même 
catégorie de fonctions propres angulaires utilisées, les sinus et consimus mais avec des valeurs 
propres discrètes réelles. 


Il est important de noter que l'intervalle en coordonnée radiale a des formes limites : 


Ll ou |l ,+œ| ou |O,l : : ; 
ly. p2 ly» | l i rol, Dans ces trois cas on peut assurer par construction séparable 


la finitude de la solution dans le domaine, ce qui n'est pas le cas pour l'intervalle complet 
n€[0,+2°]. Dans ce dernier cas à la section suivante, comme cas particuliers, on aborde la 
construction par des transformations intégrales dont les noyaux sont justement issue des fonctions 
propres de séparation. Dans la section suivante les autres cas abordent des conditions homogènes 
radiales et non angulaire (on passe alors des fonctions radiales toroïdales aux fonctions propres 
coniques de Mehler). 


Quelques illustrations géométriques des domaines à l'étude 





Par exemple pour Q = (8.o) € [1,1,0 ē E | z bo € bar] 





1.0 ii 








0.0 0.5 
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o= {non eoo e| $E] pelz] 
Par exemple pour 4 4 





Q= (0.o) efl,o}0 e Bag c bar] 


Par exemple pour 
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2 


Q= (0.o) e[0.5,218 € Ë Tlo € bar] 
Par exemple pour 4 





Les conditions aux limites de Dirichlet ou de Neumann sont homogènes en angle et fixées 
radialement en paroi intérieure et extérieure à des valeur différentes : 


AT(N,0)=0 (1,0) €Q =fr e |La1,b0 € [0,0] 








TO) ga TOY g] 
na T HUE 
rao ST? orao] | 
ou —T (n, ou —T (n, 
00 i nel; 1: >| 00 g 


nelinin2 ] 
0=02 


= f,(0) 


0-[h.02] 


TO) E so) T), 2 


0h02] 


Pour des raisons de simplification de l'exposé, on se restreint aux conditions homogènes de 
Dirichlet sur les surfaces à distance fini et sans dépendance à l'angle azimutale ©. 
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Dans ce type de problème les solutions peuvent se développer sous forme de série, mais avec des 
valeurs propres angulaires non entières. 


Solutions sur les domaines radiaux n€[l,o,+°°] 
AT(n,0)=0 (n0)eQ=hell,,;+0]0 €[6,,6, |; 


TO TD ,=0 


adno: +] 
0=0 





T0) = f(0) Tr) finie 


sa 2] 
La solution se développe en série sous la forme suivante, compte tenu des conditions aux limites 
homogènes et de la condition de finitude : 


n rA 2 NROS /Cosh(n) - Cos (0) 2 49, ne CS) 


l'application des conditions aux limites nee ne les calculs suivants : 






































T0), , = f(6)= JCosh(,5)- Cos(6) > B,Sin(, 90, a (Cosh@,) 
S 2 _0,-6, o2 Po Sin(2,9) 
OO CUITE rl [JCosh(l a) - Cos(0) 
Q (Coshn)) 
__2,/Cosh) — Cos(0) al | f (0) Sin(A, (6 —-0,) 
T(n.0) = Re 24 0 AC) de (8-8) 4, - fas JCosh(,,)- Cos (0) 


Solutions sur les domaines radiaux n€[0,l,o] 


AT(n,0)=0 (n,0)eQ= in € 0,2, b0 € [00,]} 
TONON iz =0 T(7.0) 


neoo] 


0=02 


T(,0),.. =f0) T(7.0)] finie 
La a développe en série avec les fonctions toroïdales de première espèce: 
ra er T(n,0) = JCosh(m) - Cos@) D AP P 1 (Cosh@m))Sin(a, (8-0,) 
l'application des conditions aux tes homogène donne les calculs suivants : 
Ar 2 
2 À, P dean at, (8-8) 4 = fas cet 
+ 2 














2.,/Cosh(n) - Cos(0) 


T(n,0) = 
(7,0) 0,0 
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Solutions sur les domaines radiaux n€fl,:,l,2] 


AT(n,0)=0 (n60)ea=hel,,1,,/0 [6,0] 
T(7,0) T(7,0) 








zeloo] = 
= fi(0) T0) 


de Je @] 


yeeo] 


T0) = f,(0) 


de Je @] 


La solution se développe en série avec les fonctions toroïdales de première et deuxième espèce: 





nT 


T(n,0) = /Cosh(n) — Cos (0) x 9—0 


À, = 





afe  ı(Cosh(n))o. i (Cosh(l,2))- 


a sfe 1 (Cosh(n)o, (Cosh ))- 


Sin(2,(0 -0,)) 


P (Cosh(, DO, CAC) )} 


3 i (Cosh(, DO, : 1 (Cosh@n) ) 








f Fe a (CosaG,))0, i (Cosh(l,»))- 
s 1-3 


P F (Cosh(l,, 0, i 





ı(Cosh(l,)) 


L'application des conditions aux limites inhomogènes donne le résultat suivant : 





T(n,0) = JCosh(n)-Cos(0) X. 


n=1,+00 


ae a (Cosh(m)o, (Cosh(l,:))- P 


-a[r (Cosh) 


À 


(Cosh, ia em) 


(CosA(L,,))- P (Cosh, ED) 


1 
Eg 


Sin(2,(0 -0 ) 








x 
P 
An 


g 1O Sin(2,(0-0,) 





(Cosh), | 


(Cost, )}P P :(Cosh,)0, (Cosh) 


2 











4 -f 4 JCosh(l,)- 2) 


nl 


B, = fao FOR „(0 -0,) a, = "7 
JCosh(1,,)- Cos(0) 


0, a 0, 
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Solutions des problèmes aux limites inhomogènes angulaires 9 sur des sections « coniques » 
toroïdales, conditions homogène radiales 


Il s'agit d'aborder des problèmes sur des sections coniques (soit des coupoles sphériques) coupées 
par une au moins des iso-surfaces radiales. La géométrie des domaines est identique au cas 
précédent, seul l'homogénéité des conditions aux limites change. Cela entraîne le choix des 
fonctions sinus et cosinus hyperboliques angulaires pour respecter l'inhomogénéité des conditions 
aux limites angulaires. Cela entraine l'utilisation des fonctions coniques de Mehler en dimension 
radiale. La question est donc : peut-on construire un système de valeurs propres uniquement 
discrètes et de fonctions propres de type conique de Mehler respectant les conditions aux limites 
radiales homogènes imposées sur l'une au moins des iso-surfaces radiales. Cette discussion est tout 
à fait similaire au cas des coordonnées sphéroïdales allongées ou aplaties. 


Les conditions aux limites de Dirichlet ou de Neumann sont homogènes radialement et 
inhomogènes en angle à des valeur différentes. Pour des raisons de simplification de l'exposé, on 
se restreint aux conditions homogènes ou inhomogènes de Dirichlet sur les surfaces à distance fini 
et sans dépendance à l'angle azimutale œ. Trois problèmes peuvent être envisagés : 








AT(n,0)=0 
nel] 
(7,0)eQ ee | CRC Le apea= [neo] 
T(n,8) Si 0 e[6,,6, | 0 €[0,,6,] 
acao] T(n,0) = 0 T(n,0) =0 
T(n,0)  =0 ou otha ou Phs 

adaon] T(n,0) es = fi(n) T(n,0) ae = fi(n) 
TO = An) a e 

ap T0) =A) T.O) =A) 
T0) = fln) R ma 








DS 


Eli 2] 
=02 


Les solutions par séparation sont construites à partir des fonctions suivantes : 


T(n,0)= /Cosh(n)- CO} 4 Pi, (Cosh(n))+ B 01. (Costa) fe Cosh(r0)+ D Sinh(r@)) 





iT iT 


ou encore 





T(n,0)= JCosh(n)- QUE P, (Cosh(n))+ B Qi. (cantn) fe Sinh(r(0 — 6,))+ D Sinh(r(0, -0))) 


iT 3 iT 
La finitude des solutions est gouvernée par les propriétés asymptotiques des fonctions coniques de 


Mehler), à savoir notamment la valeur en z=1, et z=% . Ces deux cas limites ont la conséquence 
cruciale d'éliminer la condition aux limites homogènes correspondantes : 


Log(z — 1) -o A. Cv t 1) 


gels- Srv TED 

y Constante d'Euler w(v) fonction Digamma 

= LimQ, (Cosh(n))=x et P, (1)=1 LimQ,(z)=0 et LimP, (z)=0 
nl ont FAT Z—>+00 i 








+ O(z 1) zel 
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Discussion sur le cas de l'intervalle n€[0,1l,.], z=Cosh(n)->z€[1,z,] (voir système sphéroïdal): 





La condition aux limites homogène en z=1 n'est plus de mise, remplacée par la condition de 
finitude. On peut donc écarter les fonctions de coniques de Mehler de deuxième espèce 
divergente en z=1, il reste à construire les solutions avec celles de première espèce. On a vu dans 
le système sphéroïdal qu'il s'agissait d'un problème de Sturm-Liouville singulier mais qu'il possédait 
les bonnes propriétés pour assurer l'existence d'un système de fonctions propres complet et 
orthogonal. Par construction on peut vérifier numériquement que les solutions de l'équation 
transcendantale suivante existent : T, tel que P T (z,)=0 z,= Coshl,.) 
LL 


Par exemple pour une valeur z0=Cosh(1), il vient les valeurs propres suivantes , t,=2.42138, 
5.52719, 8.65826, 11.7949, 14.9335, 18.0732, 21.2135, 24.3541, 27.4949, 30.6359, … 


On peut illustrer graphiquement l'existence des solutions de l'équation transcendantale en traçant 
la fonction conique de Mehler de première espèce en fonction du paramètre T. 


| | LIT IlLIlI | | | | | | | | | | ||. Re(P°, , (cosh(no)) 























Discussion sur le cas de l'intervalle n€[l,_ln2], z=Cosh(n) -> z€[z1=Cosh(l,1),z2=Cosh(l,,)] (voir 
système sphéroïdal) : 


Dans ce cas le problème de Sturm-Liouville devient régulier avec des valeurs propres négatives 
(entraînant des paramètres imaginaires pour les fonctions propres), il existe un spectre discret de 
valeurs propres et un système orthogonal et complet de fonctions propres . Les valeurs propres 
peut être construites par la résolution d'une équation transcendantale comme suit pour un 
problème homogène de Dirichlet en ly: et l,2: 


Oo ordo, Œ) T, telqueP; (Coshli,, Joe. (Coshll,2))-P ; (Cosh, Jo"! (Coshll,) 


E ——+it, tn To tn To tn 


On peut démontrer que l'expression de l'équation transcendantale en utilisation la fonction de 
Legendre Q est identique en raison de sa réalité : 
réel réel 
VX, X2 ous GO. G)-P 1. (eot (x)= P (x )2 l (xa)-P i (2) i (x)eR 
2 E 2 


n n n n x n n e n 


>T, tel que a (Cosal)o i (CosAlt, ; ))- Sli (Cosal,: Jo Y (Coshltn ))= 0 


—+iT ——+iT 
n n n 
2 2 
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Un autre choix possible de fonctions à valeur réelle indépendantes est donné par : 


e) (cohinh-0 , (Cosita) 


n n 


Pi, (Cosh(n)) = 


n 


n n 


Cosh\rt 
L, Contre Po, (Costa), (Costa) 
2 2 >” 2 
Voici le graphe de l'équation transcendantale en tant que fonction du paramètre t, les solutions 
étant l'ensemble des points d'intersection avec l'axe des abscisses : 


List = 


Fnoh)e 8 1, (Cosh) : - (CosAf)-P, (Cosh O ;  (CosAli,;) 


=g a =a ta 





[\ {\ Æ / \ 7 P, „ (cosh(ha)) Ref, (cosh(in)))-P°, , (eoshfin))Re( a°, (coin) 
VAV 








Pour l'exemple, avec les valeurs ] =1 Z ,=2Voici le jeu de valeurs propres Tn : 
n n 
Tn =3.1319869898174106382, 6.2779019008125322845, 9.4211834216040919190, 12.563654676027245805, 


15.705782929423131877, 18.847735498245959850,21.989586404634429313, 25.131373298661280558, 
28.273117314492634571,31.414831218507219320, 34.556523172690331259, 


Et avec les valeurs 1,:=4 l,2=6 voici le jeu de valeurs propres tn : 1.57078, 3.14158, 4.71238, 6.28318, 7.85416, 9.42466, 
10.9407 


Discussion sur le cas de l'intervalle _n€[l ọ+°°],_z=Cosh(n)->z€[z0=Cosh(l o), +22] (voir système 
sphéroïdal) : 


Cette fois-ci, les deux fonctions de coniques de Mehler de première espèce et deuxième espèce 

s'annulent en z=+c. Dans ces conditions toutes les valeurs de t sont licites pour construire une 

solution. Tout se passe donc comme si les valeurs propres possédaient un spectre continu. Il suffit 

de choisir la fonction suivante pour respecter la condition aux limites homogène en n=l et la 

condition à l'infini: 

D(7)=P 1. (Cosh(n)O : 
2 


Eu 


(CosA(l,0 ))- EN (Cosh(n))? ;  (Coshll o) = ®,(40)= 0 


iT ——+iT 
n n n 
2 


Dans ce cas il n'est pas possible de construire un système discret et complet de fonctions propres 
ayant les bonnes propriétés sur les deux bornes de l'intervalle. En s'inspirant de la forme des 
solutions préconisées par la séparation des variables, on développe une intégrale proche de la 
formulation de la transformation intégrale de Mehler-Fock. 
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Discussion sur le cas de l'intervalle n€[0, +22] (voir système sphéroïdal) : 


Même si ce cas n'est abordé dans ce chapitre, il n'est pas possible de construire un système discret 
et complet de fonctions propres ayant les bonnes propriétés sur les deux bornes de l'intervalle. On 
fait donc appel à la transformation intégrale de Mehler-Fock pour construire ces solutions. 


Problème aux limites sur l'intervalle [0,l,o] 


Revenons aux cas d'un intervalle fini du problème aux limites. La solution d'un problème aux limites 
de Dirichlet de la forme : 
AT(n,0)=0 (n.8)e0<n<l,,8 <0<6,} Lim, ,T(n.0) finie 

n—0 


T(n.0),,,=0 TO = fr) TO, =A) 


peut se développer comme suit, en intégrant au calcul les normes des fonctions propres radiales 
déjà été calculées dans le cas du système sphéroïdal, sur des problèmes aux limites de Dirichlet 
dans des sections d Se à deux A 


T(n,0)= /Cosh(n)- Cos(0 12 a _(Cosh( n(A, Sinh(r, (0, -0))+ B, Sinh(r,(0 -0, ))) 
























































f(n) 
,0 h( h(t (6, -8,))) = 
T sr © À pe (Cosh(n)\4, Sinh(r, (8, -0,))) Tin) CG) 
' nm Sinhln)A(n)P , _ (Cosh(n) 
>Á = dn 2 
lo, (nf Sina(r, (0, -0,)) % VCosh(n)-Cos(8,) 
| n Sinh(n)£(n)P,. (Cosh(n)) 
FX > 0) =f, >B, = d 2— 
PO) gyra TA jo, (nf simnl, 0-0) VCosain)- Costo) 
le, m) = a- air) Ce (Cosh(s,)P (Coshll ) 
TP (Cosh(n)) 
= T(n,9)=4/Coshfn)-Cos(8) X, {(1-2ir,)P, (Coshir) ZP, - (Cosh(l o )SinA(r, (8, -6,)) 
n#0,+00 Eu T ne 
x(4, Sinh(r, (0, -0))+ 4, Sinh(r, (0 —-6,))) 
n Sinh(n)A(n)P 1. (Cosh(n)) Sn) (1). (Cosk(n) 
AT a a LT on au 


Prenons le cas où les deux fonctions limites sont constantes, la solution du problème se développe alors comme suit : 


T,P ı „ (Cosh(n)) 
g a 





x 





n#0,+00 


T(n,9)=44/Cosh{n)-Cos(e) 5 (i-2ir,)P, (cos) P, (Coshll,o ))SinA(r, (6, -0,)) 


x(4,, Sinh(r,(0, -0))+ 4, Sinh(r,(0 -0,)) 


ln 


$ osh\l o P osano P | 
= À, =T m 7 nr . et À,,=T. | f i D 2 
Si fa(n)=T, aia 1 Jz-Cos(0,) UE A z-Cos(6,) 
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Problème aux limites sur l'intervalle [l,;,1,] 


La solution d'un problème aux limites de Dirichlet de la forme : 
AT(n,0)=0 (n,0)el, <n<l,,,0, <0<6,} 


n 
TOO), =O TOO a, = TOO = fn) T, =A) 
est alors la suivante : 
Pi (Cosh(n)) Q (Cosh(n)) 
itn 


12e 
Hit, 


Fa (CosAl,,) ” | (CosA( An j P 
2 —+HiT 


n 





D. (7)= 





m) = Jan sinnn) (no, (n)} 
a. fan Sinh(n) f(n). (7 a. fan Sinh(n)f.(n}.. (n) 


1 Cosh n)- 1 Toha )- Cos(0,) 
T(n,0)= JCoshln)-Cos@) Y. (4, Sinn(e (8, -0))+ 4, Sinh(e, (0-8) à (n) 


De a) Sinh(t,(0, —0,)) 


La norme des fonctions a déjà été calculée dans le cas de l'hyperboloïde à deux nappes, il vient : 
= Cosh(l,,) Z, = Cosh(l,.) 


0,0%, rdo, ©) o O-o, 6) 




















n#0,+00 | 
































0,0%, rdo, ©) 0%, O-o., 6) 
Eu 3 tn soi mii =a ia 
CCI T Oa a 
Tn \N2 DAA = 7 pride, = 2 = qu = 2 
dr A ie (z1) Oi (z) atii (z) Qi 2) HIT (£) 
2 1 Da 2 2 
lo, (7) De: | ” 
À a ta O BEN [| O o O pee e) 
Faid 2 = +7,92 2 = 
dr 2 Fire 1) Oie (z) Oiz (z) Dis (z) ——+iT ( 1) 
2 2 i 2 
ôP La (z) 2, (z) aP, (z) 60 Eog (a) 
jee 0®., LE dr z dr 0®., Gi dr _ dt 
dr cie (a) Of. 2) dr a (z,) re (zi) 
2 i 27 
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Ce qui donne formellement la solution du problème aux limites pour une fonction limite constante, 
comme suit : 


P, _(Cosh(n)) O, (Cosh(n)) 


Pi M (z,) O1, (2) 
ztn ztn 
T, telque Pir 0. (2)=0:. LAS A (2) 
7 En 77 LL 


Z= Cosh(1,,) Z, = Coshl!,) D, (n)= 





2, O-o, rdo. ( 


75 ia 








0.00%, rdo, O) o e-o, () 
) 
) 






































aD, (la) 1 a N TE 6) pe 
Tn \ M2 2 _2 trader i 2 +— 
dT 2 Pas (2) Q (z, Qi (z,) 0, (2) Pi. (2) 
2 1 2 HT, tit -gtin z 
je., r) i 2t D D imag imag imag 
| op ( ) 1 Paie () Os. (a) CEA (2) je (a) re (a) 
Ta \71 D _2 +T, Z = E 2 ï p2 
dr 2 Aa () Gus (a) he (a) PATES (z) “he (a) 
Pi 2 $ 
ôP (2) oQ la (z) oP (2) 0Q be (2) 
Ta 
Fe ôP., Ca) dr dr 0®, (i dr _ dr 
dt Pi, (2) Qi (2) dt P Ey (z,) Ori ( 1) 
B C B C 
f(n)=T. Aa) T = À =T In — In A, =T. 2n = 2n 
i | . Fire. ( 1) Qı N (z,) ? : Pi (z,) Cas ( k 
2 í 2 
2 de a 7 Qu. (2) he ( ) 2 o L +iT () 
B,,= |d -=== C, = |d -— | 


——— d 
A Jz- Cos(0,) 7 Jz-Cos(0,) = WE Tos, ) fe AZS TE ) 





























B C 
T, Lu = Sinh(r,(0, -0))+ 
| Pi, (a) Gi (a) l 6. ) 
2 2 
PC) 
T(7,0)= JCoshfn)- Cos@) 5 1. Gi) 0, @ 
n#0,+00 i 2 x 
KOG J Sinh(r, (6, -0,)) 
P, (Cosh) à, (Coshn) 
x 2i 
Py (a) Qi, (a) 
Z 
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Problème aux limites sur l'intervalle [l,5,+°] 


Comme on l'a vu dans la discussion précédente, la solution d'un problème aux limites de Dirichlet 
de la forme : 


AT(n,0)=0 (n,0)e{0<l p <n<+0,0, <0<06,} Lim, 5T(n,0) finie 
n—0 


TAO) a, =0 T0). = A) TO, = A0) 


est construite à partir d'une représentation intégrale de type de Mehler-Fock (voir plus loin). 
Problème aux limites sur l'intervalle n€[0, +22], z=Cosh(n)->z£[z0=1,+c°] 


La solution d'un problème aux limites de Dirichlet de la forme : 
AT(n,0)=0 (n,0)e{0<m<+o , 0 <0<0,} Lim, ,T(n,0) finie 


ou 
7x 
ou 
n—0 


T0) = fn) 


est construite à partir d'une représentation intégrale de type de Mehler-Fock (voir plus loin). 
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Solutions des problèmes aux limites sur des domaines union, intersection et différence de 
sphères, lentilles sphériques 


Le système de coordonnées toroïdales a la particularité de pouvoir modéliser des opérations 
géométriques particulières sur les sphères notamment les unions, intersections et différences de 
sphères, en utilisant les iso-surfaces Ÿ=Cste. Dans la suite de l'exposé on utilise les coordonnées Ÿ 
cette fois-ci dans l'intervalle [0,27] plutôt que dans l'intervalle [-n,x}, car la manipulation d'angle 
négatif introduit dans les raisonnements des intervalles non contigües moins commodes à 
manipuler. On se conforme également en cela à la présentation des solutions des problèmes aux 
limites telles que celles de l'ouvrage de N.Lebedev : « Special Functions and their applications » 


Comme indiquer en tête de chapitre sur les coordonnées toroïdales, on remarque que pour les Ÿ 
croissants positifs les iso-surfaces sont successivement englobées par la précédente. En partant 
pour la valeur 9=0 d'une sphère centrée à +2 de rayon infini (à l'exclusion du disque de rayon a 
centré à l'origine dans le plan (x,y)) jusqu'à devenir ce même disque de rayon a pour la valeur Ÿ=rr. 
Pour les Ÿ > n les iso-surfaces cette fois englobent successivement les précédentes jusqu'à devenir 
une sphère de centre tendant vers -> et de rayon infini, à l'exclusion du disque de rayon centré à 
l'origine dans le plan (x,y) pour la valeur ÿ=27. 


Il apparait donc pour une iso-surface Ÿ= Ÿ; d,€[0,n], que la zone convexe qualifiée d'intérieure est 
VE[Ÿ. ,n] et par opposition la zone concave, dites extérieure est ŸE[n,2r+0.]. Pour une iso-surface 
= Ü1, 0,€[0,n] la zone convexe qualifiée d'intérieure est ŸE[n,Ÿ] et par opposition la zone concave, 
dîtes extérieure est ÜE[D:,2n+0:] 


Modélisation d'intersections de sphères (lentilles sphérique) 


Les domaines définis par les intervalles de coordonnées toroïdales suivant, construisent des 
intersection de sphères, que l'on peut plus communément dénommées lentilles sphériques 


5D = T for asbsr)uossxé 08) 4 (a gs ET a Ka 
Le domaine intérieur est défini comme suit : 

Dyu = CEE [0,+]x[9,3]x[0.27] tą (a e [Zz] et 9, € [2x ])} 

Le domaine extérieur est défini comme suit : 

Dye fons.o) e [0 +0]x [3,27 +3]<[0.27] tą (a è [Za] be [2x ])} 
Exemple de lentille sphérique avec : 


4 


=T—@; 
9 = 9, =7<e = 
4 





4 4 & = 
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Modélisation d'unions de sphères 
Les domaines définis par les intervalles de coordonnées toroïdales suivant construisent des unions 


oz] et Q, € [z221] 





de sphères : 
(n,92,p) = [0,+0]x 9 x[0,27]L[0,+0]x 9 x[0,27] tą (a E 


ôD = 
ou (a E [o.z] et 9, € ar) 


Le domaine intérieur est définit comme suit : 
(n7,9,p)e[0,+0]x[8,8,]x[0,27] tą (a € [oz] et Q e Irar]) 





Dix = 7 
ou (a ef0,r]er & € É ) 


Et le domaine extérieur est définit ainsi : 
(n,3,p) e [0,+0]x [9,2m + 8,]x[0,27] tą (a € [oz] et Q € War) 


Ds = 
ou (a ef[o,rler 3 € er) 


Exemple d'union de sphères avec 





Exemple d'union de sphères avec : 
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Modélisation de différences de sphères 


Les domaines définis par les intervalles de coordonnées toroïdales suivant construisent des 
différences de sphères : 


> — ho e([0,+01]x 9 x[0,27]) o ([0,+0]x 9 x[0,27]) t4 <93 et (9 el0,rler 8, [0x] l 
ou (9, efx,2r]et 9, efx,2x|) 
Dans les deux cas, le domaine intérieur est définit comme suit : 
na o [0o +o]x [9,9 ]x[0,27] t4 <9 et (9 =[0,r]et 9, e[0,x]) | 
M lou (9 efx,2rler 8, efx,2x|) 
Et le domaine extérieure comme suit : 
À — €[0,+]x [9,27 +89 [x[0,27| 19 G<3 et (9 efl0,rler 8, ef0,x)) | 
~ lou (9 efr,2r]et 9, efr,2x]) 
Exemple de différences de sphères avec 


+7 





Exemple : solution du problème de Dirichlet sur le domaine intérieur de la différence de deux 
sphères 


L'exemple suivant résout le problème de Dirichlet dans l'espace intérieure d'une différence de 
sphères. Il est donné par à N.N Lebedev dans son ouvrage « Special Functions and their 
applications », page 227, section 8.12. 

AT(n,9,p)=0 (n,9,p)e D, 


Tn, LP); = fine) T(n,8,p),., = (1.6) 


D= o [0o +o]x [9,9 ]x[02z] tq 9 <9 er (9 e[0,xler 9 e[0,x]) l 
M lou (9 efx,2xl]er 9, ef[x,2x)) 


Rappelons que les solutions d'un tel problème de Dirichlet inhomogène en à, sont de la forme : 





ne I [Cosh -Cos 4 P% (Cosh(n))+ B 0‘. (Cosn(r)) (4Coshlz)+ 8 Sinh(t0))x 
n, >P = Eo Zoe 


x (A Cos(m p)+ B Sin(mœp)) 
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Supposons que le problème ne dépendent plus de l'angle azimutale ®, alors la solution se construit 
à partir d'une superposition des fonctions suivantes : 





U(7,0) = Cosh(n) E CO) 4 Pi. (Cosh(n)+B0Q, | (Cas) 4cosi(r0) + B Sinh(10)) 
=i Forte , La 
fonction de Mehler de deuxième espèce est divergente en n=0, soit en z=Cosh(n)=1. Cette valeur 
correspond à l'axe z en coordonnées cartésiennes. Il convient donc pour resptecter la finitude de la 
solution de poser B=0. La solution est donc une superposition de fonctions de la forme : 
U(n,0) = JCosh(n) — Cos(0) (ACosh(r0)+ B Sinh(t0))P , (Cosh(n)) 
So 





Faisons maintenant intervenir le théorème de Mehler-Fock qui permet de construire les 
transformations intégrales du même nom : 


Théorème de Mehler-Fock et transformation de Mehler-Fock 


Toute fonction f(n) peut se représenter sous forme intégrale d'une fonction F(t), par la transformée 
inverse de Mehler-Fock : 
+00 
f(n)= [ar T Tanh(xr)F(r)P i (Cosh(n)) 
Hit 
0 2 
et la fonction « transformée inverse » s'écrit comme représentation intégrale, transformée de 
Mehler-Fock, comme suit : 
+00 
F(x)= | dn Sinh(n) fP ,  (Cosh(n)) 
0 2 
Sous les conditions suivantes des fonctions f(n) : 
- la fonction f(n) doit être continue par morceau, et bornée dans chaque sous-intervalle [In1,1n2] 
de l'intervalle d'intégration [0,+2<] 
- les intégrales suivantes doivent être définies : 


fanfon- et fan fay Log(n) 


quelque soit la valeur a > 1 pris dans l'intervalle d'intégration. 
A ces condition on peut encore écrire : 


f(n) = far T Tanh(zt)P , (Cosh(n)) f dx Sinh(x) fP ; (Cosh(x)) 


+iT iT 


Petite nuance, au point de discontinuité, on écrira : 





LOED) FOW = far T Tanh(zt)P | (Cosh()) [ dx Sinh(x) fCDP , (Cosh(x)) 
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Illustrons cette construction pour le problème aux limites défini auparavant. Dans ce cas la 
superposition peut être vu comme une intégration continue sur le paramètre T : 





T(n,0) = /Cosh(n) - Cos(6) [dr (A(r)Cosh(r0)+ B(r)Sinh(r0))P , (Cosh(n)) 


i 
ou de manière équivalente avec les valeurs limites de 01 et 02 : 
(t)Sinh(t(0 —-0,))+ B(t)Sinh(t(0, -0))) 





P, (Cosh(n)) 





T(n,0) = /Cosh(n) — Cos@) | dr (4 




















Sinh(r(0, — 6,)) vie 
Les conditions limites s'expriment alors sous la forme : 
a - +00 fn) = 
To 
o = i h) E e 
fn) = VCosh(n) Cos(0,) Le AE ICE) > Toha Cose Jar om) 


Or selon les termes du théorème de Mehler-Fock les deux intégrales précédentes représentent les 
des transformations et transformations inverses de Mehler-Fock comme suit : 


z A0) 
f) Cosh(n) — Cos(0,) 


ha) 
ĮJCosh(n) — Cos(0,) 








t Tanh(at)F (t) = B(T) 


t Tanh(zt)F,(t)= A(t) E 


Posons | 





fa) = 











e D z5 farBCE)P, (Coshfn) | f,6n)= [arr Tanh(r)F EP, (Cosh(n) 
1 0 > 0 2 


Cosk ET Ta far Al)P, (Coshfn)) |£) = fdrrTanh(ar)F, EP, (Cosh(n)) 














F(r)= Í dn Sinh(n) J OECS P To 
> 








ne f) 
Res l dn Sinh(n) A CHU) 2 (GORE) 








T nh(n) 
B(r)= cTanh(rr) | dn AUNOM Pi, (Cosh(n)) 
J o  y„Cosh(n)-Cos(0,) -; 








+00 . h 
A(r)= c Tanh(rr) | dn ALN, ) P, __(Cosh(n)) 


La solution formelle s'écrit donc : 





( mt) Bt -0))p 1, (Ces) 





T(n,0) = /Cosh(n) - Cos@) | dr (4 








B(r)=7 Tanh(zr) | dn AS or (Cosh(n)) 


avec 








A(c)=7 Tanh(zr) | dn re Pi. (Cosk(n) 
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Les fonctions A(t) et B(t) peuvent se calculer explicitement lorsque les fonctions limites sont 
constantes. En effet la formule suivante sur les fonctions coniques a été donnée auparavant dans la 
section sur les problèmes en coordonnées sphériques : 
+00 

l = Ter GS)» (cosh(a) 
J2Cosh(a)-2Cos(p) 3 Cosh(xr) Us 
En appliquant le théorème de Mehler-Fock, on peut déduire une autre formule impliquant les 
fonctions de Mehler : 











1 +00 

= [dr C(r)P, (Cosh 

4/Cosh(r) - Cos(6,) l rC(r) a osh(n)) S ce 
= il dr Cosh(r(x -0, JE Cosh(xr) 

2Cosh(n)- 2Cos(0,) à Cosh(xr) tir 


























Comme  C(r)=r Tanh(xr) | dn ÈS Pi. (Cosh(n)) 


Sinh(n) 
P, 
ĮJCosh(n)— Cos(0,) -3+ 











(Cosh(n)= 2 Cosh(r(x -0,)) 


Ti 
>T anae dn Coa 


n Sinh(n) : Cosh(r(x -0,)) 
l í Į[Cosh(n) — Cos (0,) EM Q) 5 (a Sinh(xr) 











Ce qui donne les expressions suivantes pour A(t) et B(t) lorsque les fonctions limites sont 
constantes : 


SG) =T fG)=T 




















B(r)=T. t Tanh(ar) | dn t y A P, (Cosk(n) =T, V2 cree )) 
A2 anne) dE Post) 7,42 CE 0) 


Soit une solution pour le problème intérieur de la différence de sphères : 
(T, Cosh(r(x — 0, ))Sinh(r(0 -0,))+T, Cosh(r(x — 0,))Sinh(r(0, - 0))) 
T(n.0)= J2Cosh(m —2Cos(6) | dr Cosh(rr)Sinh(r(6, — 0,)) 
0 xP, (Cosh(n)) 


——+iT 


2 





x 





On revérifie que la solution respecte bien les conditions aux limites : 
Cosh(r(x — 6, ) 
Cosh(xr) 
Cosh(r(x -0,)) 
P 1 
Cosh(xr) He 








T(n,0,) = T J2Cosh(1) - 2Cos(0,) [dr P, (Cosh(n)}=7 
0 2 





(Cosh(n))=7, 





T(n,0,)= T, J2Cosh(n) - 2Cos(0,) Í dr 
0 
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Lorsque les deux valeurs constantes sont identiques, on doit normalement retrouver une valeur 
constante. En effet en développant le produit des sinus et cosinus hyperbolique: 


Totan cone mhb) Sinh(a + b)- Sinh(a —b) 











2 
Cosh(r(x - 0,))Sinh(r(0 — 6,))+ Cosh(r(x — 0,))Sinh(r(0, - 8)) % 
T(N,0) = T, J2Cosh(1) —2Cos() | dr Cosh(rr )Sinh(r(0, -0,)) 
o [xP „(Cosh(n)) 


a — 0,) + r(0 — 6, ))- Sinh(r(x — 6,)- r(0 — 0,)) | 
so | [+ Sinh(r(x -8,)+7(06, -6))-Sinh(r(x -0,)-7(6, -0)) 
= T, /2Cosh(n) —2Cos(0) Í dr 2Cosh(xr )Sinh(:(0, — 0, ) 

i x P, __(Cosh(n)) 








x 











CO RC | ae EEA e EE De, oo) 
Cosh(r(x - 8) 
Cosh(xr) 
Pour un bol sphérique à fond plat : 

Ici 0=0,. 0 =x 
Si TET, > 








= T, J2Cosh(n) - 2Cos(0) Í dr P, (Cosh(n))=T, 
0 AS NE 


(T, Sinh(t(0 0, ))+ T Cosh(r(x = 0,))Sinh(r(x — e))) 


T(n,0) = J2Cosh(n) - 2Cos(0) | Cosh(rr )Sinh(r(x — 0,)) 
0 XP: (Cosh(n)) 


x 











Cosh(r(x — 0, ))Sinh(r(x = O)) 
Cosh(rr )Sinh(r(x — 0,)) E (Cosh(n)) 


Sinh(t(0 -0,)) 
Coie smaila- - 1, _(Cosh(n)) 


2 





Si T, =0 => T(N,0) =T, /2Cosh(m)- 2Cos(0) Í dr 
0 








Si T, =0=> T(N,0) =T, J2Cosh(m)- 2Cos(0) Í dt 
0 


Si D=T,=T, 
. (Sinh(r(0 —-0,)- at + xr)+ Cosh(r(x — 0, ))Sinh(r(x — 8) 
= T(N,0) = T, J2Cosh(n) -2Cos(0) | dr Cosh(rr )Sinh(r(x -0,)) 
0 x P, (Cosh(n )) 
(Sinh(t(0 - x)+Tr(x -0,))+ Cosh(r(x — 0, ))Sinh(tr(z -8))) 


= T, /2Cosh(n) - 2Cos(0) | de Cosh{xr )Sinh(r(x — 6, )) 
i xP, (Cosh(n)) 


——+iT 
2 


B +2 | Cosh(r(0 — x))Sinh(r(x — 0,)) 
= T, J2Cosh(m) - 2Cos(0) l CS rs ce) F s (Ch) 


Cosh(t(0 - x)) 
Cosh(xr) 





x 








x 




















= T, /2Cosh(n) - 2Cos(0) Í dr P, (Cosh(n))=T, 
0 SE 
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Exemple : solution du problème de Dirichlet sur le domaine intérieur d'un hémisphère plein 


Pour un bol hémisphérique à fond plat, l'angle d'ouverture est de n/2 avec des conditions aux 


limites constantes. 


Kci 05 b =r Si T=T, =T, >T,0)=T, 


Si TET, > 





Cosh z JE Sin +0 = z) wi Cos z Z ]simntete z o) 





T(n,0) = 4,/2Cosh(n) - 2Co5(0) | dr 


x P 1. (Cosh(n)) 


——+iT 
2 





-L : À 4 2(7, -TN2Cosh(m) - 2Cos(6) | dr 
0 





-b 5 + 2(T, - T,)/2Cosh(1) - 2Cos() | dr 
0 











Cosh| t £ 


Sinh(2xr) 


Smi -0 Sai z) i 


- Cosh z 7 Jsme (x-0) 





x 


Sinh(2xr) 


XP. (Cosh(n)) 


——+iT 


Cosh z z) 
2 


2Sinh(r(0 - z)Cos z z) + 


Cos(r(0 - x Sim z z) 











Si T, =0= T (1,0) = AT; J2Cosh(n)-2Cos(0) | dr 
0 


Sinh(2xt) 








Si T=0=7T(7,0)=AT, J2Cosh(n)-2Cos(0) | dr 
0 


En coordonnées toroïdales l'hémisphère (iso-surface 6=11/2 ) a pour rayon le coefficient d'échelle 


du système de coordonnées toroïdales : 


x=alanh(n)Cos(o) y=aTanh(n)Sin(®) z= 


1 
xX +y +z =a’) Tank? + — 
y | (n) Cosh 





XP: (Cosh(n)) 
=  Cosh|T Z \simle( — 0) 
- B 
Sinh(2xr ) -5+it 
To Sins 0 — DIE z) 
2 2 p 
Sinh(2xr) tir 


a 
Cosh(n) 





J= = ravon a 





x 
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Le fond plat est le disque de rayon a (iso-surface 6=11/2 ) : 


NDS E Ut 
Cosh(n) +1 Cosh(n) +1 
… 12 Lim p =0 
ne PEER Le x + y > 70 
(Cosh(n) + 1) Lim p = a 
7— +0 


Regardons le profil de la solution sur l'axe z : 
Si T, AT, > axe z intérieur de l'hémisphère n=0 et 0€ [0,x| Pipu (1) 
SAT 


wa Cosh z JE Sin [0 — z) +T Cosh z 7 Jsime( — o)) 
=> T(n,0)=4V2/1- Cos(0) Í dt 


Sinh(2xr) 


B (z Cost 7 Jsimn(e (z-0))+T, Cos z 2 )sm 0 z =) 


Too) =8si Ji [dr Le 








Il reste à calculer les intégrales suivantes : 


Comme Sinh(a)Cosh(b)= = (Sinh(a +b)+ Sinh(a — b) 2Cosh (4) =1+ Cosh(a) 


> 





fa Sin s[0- 7 \cosnf F y r so (C0) ae 0) 


A Sinh(2xr) Sinh(2rr) À Sinh(2xr 








a Sinhala — 0))Cos? 17) [far Sinhlrlz —0)Coshlen) | F P D 


d = 
l i Sinh(27r) 2N Sinh(2xr) Sinh(2xr) 


0 








. : 2 m 
aia Ro ( a (Fe SO), 1 F q, SC 0) 1 F gy SAE O) 
À Sinh(2xr ) 2e Sinh(2xt) 2e Sinh(2xt) 2. Sinh(2xr) 
De plus 
far Sinh(at)_ x Tan Z2) a far Cosh(at)_ x 1 
> Sinh(br) 2b * Cosh(br) 2b Co z4) 
2b 
far Sinh(r 0) _ Ta K | far Sinh(r (2x -0)) _ x Coran $) 


A aoa A Sinh(2xr) 4 


7 fa Sinh(r(x —0)) 


= Lin 
ptk = 4 4 
Sinh(r(x- 0))Cosh(z z) 
Terme T, = Í dt - 2.21 i Coran $) — Tan +) Fa Ta =) 
Sinh(2xt) 8(2 4) 2 4 4 


E Sin 20 = 7 fees z) i 9 9 
TermeT, = far = (ra | Tan = ) 
à 8 4 4 


Sinh(2xt) 
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Ce qui donne la solution du profil sur l'axe z : 


T(8)= sn(® [r[ 3 coun(?) — 3 Tant à) + Tarf 7 — 2) + nra $) Tan Z — 2))) 
2 2 4) 2 4 4 4 4 
>T(0)= ash (5 ALAE Coran S) -5ta £ + Tarf 7°) Ta + Tan 2) 
2 2 2 2 4) 2 4 4 4 A 
se +2 Es 2 Coran 2)» l Tan $) otan £ )+ 2ranf Z £) 
2 2 2J 2 2 4 4 4 
Or sn coran $) + sran à) 
2 2 4j 2 4 
spoje 7112 + AL is aran 29 )siS) asi £ Jra 2) 
- = ) 0 


5 4 
pps a IÈCE = 2 Jsf 2) + 2cos( 2] = r) 
2 2 2 


Si T#0 T,=0— Tee Ca?) + ran 27 Jsi $) 
T 2 4 2 



























































Posons le facteur d'échelle à 1 (hémisphère de rayon unité) T=1 als cn = Coran $) 
— Cos 











l + Coran $) -1 
Cos 2) — Sin 2) 1+ Cos K — Sin g s$) : 
4 4 2 z 2 
1 


2 
s$) 
2 





Teia 2) 3) 
SE TETE TE 


+ Coran $) +1 
2 














> T(z)= 
V1+ 7° VIÉZ PERF 


Vi+z +741 


fre) Jitz? +z-1 1 | pes!) 
=> Tan à = z+ 
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En comparaison avec une autre expression (voir ci-dessus), on parvient au même résultat par un 
développement similaire : 









































2Sinh(r(0 - r Nos z z) + 















































Coshl t Z 
Ta Cos(t(0 - x Sim + z) 
T0) = E + N21, - T) J1 Cos@) Í dt 
2 A Sinh(2xr) 
„œ Sinh(r(8 — r ))2Cosh? É z) + Cos(r(8 — x Sim + z Jovi z) 
T, + 2 2 2 
= +4(T, -T )Sin 2) Í dt - 
: Sinh(2xt) 
sb + 2(T, -T Sin aji Je 2Sinh(t(0 — x )\(1 + Cosh(er)+ Cos(r(0 - x)}Sinh(rx) 
25 Sinh(2xr) 
LE + (T, -T Sin Il ia Sinh(rx )Cosh(t(0 — z))+ 2(1 + Cosh(rx))Sinh(r(0 — x)) 
25 Sinh(2xr) 
SAT + 2(T, -T Sin dl de Sinh(rx)Cosh(r(0 - x))+ 2Sinh(r(0 - z ))+ 2Cosh(rx)Sinh(r(0 - x)) 
2j Sinh(2xt) 
_I+ +2(T, -T,)Sin aji pA Sinh(t0)+ 2Sinh(r(0 z z ))+ Cosh(rx )Sinh(r(0 — x)) 
2e Sinh(2xr) 
_TL+T + -T in? J fe 2 Sinh(r0)+ 4Sinh(r (0 — a+ 2Sinh(r(0 - x))Cosh(rx) 
A Sinh(2xr) 
Comme Sinh(a)Cosh(b)= = (Sinh(a +b)+ Sinh(a — b)) 
PERS (T, -T,)Sin 2) Tac! a (3Sinh(r0)+ 4Sinh(r(0 - x))+ Sinh(r(0 - 2x ))) 
2 DIR Sinh(2xr) 
PR +(T, -T )Sin 2) Tan +) + Tan 277) — à con ©] 
2 2 (4 4 4 4 4 
Ra + (T, T, )Sin 2) Tan 2) + Tarf ® — z) = £ Coran S) + run(2))] 
2 2 4 4 4 4) 4 4 
hth g Fan FES )+ raf LE sil) l 
2 4 2 4 2j 2 
JA +(T nra 225$) ! ra su )] 
2 4 ZJ 2 4 2 
TER (i — oraz = fs) +1 am (2)} 
2 2 4 2 4 
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Ce résultat est à mettre en parallèle avec celui obtenu sous forme de développement en série en 
utilisant les coordonnées sphériques, soit : 


E r) a (2m)! (4n+3) 
reo- |?) (-1) aY aD P, „„ı(Cos(0)) 





F 








! 
Axe z en coordonnées sphériques,z =r,0 =0=T(z)= D (-1) o CNE aet 
n=0,+0 2 "(n!) 2(n + 1) 


Les profils sont identiques comme on peut le voir dans le graphique suivant : 


— ProfilCartesien(z) 
— ProfilCartesienSerie(z, 20) 








0 & 


Wod u.s 


Du reste on vérifie sur les premiers termes de la série que le développement limité de la fonction : 


1 | ns 


T(2)= 





Z 
V1+ 7° V1+27° +z+1 
coïncide pour les quinze premiers termes au moins. Il est donc fort probable que : 


1 is +51 a COL (4n+3) ou 
T(c) = z = 1 zZ 
A | ur =) 2 ) 2” (nf 2(n+1) 
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Exemple : solution du problème de Dirichlet sur le domaine extérieur de la différence de deux 
sphères 


L'exemple suivant résout le problème de Dirichlet dans l'espace extérieur d'une différence de 
sphères. La solution est également donnée dans l'ouvrage de N.N Lebedev « Special Functions and 
their applications », page 229, section 8.12. 
AT(n,9,p)=0 (n,9,p)eD,, T: 2P) =f(n.p) T(n.9,o),, = A0.) 
ne dr. € [0,+0]x [3,27 +9 ]x [0.27] tq &<2 et (Qe [o.z] etp € [0,7 |) | 

“lou (9 efr,2rler 9, efr,2x|) 
Sans dépendance à p => T(n.9),., = fi(n) T(n,9),_.. = f(n) 


L'astuce consiste à considérer le domaine extérieur du problème comme le domaine intérieur sur 
une section angulaire : [8.8 +27 Į; Dans ces conditions en appliquant le résultat du problème 
intérieur précédent, la solution s'écrit : 

(x )Sinh(r(0 — 9, - 27 ))+ B(t )Sinh(t(0, -0))) 








T.(n,0) = Cosh) - Cos@) | Pa P, (Cosh(n)) 








Sinh(r(@, — 9, -2x)) -5+ir 
2 : ©, (A(c)Sinh(r(2x + 9 -0))+ B(r)Sinh(r(0 -0,))) 
= T,(n,0) = JCosh(n) - Cos(0) l dr E 0) Le (Cosh(n)) 








H Sinh(n) 
B(t)= t Tanh(ar) Í dn A) P, | (Cosh(n)) 
avec 0€ [6,,2x + 9] et 0 /Cosh(m) — Cos(@) -3 








A(c)=7 Tanh(rr) | dn ALU Le (CosA(n)) 


Si les deux fonctions limites sont constantes de valeurs respectives T; et T, il vient : 
T,Cosh(r(x — 6, ))Sinh(r(2x + 9, -0))+ 

+o - TCosh(r(x — 6,))Sinh(r(0 -0,)) ) 

T(n.0)= J2Cosh(n)- 2Cos(0) f dr Sinh(r(2x + 9—0, ))Cosh(art) 


© xP, (Cos(n) 
arr 








x 
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Et si les deux fonctions limites sont constantes et identiques : 


T=T=T Cosha)Sinh(o)= Sa + H- Sinh(a —b) 


Cosh(r(x — 0, ))Sinh(r(2x +9, -0))+ 
+0 o -0,))Sinh(t(0 — 8,)) ) 
T(n,0) = T, 2Cosh(n) - 2Cos(0) f dr Sinh(r(2x + 9, -0,))Cosh(rr) 
"xp, (Cosh) 








x 








a —0,)+r(2x +9 -0))- Sinh(r(z -0,)-1(2x +9 -6)) +) 
+00 Sinh(r(x -0,)+1(0 -0,))- Sinh(r(x -8,)-r(0 - 0, )) 

















= T, J2Cosh() -2Cos(8) f dr 2Sinh(r(2x + 9, —-0,))Cosh(xr) 
i x P 1. (Cosh(n )) 
Sinh(r(3x -0 +9, —0,))- Sinh(r(0 -z - 9, - 0, ))+ 
+o + Sinh(r(x +0 -0, -0,))- Sinh(r(x - 0-0, +0, )) ) 
T(n,0)= T, 42Cosh(n) - 2Cos(0) f dr 2Sinh(r(2x + 9, —0,))Cosh(xr) i 
° x P 1. (CosA(n)) 





Etudions maintenant le cas particulier où l'angle limite 6; est n, alors la forme géométrique est une 
sphère coupée posé à plat sur le disque de rayon a centrée à l'origine. Il vient pour la solution du 


problème extérieur pour un bol sphérique à fond plat : 

Ici 0 =0, 0 =x 

Si TET, > 

(T,Sinh(r(2x + 9, -0))+ TCosh(r(x — 0, ))Sinh(r(0 - z ))) 
Sinh(r(x + 9,))Cosh(xr) 


x 








T (1,0) = 2Cosh(n) - 2Co5(0) Í dr 








Te S = ii Sinh(r(2x + 9, -0)) 
Si T=0=T(n,0)=T,;/2Cosh(n) - 2Cos(6) l dr Ie I Naha D (Cosh(n)) 








: E : — Wi Cosh(r (x- 0, Sinh (t (0 — z)) 
Si T,=0=T(n,0)= T/2Cosh(n) - 2Cos(6) l dr Sane + 8, NORT] P Lo) 


Si ===> 

(Sinh(r(2x + 9, -0))+ Cosh(r(x -0,))Sinh(r(0 - x))) s 
Sinh(r(x + 2 ))Cosh(at) 

x Pid (Cosh(n)) 








T(n,0)=T, J2Cosh(n) - 2Co5(0) | dr 


x 
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Exemple : solution du problème de Dirichlet sur le domaine extérieur d'un hémisphère plein 


Pour un bol hémisphérique à fond plat, l'angle d'ouverture est de 7/2, si les conditions aux limites 


sont constantes sur chacune des faces extérieures, il vient : 


Ici 0-0, = 0,=7 
rsim r{ — o) + 


4 Ticosi(z 7 Jin 22) 








Si T, +T, = T(n,0) = /2Cosh(n)-2Cos(6) Í dr 3 
0 Sin z JCoshtrr) E 














4 ST 
E CE — o) 
Si T=0=T(n,0)=T,;/2Cosh(n) - 2Cos(0) Í dr P, (Cosh(n)) 
: 3x 
o Sin 7 JCoshtrr) 





yo Cosh z 7 Jim — x)) 
Si T,=0= T(7,0)=T; J2Cosh(n)-2Cos(6) | dr - 
Sin 7 JCoshtrr) ra 


E ÉTÉ x o) + Cos z Z Jsimnte(o - =) 


T(n,0) = T, J2Cosh(n)—2Cos(6) | dr = A 
0 Sin Z Jcoshtrr) 2 





w erens 
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Prenons deux constantes limites différentes, la solution s'exprime avec les deux intégrales : 


TA sini s[ 37-0) P, (Cosh(n))+ 


0 Sini( + 77 \Cosnlar) 


m (cosi: Z Jsimn(e(o = x) 








T(n,0)= J2Cosh(n) - 2Cos(0) 














+T, Í dt 3x P, (Cosk(n) 
4 Sin 7 Jovhtre) 2 
_ Sins = o) Lu (cosi: 7 Jsi — x) 
Msez n=0=70)=25n( S) r, f az - +7 | dr = 
0 Sin Z JCosh(rr) û Sin Z JCoshtrr) 


Le long de l'axe z, il faut prendre les valeurs liées à l'angle O suivantes : 


Domaine intérieur — 0 € Ed >z= Coran 2) e [0.1] 


Oelr,2r|- z=27+ Coran 2) e[-,0] 
Domaine extérieur — j 
0e m —z=27+ Cotan $) € [1,+00 | 
2 2 / 


Voici le graphe combinant les profils intérieur et extérieur de l'hémisphère plein le long de l'axe z. 
On note que les discontinuités visibles dans les dérivées ne sont pas réellement physique quoique 
la solution mathématique des deux problèmes intérieur et extérieur de Dirichlet soit a priori juste : 











La solution physiquement vraisemblable n'est donc certainement pas dans la catégorie des 
problèmes de Dirichlet, devant faire intervenir également la continuité des dérivées premières à lal 
normale des frontières du domaine de l'hémisphère. 
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Exemple : solution du problème de Dirichlet électrostatique sur le domaine extérieur de la 
différence de sphères : potentiel dû à une charge placée à l'origine des coordonnées 





Selon le procédé usuel, nous allons retrancher de la solution recherchée le potentiel dû à la charge 
placé à l'origine des coordonnées toroïdales (également origine des coordonnées cartésiennes) 
afin de transformer le problème homogène de Dirichlet sur les deux faces extérieures de la 
différence de sphères au potentiel nul en un problème de Dirichlet avec une valeur limite égale à 
l'opposée du potentiel crée par la charge sur chacune des faces. Cela donne 





T(n,0)= +7 (7,0) > fm- aS = {0 =0,} U {0 =0,} 
x +y +z x +y +z M 
1 B 1 _ Cosh(ņn)—-Cos(0) _ Cosh(ņ)-— Cos(0) 
Jetty tz A | Sinh°(n)+ Sin? (0) Sinh?) + Sin? (0) cf Cosh? (n) -— Cos? (0) 
(Cosh(n)- Cos(0)Ÿ 
_ VCosh(n)- Cos(0) 
VCosh(n) + Cos(@) 
































T, JCosh(n) - Cos(0,) 
c /Cosh(n) + Cos(0,) 


T, JCosh(n)- Cos(6,) 
c /Cosh(n) + Cos(0,) 


T JSinh{r(2r + 9, -0))+ Bi )Sinh(r(0 —- 0, 
Sn Ge) 9 ) el De, (Costa) 








> fi(n)= 





207) = 

















T.(n,0)= Cosh(n) - Cos(0) | dr (4 








i nh(n) 
B(r) =T Tanh(rr) | dn À GDS P, | (Cosh(n)) 
avec 0€[6,,27+9] et o  yCosh(n)- Cos(0,) -> 








Ala)=r Tanh(ar) | an- ARS) p, (Coshn)) 
o  Cosh(n)-— Cos(0,) -3% 
T +00 Pi (x) T +00 Pi (x) 
> Alt )= -r Tanh(at) | dx 22  B(r)-=-7Tanh{xt) | dx 2 
() c l 1] x + Cos(6,) ©) c l 1] q/x + Cos(0,) 
On connait la valeur de cette dernière intégrale à partir d'un développement intégrale connu que 
l'on peut identifier avec les termes de la transformation intégrale de Mehler-Fock. Il vient : 











1 ‘©,  Cosh(t0) Cosh(t0) i 1 
V2,/x + Cos(0) l 3 Cosh(rx) ni Cosh(rx) HRR Go i J2.[x+ Cos(@ tir (x) 
+00 P 1 (x) 


Dai Cosh(t0) Cosh(t0) 
=> Í dx = = — 
1 V2/x+Cos(6) Cosh(rx} Tanh(rr) Tr Sinh(rx) 
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La solution s'écrit alors : 




















P, 

a n -3*7 + Cosh(r6,) 5 T, Cosh(r6,) 
A(r)= T de dx T + CG) T Tanh(rr)- Sinnen) V2 cube) V2 

R Cosh(t0,) 
AE c Cosh(rr) < 

(Cosh(t0, )Sinh(r(2x + 9, —-0))+ Cosh(r6, )Sinh(t(0 — 0, ))) ÿ 
T (1,0) = -2o Ja [Cosh — Cos(6) f dı Cosh(rx )Sinh(r(2x + 9, — 0, )) 
c 0 x P 


i (Cosh(n)) 


Cosh(t0, )Sinh(r(2x + 9, —0))+ Cosh(t0,)Sinh(r(0 -0,))= 
Sinh(r(0 +0, —-0,))+ Sinh(r(0 — 6, - 0, ))- Sinh(r(0 - 2x — 9, —-0,))- Sinh(r(0 +0, - 9, -2x)) 















































2 
Sinh(r(0 +0, -0,))+ Sinh(r(8 — 0, -8,))- 
F +o k Sinh(t(0 — 9, —-0,-2x))- Sinh(r(0 +0, — 9, — >) 
> T,(7,0)=- =D VCosh(n)- Cos(0) f dr Cosh(rx)Sinh(r(2x + 9, —0,)) S 
i x P 1. (Cosh(n)) 
\/Cosh(n) - Cos(0) 
V/Cosh(n) + Cos(0) 
T | (Sinh(r(0 +0, -0,))+ Sinh(r(0 - 0, - 0,))- 
T(n,0) = -= OEO. k Sinh(r(0 - 9, - 0, -2x))- Sinh(r(0 +0, — 9 - ) ; 
J2 | dr Cosh(rx )Sinh(r(2x + 9-0, )) 
x P int (Cosh(n ) 





On peut également montrer que cette solution s'écrit aussi en utilisant les identités sur les 
fonctions hyperboliques : 

Sinh(a + 2i Sinh(a—b) _ Sinh(a)}Cosh(b) Cosh(a + n Cosh(a — b) 
Sinh(r(0 — 6, -0,))- Sinh(r(8 - 8, -0, -2x))= 2Sinh(rx )Cosh(r(x — 0 + 9, +0,)) 
Sinh(r(0 +0, -0,))- Sinh(r(0 +0, - 9, -2x))= 2Sinh(r(x + 9, -0,))Cosh(r(x -0)) 
/Cosh(n) - Cos(@) 

/Cosh(n) + Cos(0) 








= Cosh(a)Cosh(b) 











T(n.0) = 


Sinh(rx )Cosh(r(x -0 + 9, +0, ))+ 
+0 Ë Sinh(r(x +9, —-0,))Cosh(r(x — n 
- V2/Cosh(n)- Cos (0) f dr Cosh(tr )Sinh(r(2x +9,- 0, )) 
"fap, (Cosh) 





x 
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Exemple : solution du problème de Dirichlet électrostatique sur le domaine extérieur de la 
différence de sphères : potentiel extérieur d'un corps au potentiel nul et placé dans un cham 


électrostatique uniforme orienté selon l'axe z négatif 


D'après les lois physiques de l'électrostatique on sait que le potentiel à l'infini reprend la forme du 
potentiel dans un champ électrostatique uniforme, soit zEo. Comme le potentiel est nul à la surface 
du conducteur parfait, il le demeure placé dans le champ électrique puisque la surface constitue 
une courbe d'iso-potentiel. On recherche donc la solution du problème sous la forme d'une 
somme. Le problème qui en découle est alors un problème de Dirichlet sur la surface du corps avec 
convergence en 0 de la solution à l'infini : 

T(7,0)=Ez+T.(n.0) = f(n)=-E4, 0$={0-06}0{6-06;,} 


m = a Sin(0) Sin(0,) E Sin(0,) 
osh(n) — Cos(0) Cosh(n) — Cos(6,) Cosh(n) — Cos(6,) 
(A )Sin(@,)Sinh(r(27 + 9 — 0))+ BC )Sin(8,)Sinh(r(0 -0,))) 
T, (1,6) = -Ea Cosh(m)- Cos(®) | dr Sinh(r(2x +9, —0,)) 














> fin) =—Ea 








ao P, @ . P, (9 
avec 0eļ0,, 2+9] et Aļr)=rtTanh(zrt) | dx — B(r)=rTanh(xr) gema _— 
1 (x-Cos(6,)} 1 (x-Cos(6,)} 
On connait la valeur des dernières intégrales, donner dans l'ouvrage de V.Ditkine et A.Proudnikov 
aux éditions Mir : « Transformations intégrales et calcul opérationnel », où figure une table des 
transformations de Mehler-Fock en page 416, formule n°12.10 : 











a Pit) oda snee T noo C) 5 smi(r-6) 
1 (x Cos(0,)} Sin(0,)  Sinh(xr) n = Cos(0)} Sin(0,)  Sinh(xr) 
. E Sinh(r(x — 0, )) E Sinh(r(x —0,)) 
AE 22: Cosh(xr) Dis ne Cosh(xt) 
T Pi. (Cosh(n)}x 





=> T.(n,0) = -Eça242 /Cosh(n)- Cos(0) [ drd (Sinh(r(x —-0,))Sinh(r(2x + 3 -0))+ 
0 + Sinh(r(x — 6, ))Sinh(r(0 — 6,)) 
Cosh(xr )Sinh(r (2x + 9—0, )) 





Sin(0) 
Cosh(n) - Cos(0) 





T Pi (Cosh(n ))x 
—242,/Cosh(n) - Cos(0) f dt R Pae e E à 
0 + Sinh(r(x - 6, ))Sinh(r(0 — 0, )) 
Cosh(xr )Sinh(r(2x + 9, — 0,)) 


= T(n,0)= Ea 














avec 0 €[0,,2r +3,] 


La linéarisation des produits de sinus hypergéométrique n'apporte pas de simplification notable à 
l'expression obtenue. 
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Simplifions le problème en supposant que le corps est une lentille sphérique symétrique (on 
anticipe un peu sur l'exemple suivant) : 
DT, PETHA 

TP, (Cosh(n))x 


——+iT 
2 


> T(n,0)=-2V2E,ça Cosh(n) - Cos(0) Í dr Sinh(tæ,) Sinh(t(0 - x - a,))- Sinh(r(3x - a, -0) 
° ü Cosh(xr)Sinh(r(2x - 2a,)) 











De plus 
Sinh(r(0 - x - à,))- Sinh(r(3x -æ -0))=-Sinh(r(2x - 0 -(x - &,)))- Sinh(r(2x -0 + x -«,)) 
: E a > Sinh(r(0 - z -a,))- Sinh(r(8x - a, -0))= -(Sinh(4+ B)+ Sinh(A- B)) 


Et Sinh(A+ B)+ Sinh(A — B)= 2Sinh(A)Cosh(b) 

> Sinh(r(0 - x -a,))- Sinh(r(3x -æ -0))=-2Sinh(r(27 —-0))Cosh(r(x - «,)) 

Sinh(ra, )Sinh(r(2x - 8) 

Cosh(xr )Sinh(r(x - a) ie CL n) 
Sinh(ra, )Sinh(r(2x -0)) 3 


2 Sin(0) Cosh(xr )Sinh(r(x -a,)) 
EE Cosh(n) — Cos(0) D on Í ii XP; ip C ) 








> T(n,0)=2V2E,a Cosh(n) - Cos@) | drt 











avec 0 € [x + 0,37 — CA 
La solution donnée dans l'ouvrage de N.N.Lebedev, I.P.Skalkaya, Y.S.Ufliand « Problems of 
Mathematical physics » (problème n° 505 page 241) donne la solution en utilisant des variables 
angulaires modifiées : 
B=r-d=a,=rx-$f;" f=0-2r 

Sin(0) 
Cosh(n) -— Cos(0) 


i z W3 Sinh(r(x - B,))Sinh(rB) 
2V2E,a Cosh(n) Cos(0) l dTT E EA Pi 


2 





0 


=> T(n,0)= 








avec B € [7 +5, -@]= [- Bo: Bo] 
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Exemple : solution du problème de Dirichlet sur le domaine extérieur d'une lentille sphérique, ou 
également d'une union sphérique 


L'astuce consiste à considérer que la solution du problème extérieur sur la différence de sphère est 


applicable avec les deux angles : 


m 
Q =r- à, elo 


z L=n+a, Q, E€ fo, z| < Lentille sphérique 


Q =m- a e| Zz] 4 =7+a, a, el0,r| 
< Union sphérique 
Q =m- a e [0,7] 4, =7+a, a,e[%4r| 


La solution formelle s'écrit alors : 








c)Sinh(r(8x — a, -0))+ B(c)SmA(c (0-7 -a.)), i (Cosh(n)) 


S - TA 
T(n,9)= JCosh(n)- Cos(0) | dr AE ET Da 








n Sinh 
B(r)=7 Tanh(rr) | dn TT PH B (CosA(n)) 
avec Oelr+a,,3x - a] et 9 1/2 








+00 Si h 
A(c)=7 2) dn qan m) Pia (Cosh(n)) 


Si les deux fonctions limites sont constantes, il vient : 
Cosh\r{(x — 0 Cosh(ra 
A =T= B(r)=7, CEE) _ 7 Cou) 














Cosh(xr) ! Cosh(xr) 
: £ Cosh(t(z -0,)) B Cosh(ra.) 
RDF BAC) T2 Cosh(xr) T2 Cosh(xr) 
(T,Cosh(ra, )Sinh(r(3x -a —- 8))+ TCosh(ra, )Sinh(r(0 - x - &,))) z 
T(n.0) = V2/Cosh(n) - Cos(0) f PA Sinh(r(2x - a, -a,))Cosh(xt) 
0 xP, (Cosh(n)) 


——+iT 


2 
avec 0€[0,,2r +0] 0 er +a,,3x - a] 
Si la lentille est symétrique alors la solution se présente comme suit : 
Q=, = et 0e [r +a,,37 -al 
Cosh(ra, (T,Sinh(r Br — à, — 0))+ T,Sinh(t(0 - x - œ, ))) 
T(n,0) = V2 /Cosh(m) — Cos(6) f J: Sinh(r(2x — 2a,))Cosh(xr) 
0 xP, (Cosh(n)) 


——+iT 
2 





x 
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Pour un problème totalement symétrique en condition aux limites, il vient : 
Oelr +5,37 —@, | T, =T, =T, 
Cosh(ra, XSinh(r(3x — à, —-0))+ Sinh(r(0 - x - «,))) 
Sinh(r(2x - 2a,))Cosh(xt) 
x P 1. (Cosh(n)) 
2 


x 








T (1,0) = V2T, JCosh(n) - Cos@O | dr 


. 7 2 Q 


2 — Sinh(r(2x —-0-(x-a, J 
T(n,0)= V2T,;/Cosh(n) — Cos(0) l dr a P 

















: L ©, Cosh(ra,)Cosh(r(2x —- 0))Sinh(r(x - à,)) 

FANER CAC) Le l k 2Sinh(r(x — Qo )Cosh(r(x — &))Cosh(xt) Fe (Costa) 
z = ©, Cosh(ra,)Cosh(r(2x -6)) 

T(n,0)= V2T, JCosh(n) — Cos(0) l dr Coshfr(r — a, )}Cosh(r:) se Li) 


Exemple : solution du problème de Dirichlet sur le domaine intérieur d'une lentille ou d'une 
union sphérique 


On choisit tout d'abord le cas du domaine intérieur d'une lentille sphérique dont les conditions aux 
limites sont sans dépendance azimutale : 
AT(m9)=0 TOI) p = fn) T9), =/(n) 


m 
& € [Za] Q, € [o.z] 
& € Lo z| Q, € [o z| < Lentille sphérique <— Union sphérique 


@, € [o.z] a; er 


(n,9)€ Da Du = fens) e[0+o]x (z -az+a,) tq ace ez) 


La solution intérieure se développe ainsi : 
t )Sinh(t(0 - x +a ))+ B(r )Sinh(r(x +æ, —8))) 
Sinh(t(aæ, +a )) Re 








T(n,0) = /Cosh(n) - Cos@) | dr (A 








$e (m) Sinh 
B(r)= Fran) dn ES LU P t COSA) 
avec 








F fM Sinh(n) 
Alt )=T Tanh d P Cosh 
(E) = Tanlar) f dn EP, (Conta) 
Dans le cas où les fonctions limites sont constantes, on aboutit à une solution pour le problème 
intérieur : 
(T, Cosh(t æ, )Sinh(r(0 — x + a ))+ T, Cosh(r æ )Sinh(r(x + a, -0))) 3 
Cosh(azr )Sinh(t(a, + æ, )) 








T (1,0) = /2Cosh(n) - 2Cos(0) | dr 
0 xP, (Cosh(n)) 


——+iT 
2 


Problème aux limites de Laplace — système de coordonnées toroïdales — p 698 


Lorsque les deux valeurs constantes sont identiques, on doit normalement retrouver une valeur 
constante. En effet en développant le produit des sinus et cosinus hyperbolique: 


7 — E _ i Cosh(a)Sinh(b) - Sinh(a + b)- Sinh(a = b) Sinh(a)Cosh(b) _ Sinh(a SE b)+ Sinh(a = b) 














2 2 
( Cosh(r æ, }Sinh(r(0 — x + a, ))+ Cosh(r a, )Sinh(r(x + a, - 6))) s 
T(n,0) = T, J2Coshm) -2Cos(0) | dr Cosh(xt )Sinh(r(a, + a) 
0 RE. (Cosh(n)) 


Sinh(r œ, +1(0 -r +a ))- Sinh(r(æ, -a)-T(0 -x)) 
+o ; Sinh(r a, +T(x +a,-—0))-Sinh(r(a, -a,)-T(x - A 
= T, V2Cosh() -2Cos(8) f dr 2Cosh(xr)Sinh(r(a, +æ, )) 
"kr, (Con) 











Sinh(r la, + æ, )+ r(0 — z))+ Sinh(r (or + a,)- r0 — x)» 
2Cosh(rr)Sinh(r(a, + &,)) ir 
Sinh(r (or +, ))Cosh(r(0 — z)) 
Cosh(rr)Sinh(r(a, + &,)) P ip (Cant )) 


es (CosA(n)) 








= T, /2Cosh(n) - 2Cos(0) Í dt 
0 








= T, {2Cosh(n) - 2Cos(0) Í dt 
0 





= T, J2Cosh(n) - 2Cos(0) Í dt 
0 


=T, 
On peut noter le cas particulier d'une lentille ou d'une union sphérique symétrique % = %2 7%, 
tout d'abord avec des fonctions limites symétriques soit de valeur identique f@)= f@)= fC): 


Sinh(a)Cosh(b) : Sinh(a + 2l Sinh(a — b) 








A(cXSinh(r(a, + (0 - z )))+ Sinh(t(æ, - (0 - x)))) 





T (1,0) = /Cosh(n) - Cos@) | dt P, (Cosh(n)) 














Sinh(2ra, ) Tone 
= T(7,0)= JCosh(m - Cos(@) | ar ne co 7) p 1, (Cosh(n)) 
=> T(N,0) = Cosh(m) - Cos(@) f dr sekt z) p 1, _(Cosh(n)) 








avec  A(r)=TT. anh(ar) [ dn aT an ru Pa „„(Cosh(n)) 
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Supposons maintenant que les deux fonctions limites soient de valeur opposée: 
f@)=-f@)= f(n). Alors la solution du problème doit exhiber la propriété de parité sur la 
variable 8 7(,._9)--7(7,9), il vient donc : 











Cosh(a)Sinh(b) = Sinh{a + b)— Sinh(a — b) 
2 
T(n,0)= JCosh(n) - Co | AC) CUE -2 adits, --0)) p 1, (Cosh(n) 





Cosh(ra, )Sinh(r (x = 0)) P 
Sinh(ra,)Cosh(ra,) -4+ 

Sinh(r(x — 6) 

a EEG P i 
Sinh(ra,) -5% 





=> T(n,0) = J Cosh(n) - Cos@) | dt A(t) 





=> T(n,0) = /Cosh(n) - Cos@) | dt A(t) 








avec  A(r)=TT. anh(re) [ dn te M AUN ; Pi „„(Cosh(n )) 


On vérifie immédiatement que la valeur s'annule sur la section =n, représentant le disque de 
rayon a. Lorsque la condition aux limites est constante dans ce cas antisymétrique, il vient : 


DS Aer O 














T(n,0)=T, V2 /Cosh(m) - Cos@) | qı Cosh(r a) Sinh(r(x -0)p 


Cosh(xr)  Sinh(ra, 


Cotanh(r 2e puni -0)) P de (Cosh(n)) 








T(n,0) =T, V2 /Cosh(n) - Cos@) | dr 


Reprenons la cas de la lentille sphérique symétrique avec des conditions aux limites identiques, 
pour lesquelles la valeur en n est d'abord constante en puis s'annule au delà de n= n0 (fonction en 
escalier) : 








T, pour nel0.n] 7 Sinh(n) 
= = h 
JG) k pour n€ [n+] B Al) die (er)[ an JCosh(n) + Cos (œ) EL ) 
S 


x, = Cosh(n,) x= Cosh(n) => A(t) = t Tanh(rr)| dx z 


1 x+ Cos(&) 


eheee, (costa) 








T (1,0) = /Cosh(n) - Cos@) | dr Al 
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Exemple : solution du problème de Dirichlet sur le domaine extérieur d'une coupole ou coquille 
sphérique 


Seule demeure une coupole sphérique en considérant que les deux angles du problème aux limites 
n'en sont plus qu'un seul 6=6,. La problème de Dirichlet devient donc le problème extérieure avec 
une seule condition aux limites sur l'iso-surface de l'angle 6-64. La solution devient 
immédiatement : 

4 = 49, = 8, 





A(c)Sinh(r(2x +9, — O)) + B(r )Sinh(r(@ —06, )) P 
Sinh(2x 1) tie 





T(n,0) = /Cosh(n) - Cos(@) | dr 


(Cosh(n)) 








avec 0eļl0,2r+9] et Afr)= B(r)= 5 Tanh(xr) | dn rer P, (Cosh(n)) 
a osh(n)— Cos(0,) -3% 


Dans le cas où la fonction limite est constante, la fonction A(t) se calcule directement : 


f= > A(e)=T, V2 ce o) 





Il vient la solution : 





Cosh(r(x — 6, ) (Sinh(r(2x + 9, — ))+ Sinh(r(0 — 6,)) 
Cosh(rr) Sinh(2x r) Pha (Coshr)) 
Cosh(r(x -0,)) er r(8-0,)) r 


Cosh’ (at) — Sinh(xr)Sinh(r(0 — 6, )) 





T(n,0)=T, J2Cosh(n) - 2Cos(0) | dr 





P, (Cosh(n)) 


Tone 





=T, J2Cosh(n) -2Cos(0) Í dt 
0 








T(n,0)=T, J2Cosh(n) - 2Co5(0) | dr Chlla 0, JCosh(e(r +86) p | (Cosh(n)) 
A Cosh°(xr) tir 


On peut linéariser le produit des cosinus hyperboliques : 
Cosh(r(x -6,))Cosh(r(x +0, -0))= Cosh(e(2r D Cosh(r(28, -6) 
e eE î i Coshc(2r -0))+ Cosh(e(28, -0)) P, (Cosh(n) 
Cosh (x T) a: 
UE 2Cos(0) Í 2 Cosh(r(2x - 6)) P, (Cosh(n + PA Cosh(t(20,-0)) P, (Col) 
Cosh”(xr) St Cosh’ (ar) = 
Pour la suite du calcul on utilise la formule suivante donnée dans l'ouvrage en 
russe « A.P.Prudnikov, Yu.A.Brychkov, O.I.Marichev, Integrals and Series - Volume 3 - More Special 


Functions» en page 181 formule 2.17.24.6 (la formule contient une erreur dans le deuxième terme 
à droite) : 


Tu Cosh(e b) wE 1 5 i - arc 1+ Cos(b ] 














— T(n,0)=T, 
































v Cosh (ar) -+ V2.Jc-Cos c — Cos(b) 
1+ Cos(b 
ArcCotan| | —— ~< 
+ Cosh(r b) 42 1 m 1+ Cos(b) V2 | = 2 
<> [eP | (c)= Arctan = 
2 Cosh’(rt) -4+ z Je- Cos(b)| 2 c— Cos(b) T Je- Cos(b) 


arcod, paa Q | 
f p Cosh(t b)» ce V2 1+c 


E Cosh’ (at) tie T Je- Cos(b) 
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Soit en l'appliquant aux intégrales suivantes : 























b=2x-0 c=Cosh{n) /1+Cos(b) = Vicof ?)- -Vico 2) Arctan(- x)= —Arctan(x) 
V2Cos 2 
‘©, Cosh(r(2r -0)) | à 
l í Cosh°(xr) er St n)= V2 TG — Cos(0 ” Gay JCosh(n)- Cos(0) 
20, 





b=20,-0 c=Cosh(ņn) 4/1+Cos(b) = VC! (2)- icon! 





: Conf 20, -0 ) 
f T Cosh(t(20, see (Cosh(n)) 2 l 1 - Arctan 2 


: Cosh’ (at) St V2 JCosh(n)- Cos(26, - 8) T JCosh(n)- Cos(20,-0) 

















Il vient la solution : 


i icon?) 
1+— Arctan 


V2 /Cosh(n)- Cos(8) T VCosh(n)- Cos(0) 

















ĮJ2Cosh(n) —2Cos(0) 


2 = 
VIco #8 


1 
1- — Arct 
V2 VCosh(n) — Cos(26, -8) T Fe Cosh(n) — Cos(26, -8) 


T(n,0)=T, 























0 
VCosh(n)- Cos(0) + 2 Arci icon?) R 
JCoshn)-CosCe, 6) x "| JCosha)- Coste) 














T 
TODS 





20,—0 
2 Jom- a, icon 2 ) 
T VCosh(n) — Cos(28, - 8) Cosh(n) — Cos(26, -8) 

















Qui est l'expression donnée par N.N.Lebedev dans son ouvrage « Special Functions and their 
applications », page 230, section 8.12, formule 8.12.13. En utilisant, si on le souhaite la fonction 
Arccotangente, il vient : 
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ArcCotan(- x)= - + Arctan(x) et ArcCotan(x)= - — Arctan(x) 
































Vic) 
ArcCotan Ton) co) + 
V 
T,0)=> 
T 26, -0 
na OS V2Cos| 7 ) 
JCosh(n) — Cos(28, - 8) VCosh(n) — Cos(28, - 8) 


En utilisant l'expression avec la fonction ArcCos, il vient : 








PE 
‘© Cosh(r b) V2 1+ Cosh(n) 
Í dt —— —<P; (Cosh(n))= 
1 Cosh(rr) -4+ T7 [Cosh(n)- Cosb) 








b=27r-0 c=Cosh(n 


EOE Vico ? }- 2 co $ | V1+ Coshfn) = Vcosh| 1) 











) 
b) cod?) p 
-i a =. NC a 


Cos (2) 
2 
ArcCos| — —< 


n 
= Cosh| — 
2 osy 1 


x Cosh(n)- Cos(0) 











29, 0) 











Ca 
b=20, -0 c = Cosh(n) J1+ Cos(b) = Vice?) Vico 2 +) 1+ Cos(b = 2 
2 2 1+ Cosh(n) Con! n | 
2 




















C { a- 0 | 
ArcCos rs 
Cos 2) 
Í dr Cosh(t (26, 8); : (Cosh(n))= J2 2 
are (zr) = tit T Cosh(n)- Cos(26, - 6) 
Et la solution s'écrit donc : 
0 20,—0 
Ca?) Cos 248) 
T(n,0) = L ArcCos 2 SE CONE) ArcCos RE 
I 


Cos 2) BTA — Cos(20, - 8) Cos 2) 
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Le cas limite 6,->n, donne la solution du problème aux limites sur le disque de rayon a centré à 
l'origine en coordonnées toroïdales. Il vient : 














T > Vic) ViCon| z — 2) 
9)= | 2+<A A 
TE ren le oE 
Vico $) 
2 








2 
T(n,0) =T 1+— Arct 
PES ge Da Cosh(n)- Cos(0) 


En utilisant la fonction moins usitée ArcCotangente, il vient : 


ArcCotan(x) = - = Arctan(x) => 








> icon? 
T(n,0)=T,] 2-— Arccot 
7 7. T Da VCosh(n)- Cos(0) 


Comme ArcCotan(— x) =T— ArcCotan(x) > ArcCotan(x) =m- ArcCotan(— x) 


icon?) 
VCosh(n)- Cos(0) 








2 
> T(n,0)=T| 2-—| x — Arccotan 
m 


Vic?) 


VCosh(n)- Cos(0) 








2T, 
> T(n,0) = — Arccotan 
m 


En utilisant la fonction ArcCos, il vient : 


2T, Cf) 
T(n,0) = Eu — — >< 


Cosi 2) 
2 


0 elr,3x| 
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Convertissons ce résultat en coordonnées cylindriques à l'aide de ce petit formulaire de passage. 


Formulaire de passage entre les coordonnées toroïdales et les coordonnées cylindriques 3D 
























































7,2 cCosh(n) z _ Sin(6) 
D = = 9e [0,4 0 
R Re Cosh(n)- Cos(8) p Sinh (n) el ; n] nel +] 
2 2 
Tanh(n) = — DE Z Be) 
z°+p° +c V+(p-0) 
2 2 2 2cp 
Cosh(n)= z +p +c Sinha) 
J? +(p-cf 42 +(p+e? Tae ae 
2 
Tan(0)= = a 
2 2 2 
Cos(8)= 2 - =e — 2 Sin(0)= 2 = 2 2 
Jz +(p-c) Jz +(p+c) Jz +(o-c) Jz +(p+c) 
2c° 2e 
Cosh(n)- Cos(8)= = \Cosh(n)- Cos(6) = 
J? F (o - c) Pz + (p+c) J? + (p-cŸ Jz° +(p+c) 




















Doa iie (2 +p?+c°) 11 sin? -| (z + p 2+c°) 
Fe Fr a a a 


Vicos?) Le | i J2? T- HT aR +c) sn +] J onj Jz? + = F = +c) 
gel 


0,7 | +1 0el0,r] 


; (2) -1 0 elx;,2r| f f (2) +1 0 e [7,27] 
©, = Signe| Cos| — | | = o, = Signe| Sin — | | = 
2 -1 6ef[2r,3r| 2 -1 6ef2r,3r] 


+1 0ef3x,4x| -1 0ef3x,4r] 
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Application à la solution du disque 




















0 (z? + p°—c J 
cod £) 2T. s der + c) Pz +(p+c) 
T(n,0) = = ArcCos = = ArcCos (- (e 
4 Cos 2) i jH (z? Ee de à 
; V 


7 +(o-cŸ yz +(o+c) 























2 2 2 2 SN 2 
2m cod VZ+- +de +(0-) Nz loe) 
7 rtp retr- J2+(p+c) 


2 2 
De plus AreCos(y)= AreSin(f1- ) 1— y? =1 cn Dee +7 L c} Vz° . 
PED ACH +( cf 2° +(o+c) 


2 
(V2 +lp- cF +42 lote ) arae ETAGI 
2 2c° 4c? 
Fe 2 D. à 2 2 f2 2 z 
z +p +c +yz +(p—c) Vz +(p+c) (2+06-9 +12+(0+0)) 


En e 2c 


Vz +(p-c} +42 +(p+c) 


2c 
J? +(p-cf + J? +(o+cŸ 
C'est bien l'expression donner dans le livre de J.D.Jackson, Classical Electrodynamics, page 92, 
formule n°3.178. 












































T(n,0) = — sl 
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Pour démontrer la formule donnée en référence : 


‘©, Cosh(r b) : 1 2 1+ Cos(b) || _ V2 1+ Cos(b) 
Í úi an) 0 2 Je — Cos(b) i T arctan c— mo) ze — Cos(b) | c— Cos(b) 














Avec la fonction ArcCotan Arctan(x)= = — ArcCotan(x) 








©, Cosh(r b) V2 1+ Cos(b 
=> Í dr —; ~P, c)= —— 
o Cosh (ar) -54 zc- Cos(b) c— Cos(b) 
On va d'abord commencer par démontrer la formule suivante : 
f dı Coslr a) Pie V2 AA Ie Cosh(a 
Cosh’ (xt) -tit ze — Cosh(a) 1+ Cosh(a) 


Et partons de la représentation intégrale de la fonction conique de Mehler : 


arecon| 














V2 i Cos(t0) 
P = — Cosh d0 ———= 
40) m g Gr) ce + Cosh(e) 
Il vient : 
far beare i ( jee | Cos(r a ) fao Cos(t0) 
osh? 5+ T 3 Cosh(xr) A c+ Cosh(0) 


Cos(ta)Cos(t0) 


2 fao a/c + TC ) + fa Cosh(xr) 


Or on peut utiliser la formule suivante que l'on va démontrer plus loin 


Cosh 7 |Cosh| Z a 0 
L a. Costpr)Cos(rr) 7 J ” z) f 4-Cos(a)Cos(r6) _ cosi 5 Jcosi 2 
0 Cosh(qr) Cosh(xt) Cosh(a) = Cosh(®) 











> 
q cosi 2 | + cosi =) o 
q q 





Il vient alors : 


0 
+0 $s cos à ) 
a, -cos Jfae = cos $] - /1+ Coshlo 
Cosh’ (z HE T Je+ Cosh(0) Cosh(a i Cosh(0) 2 V2 
1+ Cosik } fao 1+ Cosh(0) 2 cos £) _ V1+Cosh(a) 
(Coshla)+ Cosh(8)),fc + Cosh(6) 2 V2 
o a) V1+ cote e J1 + Cosh(0) 

> Î dr ee) 2 (c)= [ee 

: osh (nr) -3+ (Cosh(a)+ Cosh(8))/c + Cosh() 
Posons s = Cosh(0) ds = Sinh(0)d0 = A —1=dOVs" -1 
= f dı n a) Pote) 1+ Cosik l fæ 

5 osh (nr) -3+ 
Posons ç = Cosh(a)+ s 

Cos(r a) RCE 1+ Cosh(a) Tae 1 

Cos) e ada noO ede 1 coder cos) 


























(Cosh(a sNs-1Vc+s 











> far 
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La dernière intégrale peut être assimilée à une intégrale connue faisant intervenir la fonction 
hypergéométrique : 


[dss"(G+b)(G-dŸ"=d"(a+b)" B(A- u-v, u) Fa à) 
d 









































Comme B(x, y) = dre fonction Beta 
fes - mi ya uv TG -u-v (u) | A 
TE Z F| A uv: 

Jdse (ç +b) (s d) d (d+b) T(4-v) 221 sH; vV; d 
En identifiant les paramètres dans la formule précédente, il vient : 
A=1 v=-2 u= u+v=0 d=1+Cosh(a) b=c-Cosh(a) 

1 1.3 Cosh(a)-c 

E 1 . . 
fa 1 RUE í { Tendre) 
ne c4 —-1-Cosh(a).Js +c- Cosh(a) (3) 1+ Cosh(a) 
2 
| 13 care) 
pa 2h L33: 
p fas 1 $ 2’ 21+ Cosh(a) 
1+Cosh(a) Ge —1- Cosh(a).fs +c- Cosh(a) 1+ Cosh(a) 
D'autre part z(a : d z) = Arctanlz)_ p 1, l 2 c - Cosh(a) _ [1+Cosh(a Aretan € Cosh{a) 
2° 2 Z 2°2?° 1+Cosh(a) c— Cosh(a) 1+ Cosh(a) 
c—Cosh(a 
Arctan .|——< 
Wa 1 | 1+ Cosh(a) | 

> dç =2 





ne SV —1- Cosh(a),fs +c- Cosh(a) Ja + Cosh(a)Kc — Cosh(a)) 
Cette intégrale introduite dans le calcul précédent donne : 
c—Cosh(a | 
+00 1 
f dı Cos(t a) p lee J1+ Cosh(a) > + Cosh(a) 


Cosh(rr) -3*7 zV2 “J+ Coshla)c- Coshla)) 
=> Î dt _Cos(r a) p = v2 | one) 


arctan 














A Cosh’ (azt) -+i i z(c- Cosh(a)) 
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AU passage la formule intermédiaire que nous avons introduite précédemment: 
far Cos(pt) = rt) =; jee Cos((p + r}r)+ Cos((p - r}r) 
Cosh(qt) Cosh(qt) 


=; [ar Cosh(i(p + r}r)+ Cosh(i(p — r}r) ss far Cosh(at) ~x 1 
Cosh(qt) Cosh(br) 2b & (z) 
OS 2b 
1 



































far Cosh(i(p +r)_ 1 a far Cosh(ip-r})_ x 

5 Cosh(qr) a Co q: zlet 2) A Cosh(qt) 2q Co de 2) 

Las A ec 

MCE: 

x cosy 70- 2) + Cosy + 2) E Cost Z2 Josi Z) 
ai cos 72+ r) + zie- 2) + cos =+ r) zie- 2) q cos 7P) + Cos 7 





T ar-Colpcostr) z cos 2 Jcosi[ z) 


Coshlqz) q cos 2 ) + cos 7 
q q 


La formule précédente : 
f Ji aal a) Peje V2 c—Cosh(a 
y Cosh (ar) =i ze — Cosh(a) 1+ Cosh(a) 
ayant été démontrée, il suffit de poser a=ib. Il vient alors : 
f da Cosi b) P, (e)= V2 c — Cos(b 
Cosh?’ (zr) -5+ mfe- Cos(b) 1+ Cos(b) 


De plus Arctan € Cos(b) As 1+ Cos(b LT 
1+ Cos(b) c— Cos(b) 2 


ma Cos ilb) p, Ge V2 5 i 2 art 1+ Cos(b ] 


Cosh”(xr) RS 2 mfe- Cos c— Cos(b) 








artan 








arctan 


0 

















‘©, Cosh(rb) B 1 2 1+ Cos(b 
| € Cosh’ re) un) V2 Je — Cos(b) i T arctan c— Cos(b) ] 


0 2 
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Démontrons également la formule suivante : 















































P Si ==) 
far peen - soo. ž q pelia 
à inhlgr a Cos =” helg 
q 
far Cosh(pr )Sin(rt) _ far Cosh(pr) ai irt) -L far Sinh((p + irj )- Sinh((p — ir) 
5 Sinh(qr) is Sinh(qt) A Sinh(qt) 
Comme far Sinh(ar) _ Z Tan Z2) b > |a| 
o Sinh(bt) 2b 2b 
f dè Sinh((p +ir}r) -Z Tan a +ir) 3 J Sinh((p + ire) _ T Tan z(p-ir) 
5 Sinh(qt) 2q Sinh (gr) 2q 2q 
| T Cosh(pr }Sin(rr) BE PE + $ Taf =% = D) 
: Sinh(qt) 2qi 2q 2q 
sin" + =) si" = =) 
Om 2q 2q 
2qi Cod ZE + x) Co Z = x) 
2q 2q 
Sin z(p + ir) Eos z(p T ir) Cos z(p + ir) Sin z(p — ir) 
- 2q 2q 2q 2q 
2qi Co 4 + Deos zi = ) 
2q 2q 
Cos(a)Cos(b)=  (Cos(a +b)+ Cos(a —b)) 
sf) Sim =) 
T q OT 
a Cos ==) + co Z2) E Cosi Z r) + Cof ZP) 
q 
9 
Da Cosh(br)Sin(9r)_1 Sim $) 
>r=9 | far x 9 p 
p=b 0 . (27r) t cos 2) + Cod?) 
Une autre démonstration de cette formule est proposée par N.N.Lebedev : 
©, Cosh(r b) 1 2 1+ Cos(b) 
dt —— © P = 1 Arct = Cosh 
l É Cosh’ (xt) 6) V2 Je — nl T 5 an c— Cos(b) re) 


En partant de la représentation intégrale de la fonction conique de Mehler : 
Sin(r@) 
\Cosh(0)- Cosh(a) 








P, (Cosh(æ))= V2 Cotanh(rr)[d0 
tir T % 
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Il vient en injectant cette représentation dans l'intégrale initiale : 














far Canh b) : (Cone)= 24/2 far Cosh(r b) jja Sin(t0) 
Cosh?’ (at) -5# 2Cosh(rr) m (zt) JCosh(0)- Coshla) 


Inversion de l'ordre d'intégration 


0 








Cosh(r b) Sin(r@) 
a2 dr f do 
Sinh(2xr ) Cosh(o)- Cosh(a) 
a2 ao | p Cosh(t b) Sin(t0) 
Sinh(2xr ) JCosh(o)- Cosh(a) 
22 f ao i ? Cosh(t b)Sin(t0) 
A Cosh(0 EC Je Sinh(2xr) 




















. (+ 
Cosh(r b}Sin(r0) _ 1 Sim $) 


Sinh(2xr) $ cos 2) + co? 





Comme | dt 
0 











n s Coheb) ÈS, a 
EN l a a TRQE 2 l do CTÉI Ca? Jar Cosh(a) 


Posons s = Cosh(0) ds = Sinh(0 40 cos 2) = L> Sim 2) = 


ta 
NI 
= 








Sin $ 
æ mie (CHa)= TE | a | , 
S OS TT 2 T Cosh(a) | Hs + Co] }sm() S — Cosh(a) 
fx l ess d=? 


2V2 4 Š í 
T22 canla) (EO) = s — Cosh(a) 


1+ Cosh(æ 
Js- Cosh(a) = Se - + costa) 


ja 1+ Cosh(a) 
1o 2 


on aa 1 Z O b 
T1} Cosh(a) + Cos z 
Ce) A ç + Cos E) }  |ø=cos(Ż] 
o1 
m 

















fi CA du = te 











C1 
27 | 4 Vs = 
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Si l'on calcule la dernière intégrale, il vient : 


co? 
ArcCos 


: : areco À) EEE 
Calculons Í z £ 


> (u+ BW -e ep EO cos{?] 











1+ Cos(b) 
NV 2 _ 
1+ Cosh(a) 


Co 2)- ORN] du D 2 
2 2 (u+ pW -e e 1+ Cos(b) 


| 2 2 


ArcCos 














arc] 1+ Cos(b | 


n du Pt: 1+ Cosh(a) 
> (u+ BW -e JCoshla)- Cos(b) 


areco 











1+ Cos(b 
J2 1+ Cosh(æ) 
m JCosh(a)-Cos(b) 


Cosh(t b) 
Cosh’ (zt) 








=> far Pi (Cosh(æ))= 
0 DE 








à Arccata | A Cos(b) | 


Cosh(r b) (Ge) Cosh(a)- Cos(b) 








A l l d 
comparer avec la formule l T Core): z TCosh BEA 


Or arc] ren) = arc | l Cos(b) | 








1+ Cosh(a Cosh(a)- Cos(b) 








Puisque ArcCos(t) = AreCotn| 


t t 
< ArcCos = ArcCotan(t) c.q.f.d 
=) =) 
1 t 
Sin(ArcTan(t)) = 
e de 


De plus Cos(ArcTan(t)) = 








ArcCotan(t) = = — ArcTan(t) 
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Exemple : solution du problème de Dirichlet électrostatique sur le domaine extérieur d'une 
coupole sphérique : potentiel dû à une charge placé à l'origine des coordonnées 


Selon le procédé usuel, nous allons retrancher de la solution recherchée le potentiel dû à la charge 
placé à l'origine des coordonnées toroïdales (également origine des coordonnées cartésiennes) 
afin de transformer le problème homogène de Dirichlet sur la surface de la coupole : potentiel nul 
en un problème de Dirichlet avec une valeur limite égale à l'opposée du potentiel crée par la 
charge. Cela donne 





GS = {0-06} 














T, T, 
a NON 
x +y +z x +y +Z | 


1 : 1 __ Cosh(ņ)—Cos(0) _—— Cosh(ņ)— Cos(0) 


EENE iy | Sinh? (n) + Sin? (0) Sinh? (n) + Sin? (0) cf Cosh? (m) - Cos? (0) 
(Cosh(n) - Cos(0)Ÿ 


_ VCosh(n)- Cos(6) 























__ T, JCosh(n)- Cos(6,) 




















/Cosh(n) + Cos(0) >S) c Cosh(n)+ Cos(0,) 
760,9) = JCosh) = cos) | dr Ale) AECT +9, -0)+ Sinh -0 D (cosh(n) 
A Sinh(2x r) She 








avec 0€[0,,27+9,] et Afr)=7 Tanh(ar) Í dn T sn Pi. (Cosh(n)) 
o osh(ņ)— Cos(0,) -3+ 

















sis n Sinh(n) Cosh(n) - Cos(6,) 
FAR cel PT T CRO John + CO) net ) 
+00 Si h 
> A(r)- -for Tanherr) dn T A Pa COG < x = Cosh(n) 


T +00 Ft (x) 
= A(r)=-"7 Tanh(rr) | dx ——©>— 

c 1 x+ Cos(0,) 
On connait la valeur de cette dernière intégrale à partir d'un développement intégrale connu que 
l'on peut identifier avec les termes de la transformation intégrale de Mehler-Fock. Il vient : 











1 + Cosh(r6,) Cosh(r0,) i 1 
= | dr 0 P Svo -t Tanh d P 
V2 ,/x+ Cos(0,) l i Cosh(rx) ae Cosh(rx) one Gr) “R |x + Cos(0, a 
P 
ie 10 Cosh(re,) B Cosh(t0,) 








Š l dx V2 /x + Cos(6,) | Cosh(rx } Tanh(rx) Tr Sinh(rx) 
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P, (x) 


——+iT 
2 


Jx + Cos(6,) 








Cosh(t0,) Jaa T, Cosh(r8) 5 
t Sinh(rx) c Cosh(rx) 
Cosh(r6,) (Sinh(r(2x + 9, -0))+ Sinh(r(0 — 6, ))) P, (Cosh(n) 
Cosh(rr) Sinh(2x r) Hi 
(Sinh(r(2x + 9, — 0))+ Sinh(r(0 - 0, )  (Sinh(r(2x + 9, —- 8))+ Sinh(r(8 —-6,))) 
Sinh(2x r) g 2Cosh(xt)Sinh(x r) 

: Cosh(rx)Cosh(r(0 — 0,))- Sinh(xr)Sinh(r(0 -0,)) u Cosh(r(x +0, -0)) 
n Cosh(xt) z Cosh(xt) 
Cosh(r6, )Cosh(r(x +0, -0)) p 

Cosh’ (tr) rit 
La solution du problème aux limites est Lie la suivante : 


T(1,0) = fCosh{p - ca l J2 far Cosh(r8, Cosh(r(r +0, -0)) p , (ct) 
C 


A(r)= -Er Tanhlar) f dx 


pgi t Tanh(at) 
c c 








T,01,0)= -V2 Cosh- Cos(0) Í dr 














=> T. (ņ,0) = -42 V2 Cosh(n = Cos(@) il dr 

















JCoshm) + Cos(@) ‘À Cosh’ (rx) a 
L'intégrale peut être calculée jusqu'au bout en linéarisant le produit de cosinus hyperbolique : 


Cosh(r(x + 20, -0))+ Cosh(r(x -8)) 
Cosh(r6, )Cosh(r(x + 0, -0))= 5 























f dr Cosh(t0, JCosh{r(x +0, — 0)» -Aona 
5 Cosh’ (t7) SE 
- f p3 Cosh(r(x +20, — 9) + Cosh(r(x — 0)» 1, (Cosh ) 
25 Cosh°’(rx) 

Or Eaa (e) -32 l Arctan o 

o Cosh (ar) -54 z Je- Cos(b c— eao 

1+ Cos(b 1+ Cos(b 
ArcCotan| EE ArcCos| LT 

fa Cosh(r b) „ (c)= | = | q l+c | 


à i Cosh?’ (ar) -5+ T Je- Cos(b) x Je- Cos(b) 
areco | — Cos(26, - a) 


f g Cosh(r(x +20, 0) , An 1+ Cosh(n) 
Lee Cosh? (er) Pa CA z -JCoshfn)- Cos(b) 


arc] 1= Cos(0) (0) | 


‘©, Cosh(r(x -0)) S \1+ Cosh(n) 
e a Oo 


— T(7,0)= Lo Cosh(n) a ArcCos = Cos(26, 6) + ArcCos RQ 
7 c(VCosh(n)+Cos(@) ~ \ 1+Cosh(n) 14 CosA(n) 


Expression finale qui s'avère relativement condensée. 
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Représentations intégrales des fonctions de Legendre et des fonctions associées de Legendre 


Références : 
- A.Erdelyi.H.Bateman-HIGHER_TRANSCENDENTAL_FUNCTIONS_VOL_ I Toroidal functions section 


3.13 page 173 : 

- l. S. Gradshteyn and I. M. Riyzhik, Tables of Integrals, Series and Products, Fourth (English) Edition 
section 8.85 Toroïdal Functions (page 1053) 

Valable pour u entier = 0,1,2 

w+u+1) Sinh" © Sin™” (0) 

G-u+1) Vera ; : ) 3 (Cosh) + Sinh) Cos(0)) 





PA (Cosh) = 














E T(v +1) í Cos(u0) 
2m I(v-u+1)%  (Cosh(n)+Sinh(m)Cos(0))" 
1 
r{nima,) : à 7 . 2m 
P", (Cosh(m))= 2e SIA — ja ne) m 
a r|» -m+s) 2" Jarm- >) 0 (Cosh(n)+ Sinh(n)Cos(0)) "> 
Ge 
D an+) ,, 
_ CD” 2 do Cos(m0) 


27 r|» Sm ;) o (Cosh(n)+ Sinh(n)Cos(0))*z 


1 
rfnt,) i Cosh(mr) 
dt OS 


r|» Cm >) 0 (Cosh(n) + Sinh(mCosh(n} "à 


r|» +m+ 5) Los( Coram (2)) 


=(-1)” i | dt Cosh(mt)(Cosh(n)— Sinh) Cosh] 2 n>m 
r|» + | g 
2 





o” (Cosh) )= CD” 


1 
n—— 
2 
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- Lebedev "Special Functions and their applications" 


- 1. S. Gradshteyn and I. M. Riyzhik, "Tables of Integrals, Series and Products", Fourth (English) 


Edition section, 8.712, 8.715 (page 1002) 


1 
JISinh" a)’ a cosi (v + Y) 








P" (Cosh(n))= - n>0 Re(u)<} 
Jaf r n) 0  (Cosh(n)-Cosh(t)) 7 2 
NI a) 
or (Cosh) = jZ CEN fi r n>0 Re()<> 


r{, B n) a (Cosh(t)- Cosh): 





Re(u+v)> -l1 








1 
” cosi (v + }) 
B(Cosh@m}= *2 Jr 2 
T ò (Cosh(n)-Cosh(t)} 
7 1 17 1 
ECESO 190: g (Cosh(n)+ Sinh(m)Cos(6)) "2x l j (Cosh(n) + Sinh(n)Cos(0)}" 
P,(Cosh(n))= l lfae( (Cosh(n) + Sinh(n)Cos(0)) = = fao (Cosh(n) + Sinh(n)Cos(0)) 
0 m 0 
Q,(Cosh(n)) )= fao - 


3 (Cosh(n)+ Sinh(n)Cosh(0))" 


- l. S. Gradshteyn and I. M. Riyzhik, Tables of Integrals, Series and Products, Fourth (English) Edition section 8.7145 


(page 1002) et 8.742 (page 1007) 
- Lebedev section 7.1. 





rat) — u+1)° 


pe) MEET ar fe 1Cos(t ) cos(ut) Re(z)>0 u =0,,2,.. 
zT(v +1) 


UE Conf» + +) 
u= of; }- Jr => P (Cos(0))= 2 fat 2 - 
2 T 3 (Cos(t)- Cos(0)} 








Co» + k z s)] 
Cos(vr }P,(Cos(0))- Ž Sin(vr)O, Cos(0))= ja- (Cos(0) : ( J} 
9 os — Cos(t))2 


v=n- 5 > Cos(vr)=0  Sin(vr)=-(-1)" 








P_,(Cos(@))= s fa Costa l 
a: mi (Cos(t) = Cos(0)}: 
ren Cos(n(t-x) _ 1 Cos(nt) 
Q (Cos(0))= Í dt = [dr 


2% (Cos(8)- Cos(n} 2; (Cos(8)- Cos) 


P“(:)= Tesut 1)? far ( + V7 _1Cos(0) Sin™”(t) Re(z)>0 y =0,12... 
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- Lebedev "Special Functions and their applications”, section 7.10 formule 7.10.10 

- |. S. Gradshteyn and I. M. Riyzhik, "Tables of Integrals, Series and Products", Fourth (English) 
Edition section 8.712, 8.715 

- Louis Robin 1959, Fonctions-Sphériques-de-Legendre-et-Fonctions-Sphéroidales Tome III, page 
113 


Of(z)= a r( ce +1) (z pae Re(v+u)>-1 Re(u)>-1 arg(z + 1)] <T 





2"*l ro +1) ! (z-e 
1 
= == = 
u=m v > 
z = Cosh(n) 








g Qi rismi) gu a-e 
Q” (Cosh(n)) = | i (Sinh(n)) fat 








2": (1) à (Cosh(m-1)"+" 
t = Cos(0) 
1 
jé riem) > | 3 
Q” (Cosh) = — ~ (Sinh) |do MEL — 
E 2": rfi) à (Cosh(m)-Cos(0))""=" 


imm 


y : | 1\f Cos(10) 
Q” (Cosh) = <= —(Sinh(n)) r(m+ ) de | 
L3 v27 2 l (Cosh(n)- Cos(0))"* 


Comme e™7 =(-1)" pour m entier 














ai CD" po, P , à Cos(10) 
Q” (Cosh(n)) = (Sinh) "T| m+- || 40 
E v27 É (Cosh(n)- Cos(8))"" 
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Autres représentations intégrales des fonctions de Legendre de degré demi-entier. 


Référence : 


- Lebedev “Special Functions and their applications”, section 7.10 formule 7.10.10 
- |. S. Gradshteyn and I. M. Riyzhik, "Tables of Integrals, Series and Products”, Fourth (English) 


Edition section 8.712, 8.715 


- Louis Robin 1959, Fonctions-Sphériques-de-Legendre-et-Fonctions-Sphéroidales Tome III, page 


113 


- S.C.Loh, 1959 « The Calculation of Electric Potential and the Capacity of a Tore by means of 


Toroidal functions » formule (6) page 699 





ot (2) =° ee Del [ae ee =] Re(~v +u)>-1 Re(u)>-1 larg(z + D| < 7 


(z-e 


Dv+ T(v+ 
Partons de la formule AONE en HU : 
u=m v gjan 
2 
z = Cosh(n) 


Puis réalisons une intégration par partie : 








ñ e™ (rem) , oi, (-r}> 
Q”, (Cosh (= — i (Sinh(n)) Í dt 1 
2 2 rf + ;) a (Cosh(n)-t}""2 


r{im+ z 

imm A m : 21 

O" (Cash) = E > Sinon)" | ao — OI 
R 2 2 rfi) 








seal 
0 (Cosh(n)- Cos(@)) 2 
La première intégration par partie donne : 


i (DE 1 
far i >u'= i et v=(i-#) 2 
a (Cosh(n)-1)"2" (Cosh) - t)z" 

-1 


(m DoS 1 )(Coshon =)" RS 





> U = 





+ 


(1 -t E 





e z E ja (a-e) E 


a (Cosh (madari)? (Cosh (meii )i Cosin- 





fn] a s(a] 


: 1 
a (Cosh(n)-1)"2" CEA (Cosh) -t)z 


P 


1 
1 
De même HG J=] =0 pour p=1,l-1 
-1 


=l 
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En itérant (l-1) fois de plus cette intégration par partie, et avec le résultat de la ligne précédente, il vient : 


1 a m 
jae E le 








f 1 
i (Cosh! [messe et-2) mezer)" Cosroe 
l d' ni 
1-2 p 1-2)? 
1 r_ dt A era i | 





= Í 1 I 
(m+ L+- +41-2).{m +3) “q-e E (Cosh(n) -t)z 
2 2 2 
Formule de Rodrigues pour les pôlynomes de Tchebychev 
1 l 1 
C1) F - (c 2D = (21-1)(21-3)...x1x T(t) 


Si t = Cos(9) > T,(Cos(9)) = Cos(19) 





oz a-e) , (21-1X(21-3)...x1 L dt T(t) 
(Coship? 


q (Cosh(m-1)" 3" i CE rene 
2 2 2 





x 








1 
ima [| n+tm+— i 
Q” (Cosh(m)) = — < (Sima)! — PE D da T() l 
É "3 rfn+s) CROSS OR (Cosh(n)-t}""2 





: -U p 

r r(:+2) mes fmes+t-2) mes) (6) (Cosh(m-1)""2 
Conmer(1+3)=[1 JL firl) nf )=v 
rfitmes)-(meger-ifmeser-2) (memes) 
































Il vient 
O", (Cosh) = — (simon T{m +3) Lewu 
z 2 Jr a(i} (Cosh): 
dt 
P t=Cos(9) dt =-Sin(9)d9 =-V1-t d9 = - = 49 
osons os(4) in(9) > NS 
D Cho ee (Simont m+ ] d9 T, (Cos(9)) -= T,(Cos(9))= Cos(19) 
z V27 27% (Cosh(n)-Cos(9))"" 
E DE n E ie Cos(19) 
=> Q" , (Cosh(n)) = £= (Sinh) rhm | d9 - cqfd 
z 27 s l (Cosh(n)- Cos(9))""2 


Comme e7 =(-1)" pour m entier 








z CD a ? | | Cos(19) 
Q” (Cosh(n)) = (Sinh(n))'T| m+— || d9 - 
L3 v27 2 l (Cosh(n)-Cos(9)}"*2 
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A partir de la représentation avec les polynômes de Tchebychev, on obtient également : 
Pour m entier 


a (Cosh(n)) = C 1)” (Sino rm F 1 J | dt ; : T(t) 
2 n 











2 1 el 
22 V7 “a-e E (Cosh(n)-t}"*2 
Posons t = Sin(9) dt = Cos(9)d9 = V1- ť d9 => = =d9 
l1-ź 


T,(z)= Cos(ICos™(z2)) T(z)= Cof © TL cost sim +s Ti )sint Sin (z)) 


= T(Cos(9))= Cos(19) T,,(Sin(9))=(-1) Cos(219) T,,.(Sin(9)) = (-1) Sin((21+1)9) 








Qi (Cosh(n)) = E 1)” (Sinh(n))" r|m E 2) i d9 T,(Sin(9)) 








V27 Cosh(n) — Sin(9))"*2 
7 S TENES | i Cos(219) 
Cosh(n)) = Sinh(n))"T| m+— |f d3 
=> O7 (Cosh(n) = Te (Sinh )'T| m + > | Ciona 





D” E E T Sin((21 + 1)9 
= 2,1 (Coshtn) = G c" (Sinh) r|m +) (SAGE DE) 

É o (Cak Sin(9))"* 
La formule est également valable avec les indices m entiers négatifs : 


cy G”): 
o” (CAGE Lx Gien’ | Í d9 Cos(19)Cosh(n) - Cos(9))"= 





dus 








1 
I l+—-m 
dm m T(v -m+ 1) a -m E m | ) m 
PrE EU ne) IR 1) i 1 2) 








T(v+m+1) 7 


l r{itsim) nf) 
>Q” (2)= 


faa Cos(13)(Cosh(n) - Cos(9) 
H 2x r{i+ Eo PUARE 


r( B S r(1+-m) 
OE TOEO 1)" | Es 
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A partir de la relation pour m entier : 
: =1) a . 1 Cos(19 
O” , (Cosh(m))= l coi (Sinh(n)) (m+ 3 fas ee 
z 27 0 (Cosk(m)-Cos(9))""7 


et en partant de m=0, puis par dérivation sur le paramètre n, on arrive aux résultats suivants : 














o Con- [49 — [LC 91, (Cosh(n}=-sinnem ag — 09) 
z 0 (2Cosh(m)-2Cos(9)} "2 0  (2Cosh(n)-2Cos(9)} 


dO” ;(Cosh(m)) 


> ShO, Cash) = -Simen f ao — CCE 


2 o (2Cosh(n)- 2Cos(9)}2 
=> Sinh) QÙ, (Cosh(m)) = Q" , (Cosh) 














E E E EE E a 1 at Cos(19) 
z Sinn) dn| Sinn) dn| Sin dna  (2Cosh(n) -2Cos(9)} 
+ o (Cosh) = (= 1)" (2m - uf d9 Cos(19) 


0 (2Cosh(n)— 2Cos(9))""7 i 
Ne 127 r(m+ ajas ja Cos(19) 
GU 








Vr (2Cosh(n)- 2Cos(9))"" 


os Be CAT (m+ fas EOSL) s 
2} (Cosh(n) - Cos(9))"*z 





Ym Q” DE 


se y2 o + -Jj d9 Casa) 
va 0 (2Cosh(n)-—2Cos(9))"*2 


=> Q” ,(Cosh(n))= Sinh"(n)o!")(Cosh(n)) 








On a également pour la Fo. première de la fonction toroïdale de première espèce : 


P' ,(Cosh(n))= Sinh(n )P° (EG ) 
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Autres intégrales dérivées de la formule pour m entier : 






































m Elf a | LÍ Cos(19) 
O” ,(Cosh(n)) = (Sinh(n))"T| m+— || d9 - 
Ri 27 Š l (Cosh(n)- Cos(9))""2 
Il vient : 
P C1)" ro, à | | Cos(19) 
Q", (Cosh(n)) = (Sinh(n))'T| m+— || d9 
Le 27 z | (Cosh(n) - Cos(9))"”"2 
rfm+ J- VaT (2m) E Vz (2m-1) R 
2) 2”T(m) 2” (m—1) 
o° (Cohn = fao — 509 
E aa © (2Cosh(n)-2Cos(9)}> 
OA (Cosh) = AD (simnemy 229 zÍ d9 M 
(n-1} 3 (2Cosh(n)—2Cos(9))"*z 
A Cos (L9) -= 07 Coha) fas Cos(l9) 
o (2Cosh(n)-2Cos(9)} “> o0 (2Cosh(n)+2Cos(9)}: 
Fagg — Sin(@)Sin(19) Í qg — Sin(@)Sin(18) 


d9 Sin(0)Sin(19) 


7 =lQ  (Cosh(n)) 
0 (2Cosh(n)-2Cos(0)} l-7 


Cos((1-1)}9) 


o (2Cosh(n)+ 2Cos(0)¥}: 


d9 


1) 0° , (Cosh) 


=10 _,(Cosh(n)) 





T 


_ gg Sin(O)Sin(9) 


0 (2Cosh(n)- 2Cos(0)) 


= d9 
2° NS —2Cos 





d9 Cos(0)Cos(19) 


o0 (2Cosh(n)- 2Cos(8)P 


| ae 


Cos((1 +1)9) 


OE 


lor, ı (Cosh(n))- Q' C0) 


Cos((1 +1)9) 


(2Cosh(n) - 2Cos(0)P 





ED ! [as 3 
0 (2Cosh(n)-2Cos(0)} 2 | an 0 














d9 


j Cos((1-1)9) 


(Cosh(n) - Cos(0)} 

















1 
0 (2Cosh(n)+2Cos(9)}""2 





=Q", (Cosh(n)) 
14 


1)" Sinh” (7) 2m -1) 


= jag Oct Lo. ,(Cosh@m)- 0! Cat} 
0  (2Cosh(7)-2Cos(6) ni 

A TT AL 
0 (2CosA(m)- 2Cos(9)): SARGI) o0 (2Cosh(n)+ 2Cos(9)}: Sinh(n) 
tra Cos(19) i Qı (Cosh(n)) E Cos(19) : gi (Cosh(n)) 
0  (2Cosh(7)-2Cos(9)} 3Sinh” (n) 0 (2Cosh(7)+ 2Cos(9)} 3Sinh°(n) 
f Cos(19) z” (m—1} 

d9 = Q”, (Cash) re 
o  (2Cosh(7)—2Cos(9))""7 2(-1)" Sinh” (q) (2m —1) 
f Cos(19) 2: m —1)! 

d9 (m-1) 
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Quelques développement en série utilisant les fonctions toroïdales, dans le cadre de l'inverse des 
distances en coordonnées cylindriques 3D 


Ce que nous avons nommé la « partition de l'unité 


1= 20 VCosh(m)- Cos(0) Ha Q, (Cosh(n))+ > Q à EE) 
T 


1 _2V2/f1 
oo) Cos(0) x 














10, (Cosh(m))+ D 0 : EE) 


n=1,+00 


Le développement de la solution élémentaire (une fonction de Green) pour les problèmes 
d'électrostatique : 
En coordonnées cylindriques axiales r(p,9,z) r(p',9',z') 


—r|= Je + p°+(2- 27) -2pp'Cos(9— 8) 
= p’ + p°#(z-2) 














=tr-r|=/2pp' JB-Cos(9-39) B>1 























2pp' 
1 1 
LE {0.68 (#2 X 0, (Ben Ka) 
Démonstration 
eu l Posons B = Cosh(n) 
r-r| /2pp' JB-Cos(9-3") 
1 242] l 
Tosha- Cae 20 Gon 2 9, : «(Com costo =) 
1 2V2 (1 
<> VEZORE (Ge F PA "i 








l 1 1 1 
> = = B)+2 Q : B)C (n9) 
r-r] /2pp' JB-Cos(O) xpp' foe E 220 ra | 
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Partons maintenant de la formule plus générale : 


(1) 1)" rm +j Cos(n 9) 
O Conan r fas 
g 2x 24% ( 


Cosh(n) - Cos(9))"*3 
1 


(Cosh(n)-Cos W 


= 1 » EL le, (Cosh(n))+2 Ÿ_ Q" | Cash) 
e" 


(Cosh(m)- Cos(9))" +  (Sinh(m)"T ne 








Développement en série de Fourier de 





B = Cosh(n) = Sinh(n) = JB? =1 = (Sinh(m)" = (82 1} 


2i z = | Eles, (B)+2 2 2. (Eco) 
B-Cos(9) B? -1 {m+;) 








k -r|= 290 JB -Cos(9-3) > r-r" = (20p)"2(8- Cos(a- 2)": 














1 1 1 
rt pp)" (8 _ Cos(a - 9')"*z 
= = {e: (8)+2 Z, 9° 1(B)Cos(r (8-3) 
mr! 2" (8° -r|m LEl EN T 
> -=— (e, (B)+2 X, O°:(B)Cos(n w-s)] 
or] xp | 7 ROLE 
1 1 {o 
+ —— = 9',(8)+2 $, O' ,(B)Cos(n -s)] 
r-r z |B” -1(pp ÿ n=l+o 772 





= = T z, (8)+2Y © RCE -s)] 


3 n=1,+00 
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Formules avec les puissances positives de la distance : 


a =, m = fo; (B)+2 £ Q” "iP Gasla (8-8) 
r-r r(e -Er m+ Nr(opY m 2 n=1,+00 


1 = (- 1)” { =m ' 
E = Q7(8)+2 Y 9", 1 (B)Cos(r (9-9 ) 
r-r 2” (B? SE rm No je 2 n=1,+00 








T msm] 











"O To=" 
Ij n+ 2 + m) 2 2 








2 y" (8° AP rfi) r(n+3-m] 
=>> |r Z r77 2 a o", (B) +2 a = 20 (B)Cos(n(9 = g') 
(> T mz 72 f Ș m) n=1,+00 f n= 








2" ( y 2} (3-7) {ns mm) 
mi _ 2"(op)"(8 o" (+2 D + ——<0" (B)Cos(n(9— 8) 


> |r-r 
Vr pp'T on 





2m-1 2”( y=(g 1) r(n+3-m) (34m) 
Re PS > Q",(B)Cos(n(9— 9") 























nn 1 I I E 
rtmNr erasi +m)r(-) 2 
2 2 2 
m=1>|r rt= SVP Ho, (8) 22 Gr j2 (9-9) 
2F y| (pb) 
3 _ 48(pp'}2(B"-1 -= 1 2 ' 
EET E cau rf 2a arare Pee) 
5 3 ( °, (8) 
_ p _ GApp'E(8-1p| “> 15 p 
ARE x 15 Lt 25{4n° 94n - j2 K. 
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En introduit la notation du bi-factoriel, bien la manipulation du bi-factoriel soit plus délicate en 
général que celle de la fonction Gamma, nous obtenons 


na __2"(ppY"(82 - 1} 


F-r] 





fin) {ns -m) 
9 PB) +2 a —< 9" 1(B)Cos(n(8 — 3) 
Vz pp: (3-7) 2 fit) er(ns+m) E 


(34m) Ti Cn-27) -2 (mr {> m0) En 
(m=1}t= [T(@m+1-21)=[](21-1) 


l=1 





hia 
ll 
= 


T nt y) 5 1 
= 2?" 


Lo) e ne cer) 








= 

+ 

3 N 
Ne 7 


l=1 





= 


= j [G+ mt) = 22" ] [(1+2(n-0) 
2 




















rfa + 5 + m) 1=2m I=m 1=2 
>- 4-2" ] J(n+2m+1-21)2 T (On + 2m +121) [ [(on + 2m +1- 21) 
Tnt =) {=l I=1 l=m+1 
= (2n+2m-1)! Pr = (2n-1)}! 
2n+2m+1-21)= 2n+2m+1-2/)=T [(2n+1-2/)- 
[I n+2m+ 1) n-i)! IK n+2m+ 1) [I n+ 1) Gn- 2m1) 
l=n 
s , [I (22 +1-21/) ni} 
Puisque (2n-2m-1)!= [| [(27-2m+1-21)= [Tr +1 21)= -= VE $ 
11 l=m+1 [Oeris 21) (2n+1-2/) 
{=l 1=1 
1 De) RD 
Caen) Gen -m re m) (Cm _ CD"2" 
rfn+ a ~m] (2n-2m-1)! (4 +m) 2"(2m—1)!Vr  ((2m-1)) 
Q".(B) 
= pr pm (EPY (e -| 7 XD" (my Y Lr=2m 21) Cos(n(9 - 9' 
Pol A EP A ran A eee) 
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Vérifions qu'il n'y a pas d'erreur : 





m O0": (8) 
O mi om (PPY (B-10 (2m-1)) “= (2n-2m-1)! -n = 
a x pp ((2m—1)) a D Grimm US) 








m=1=lr-r|=4 a ep" Ho, (8) >. =) o aCA as) 














m=2=>|r -r| =48 (ep (g° -1) 2 Z H25 : Q° ,(B)Cos(n(9 - 9) 
a g T 9 (4n -14m 9) 
5 3 (9°, (8) 
2 LR (op'} (8° -1} 72 15 3 q 
m=3=> r -r| =64 E o 22 ar- ie e aR 9) 
; Q",(8) 
m=4 -rf 22561008" 1) 105 
T PP' 





105 ; pe 
r2 (4n? -14n -9 fan 25 jar 49) P )Cos(a(3 - 3')) 
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Quelques propriétés des fonctions coniques de Mehler 


Fonctions coniques de Mehler de première et deuxième espèce, avec un argument purement réel 
allant de [1,+ce], fonctions intervenant sur des problèmes aux limites sur des hyperboloïdes à deux 
nappes (iso-surfaces du système sphéroïdal allongé ou prolat) 


Voir l'ouvrage « W.Magnus.F.Oberhettinger.R.PSoni-Formulas and Theorems for the Special 
Functions of Mathematical Physics », section 4.2, page 164. 


Sur l'intervalle [1,+], on ne va utiliser qu'une unique formule de liaison pour les fonctions 
coniques de Mehler, elle ne concerne plus la détermination des valeurs des fonctions pour des 
arguments inverses, puisque l'intervalle est uniquement positif. II s'agit donc plutôt de relier les 
fonctions de première et deuxième espèce. 


Pour x € [+ 1,+00| 
VueR VveC P“(x)=P# (x) 
Pour xe [+ 1,+00| 


-uri 








lu (x) = Cos(vr (v HR nr a v) (ot (x)- Qr (x)) 
v= 3 +iT => 0 = T (o 6) os] 


i Sine u +5) 4 soir) 


2 








= Py 0- E T 0-00) 


2 Sinha a+ sir ass tie) 
TT 
Sin (- 


Tv) (v+1)=- 7 sii ie (Z+ J= Z+) = 
Sin(vr) 2 2 2 “ir}r) Cosh(rx) 


De — fo, 9-01, (D) ED (01, 6) 


Sinner Z+ ir) 


2 




















NI! 








Les valeurs des fonctions coniques de première espèce sont réelles en appliquant la formule miroir 
et la première formule de liaison sur les fonctions de première espèce 


P#(x)= P” (x) i 7 s ne. 
P”(x)= PE P? lez 0-7 Pi Pi (= P’ Eo 


=P“ (x) à valeur réelle 


——+iT 
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La fonction conique de deuxième espèce quant à elle est à valeur imaginaire et sa partie 


imaginaire est proportionnelle à la fonction de première espèce : 
x Ef1,+00) 


jCotanh 
Pi, C x)= ICoIan refo MORIO) 
2 2 
Propriété miroir de Q 


A 9,06) be 9,0 = 01,0) 
(x )= comte | 9 : ()-01,5)-- 2Cotanh(tx) inf 0 3 o] 
= E mT Dr 





tir 


H „O= -F Tamer, 16 


Donc comme auparavant dans l'intervalle [-1,+1] et pour des arguments réels, le choix d'une 
fonction indépendante de deuxième espèce, se réalise par exemple à partir de la partie réelle de 
de la fonction conique de deuxième espèce. 


Certains auteurs comme par exemple dans l'article A.F.Ulitko « On a generalization of the integral 
transformation of Mehler-Fock », SOVIET APPLIED MECHANICS, Prikladnaya Mekhmnika, Vot. 3, 
No. 5, pp. 45-49, 1967», est introduit les deux fonctions de première et deuxième espèce : 
x=Cosh(n) n e€[0,+0) => xell,+) 


7,0 o, -0 O) 2,0 0, (he, o] 


Que l'on relie immédiatement aux parties réelles et imaginaires des fonctions coniques de 
deuxième espèce, et qui démontre que les deux fonctions introduites sont bien à valeur réelle. 


L So, o, o)Ra, (1) 
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Les Wronskiens de ces diverses solutions sont : 
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Développement en série et dérivées paramétriques de la fonction conique de Mehler de 
première espèce 


Le développement en série de la fonction de Legendre suivant est également utilisé pour obtenir 
celui des diverses dérivées paramétriques : 
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Ces développements en série sont valables autant dans l'intervalle [-1,+1], que sur l'intervalle 
[1,+ce]. Appliquées ces formules sur les fonctions coniques de Mehler, il vient : 
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Développement en série et dérivées paramétriques de la fonction conique de Mehler de 


deuxième espèce 


Pour les fonctions Q de Legendre, on développe les dérivées premières du paramètre réel v selon 
les deux formules suivantes : 
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Cela donne pour les fonctions coniques, les deux formules suivantes : 
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